Abstract: Structural symmetry breaking in two-dimensional materials can lead to superior physical properties and introduce an additional degree of piezoelectricity. In the present paper, we propose three structural phases (1H, 1T, and 1T') of Janus WXO (X = S, Se, and Te) monolayers and investigate their vibrational, thermal, elastic, piezoelectric, and electronic properties by using first-principles methods. Phonon spectra analysis reveals that while the 1H phase is dynamically stable, the 1T phase exhibits imaginary frequencies and transforms to the distorted 1T' phase. Ab initio molecular dynamics simulations confirm that 1H- and 1T'-WXO monolayers are thermally stable even at high temperatures without any significant structural deformations. Different from binary systems, additional Raman active modes appear upon the formation of Janus monolayers. Although the mechanical properties of 1H-WXO are found to be isotropic, they are orientation dependent for 1T'-WXO. It is also shown that 1H-WXO monolayers are indirect band-gap semiconductors and the band gap narrows down the chalcogen group. Except 1T'-WSO, 1T'-WXO monolayers have a narrow band gap correlated with the Peierls distortion. The effect of spin-orbit coupling on the band structure is also examined for both phases and the alteration in the band gap is estimated. The versatile mechanical and electronic properties of Janus WXO monolayers together with their large piezoelectric response imply that these systems are interesting for several nanoelectronic applications.

Abstract: By using density functional theory (DFT) based first-principles calculations, electronic, vibrational, piezo-electric, and optical properties of monoclinic Janus single-layer Ga2TeX (X = S or Se) are investigated. The dynamical, mechanical, and thermal stability of the proposed Janus single layers are verified by means of phonon bands, stiffness tensor, and quantum molecular dynamics simulations. The calculated vibrational spectrum reveals the either pure or coupled optical phonon branches arising from Ga-Te and Ga-X atoms. In addition to the in-plane anisotropy, single-layer Janus Ga2TeX exhibits additional out-of-plane asymmetry, which leads to important consequences for its electronic and optical properties. Electronic band dispersions indicate the direct band-gap semiconducting nature of the constructed Janus structures with energy band gaps falling into visible spectrum. Moreover, while orientation-dependent linear-elastic properties of Janus single layers indicate their strong anisotropy, the calculated in-plane stiffness values reveal the ultrasoft nature of the structures. In addition, predicted piezoelectric coefficients show that while there is a strong in-plane anisotropy between piezoelectric constants along armchair (AC) and zigzag (ZZ) directions, there exists a tiny polarization along the out-of-plane direction as a result of the formation of Janus structure. The optical response to electromagnetic radiation has been also analyzed through density functional theory by considering the independent-particle approximation. Finally, the optical spectra of Janus Ga2TeX structures is investigated and it showed a shift from the ultraviolet region to the visible region. The fact that the spectrum is between these regions will allow it to be used in solar energy and many nanoelectronics applications. The predicted monoclinic single-layer Janus Ga2TeX are relevant for promising applications in optoelectronics, optical dichroism, and anisotropic nanoelasticity.

Abstract: The effect of the spatially varying spin current on a skyrmion is numerically investigated. It is shown that an inhomogeneous current density induces an elongation of the skyrmion. This elongation can be controlled using current pulses of different strength and duration. Long current pulses lead to a splitting that forms two replicas of the initial skyrmion while for short pulses the elongated skyrmion relaxes back to its initial circular state through rotation in the MHz-GHz frequency range. The frequency is dependent on the strength of the damping coefficient.

Abstract: Thermal expansion behavior of two-dimensional (2D) nitrides and graphene were studied by ab initio molecular dynamics (MD) simulations as well as quasiharmonic approximation (QHA). Anharmonicity of the acoustic phonon modes are related to the unusual negative thermal expansion (NTE) behavior of the nitrides. Our results also hint that direct ab initio MD simulations are a more elaborate method to investigate thermal expansion behavior of 2D materials than the QHA. Nevertheless, giant NTE coefficients are found for h-GaN and h-AlN within the covered temperature range 100-600 K regardless of the chosen computational method. This unusual NTE of 2D nitrides is reasoned with the out-of-plane oscillations related to the rippling behavior of the monolayers.

Abstract: Different types of titanosilicates are synthesized, structurally characterized, and subsequently catalytically tested in the liquid-phase epoxidation of cyclohexene. The performance of three types of combined zeolitic/mesoporous materials is compared with that of widely studied Ti-grafted-MCM-41 molecular sieve and the TS-1 microporous titanosilicate. The catalytic test results are correlated with the structural characteristics of the different catalysts. Moreover, for the first time, immersion calorimetry with the same substrate molecule as in the catalytic test reaction is applied as an extra means to interpret the catalytic results. A good correlation between catalytic performance and immersion calorimetry results is found. This work points out that the combination of catalytic testing and immersion calorimetry can lead to important insights into the influence of the materials structural characteristics on catalysis. Moreover, the potential of using immersion calorimetry as a screening tool for catalysts in epoxidation reactions is shown.

Abstract: The self-assembly of nanoparticles into clusters and the effect of the different parameters of the competing interaction potential on it are investigated. For a small number of particles, the structural organization of the clusters is almost unaffected by the attractive part of the potential, and for an intermediate number of particles the configuration strongly depends on the strength of it. The cluster size is controlled by the range of the interaction potential, and the structural arrangement is guided by the strength of the potential: i.e., the self-assembled cluster transforms from a faceted configuration at low strength to a spherical shell-like structure at high strength. Nonmonotonic behavior of the cluster size is found by increasing the interaction range. An approximate analytical expression is obtained that predicts the smallest cluster for a specific set of potential parameters. A Mendeleev-like table is constructed for different values of the strength and range of the attractive part of the potential in order to understand the structural ordering of the ground-state configuration of the self-assembled clusters.

Abstract: In vitro experiments provide a solid basis for understanding the interactions between particles and biological systems. An important confounding variable for these studies is the difference between the amount of particles administered and that which reaches the surface of cells. Here, we engineer a hydrogel-based nanoparticle system and combine in situ characterization techniques, 3D-printed cell cultures, and computational modeling to evaluate and study particle cell interactions of advanced particle systems. The framework presented demonstrates how sedimentation and diffusion can explain differences in particle cell association, and provides a means to account for these effects. Finally, using in silico modeling, we predict the proportion of particles that reaches the cell surface using common experimental conditions for a wide range of inorganic and organic micro- and nanoparticles. This work can assist in the understanding and control of sedimentation and diffusion when investigating cellular interactions of engineered particles.

Abstract: We present a zero-dimensional plasma kinetics model, including both surface and gas phase kinetics, to determine the role of vibrationally excited states in plasma-catalytic ammonia synthesis. We defined a new method to systematically capture the conditions of dielectric barrier discharges (DBDs), including those found in packed bed DBDs. We included the spatial and temporal nature of such discharges by special consideration of the number of micro-discharges in the model. We introduce a parameter that assigns only a part of the plasma power to the microdischarges, to scale the model conditions from filamentary to uniform plasma. Because of the spatial and temporal behaviour of the micro-discharges, not all micro-discharges occurring in the plasma reactor during a certain gas residence time are affecting the molecules. The fraction of power considered in the model ranges from 0.005 %, for filamentary plasma, to 100 %, for uniform plasma. If vibrational excitation is included in the plasma chemistry, these different conditions, however, yield an ammonia density that is only varying within one order of magnitude. At only 0.05 % of the power put into the uniform plasma component, a model neglecting vibrational excitation clearly does not result in adequate amounts of ammonia. Thus, our new model, which accounts for the concept in which not all the power is deposited by the micro-discharges, but some part may also be distributed in between them, suggests that vibrational kinetic processes are really important in (packed bed) DBDs. Indeed, vibrational excitation takes place in both the uniform plasma between the micro-discharges and in the strong micro-discharges, and is responsible for an increased N2 dissociation rate. This is shown here for plasma-catalytic ammonia synthesis, but might also be valid for other gas conversion applications.

Abstract: Plasma has been studied for several years to convert CO2 into value-added products. If CO2 could be converted in the presence of H2O as a cheap H-source for making syngas and oxygenates, it would mimic natural photosynthesis. However, CO2/H2O plasmas have not yet been extensively studied, not by experiments, and certainly not computationally. Therefore, we present here a kinetic modelling study to obtain a greater understanding of the vibrational kinetics of a CO2/H2O microwave plasma. For this purpose, we first created an electron impact cross section set for H2O, using a swarm-derived method. We added the new cross section set and CO2/H2O-related chemistry to a pure CO2 model. While it was expected that H2O addition mainly causes quenching of the CO2 asymmetric mode vibrational levels due to the additional CO2/H2O vibrational-translational relaxation, our model shows that the modifications in the vibrational kinetics are mainly induced by the strong electron dissociative attachment to H2O molecules, causing a reduction in electron density, and the corresponding changes in the input of energy into the CO2 vibrational levels by electron impact processes.

Abstract: Nanosecond pulsed discharge plasma shows a high degree of non-equilibrium, and exhibits relatively high conversions in the dry reforming of methane. To further improve the application, a good insight of the underlying mechanisms is desired. We developed a chemical kinetics model to explore the underlying plasma chemistry in nanosecond pulsed discharge. We compared the calculated conversions and product selectivities with experimental results, and found reasonable agreement in a wide range of specific energy input. Hence, the chemical kinetics model is able to provide insight in the underlying plasma chemistry. The modeling results predict that the most important dissociation reaction of CO₂and CH₄is electron impact dissociation. C₂H₂is the most abundant hydrocarbon product, and it is mainly formed upon reaction of two CH₂radicals. Furthermore, the vibrational excitation levels of CO₂contribute for 85% to the total dissociation of CO₂.

Abstract: The use of microfluidic devices in the field of plasma-liquid interaction can unlock unique possibilities to investigate the effects of plasma-generated reactive species for environmental and biomedical applications. So far, very little simulation work has been performed on microfluidic devices in contact with a plasma source. We report on the modelling and computational simulation of physical and chemical processes taking place in a novel plasma-microfluidic platform. The main production and transport pathways of reactive species both in plasma and liquid are modelled by a novel modelling approach that combines 0D chemical kinetics and 2D transport mechanisms. This combined approach, applicable to systems where the transport of chemical species occurs in unidirectional flows at high Péclet numbers, decreases calculation times considerably compared to regular 2D simulations. It takes advantage of the low computational time of the 0D reaction models while providing spatial information through multiple plug-flow simulations to yield a quasi-2D model. The gas and liquid flow profiles are simulated entirely in 2D, together with the chemical reactions and transport of key chemical species. The model correctly predicts increased transport of hydrogen peroxide into the liquid when the microfluidic opening is placed inside the plasma effluent region, as opposed to inside the plasma region itself. Furthermore, the modelled hydrogen peroxide production and transport in the microfluidic liquid differs by less than 50% compared with experimental results. To explain this discrepancy, the limits of the 0D–2D combined approach are discussed.

Abstract: Groundwater fluxes in tidal marshes largely control key ecosystem functions and services, such as vegetation growth, soil carbon sequestration, and nutrient cycling. In tidal marshes restored on formerly embanked agricultural land, groundwater fluxes are often limited as compared to nearby natural marshes, as a result of historical agricultural soil compaction. To improve the functioning of restored tidal marshes, knowledge is needed on how much certain design options can optimize soil-groundwater interactions in future restoration projects. Based on measured data on soil properties and tidally induced groundwater dynamics, we calibrated and evaluated a 2D vertical model of a creek-marsh cross-section, accounting for both saturated and unsaturated groundwater flow and solute transport in a variably saturated groundwater flow model. We found that model simulations of common restoration practices such as soil amendments (increasing the depth of porous soil on top of the compact layer) and creek excavation (increasing the creek density) increase the soil aeration depth and time, the drainage depth and the solute flux, and decrease the residence time of solutes in the porewater. Our simulations indicate that increasing the depth to the compact layer from 20 cm to 40 cm, or increasing the creek density from 1 creek to 2 creeks along a 50 m marsh transect (while maintaining the total creek cross-sectional area), in both cases more than doubles the volume of water processed by the marsh soil. We discuss that this may stimulate nutrient cycling. As such, our study demonstrates that groundwater modelling can support the design of marsh restoration measures aiming to optimize groundwater fluxes and related ecosystem services.

Abstract: We investigate discharge contraction in a microwave plasma at sub-atmospheric pressure, operating in CO2 and CO2/CH4 mixtures. The rise of the electron number density with plasma contraction intensifies the gas heating in the core of the plasma. This, in turn, initiates fast core-periphery transport and defines the rate of thermal chemistry over plasma chemistry. In this context, power concentration describes the overall mechanism including plasma contraction and chemical kinetics. In a complex chemistry such as dry reforming of methane, transport of reactive species is essential to define the performance of the reactor and achieve the desired outputs. Thus, we couple experimental observations and thermodynamic calculations for model validation and understanding of reactor performance. Adding CH4 alters the thermodynamic properties of the mixture, especially the reactive component of the heat conductivity. The increase in reactive heat conductivity increases the pressure at which plasma contraction occurs, because higher rates of gas heating are required to reach the same temperature. In addition, we suggest that the predominance of heat conduction over convection is a key condition to observe the effect of heat conductivity on gas temperature.

Abstract: Since the onset of pattern transfer technologies for chip manufacturing, various strategies have been developed to circumvent or overcome aspect ratio dependent etching (ARDE). These methods have, however, their own limitations in terms of etch non-idealities, throughput or costs. Moreover, they have mainly been optimized for individual in-device features and die-scale patterns, while occasionally ending up with poor patterning of metrology marks, affecting the alignment and overlay in lithography. Obtaining a better understanding of the underlying mechanisms of ARDE and how to mitigate them therefore remains a relevant challenge to date, for both marks and advanced nodes. In this work, we accordingly assessed the neutral transport mechanisms in ARDE by means of experiments and multiscale modeling for SiO₂etching with CHF₃/Ar and CF₄/Ar plasmas. The experiments revealed a local maximum in the etch rate for an aspect ratio around unity, i.e. the simultaneous occurrence of regular and inverse reactive ion etching lag for a given etch condition. We were able to reproduce this ARDE trend in the simulations without taking into account charging effects and the polymer layer thickness, suggesting shadowing and diffuse reflection of neutrals as the primary underlying mechanisms. Subsequently, we explored four methods with the simulations to regulate ARDE, by varying the incident plasma species fluxes, the amount of polymer deposition, the ion energy and angular distribution and the initial hardmask sidewall angle, for which the latter was found to be promising in particular. Although our study focusses on feature dimensions characteristic to metrology marks and back-end-of-the-line integration, the obtained insights have a broader relevance, e.g. to the patterning of advanced nodes. Additionally, this work supports the insight that physisorption may be more important in plasma etching at room temperature than originally thought, in line with other recent studies, a topic on which we recommend further research.