Abstract: CsPbI3 is a promising perovskite material for photovoltaic applications in its photoactive perovskite or black phase. However, the material degrades to a photovoltaically inactive or yellow phase at room temperature. Various mitigation strategies are currently being developed to increase the lifetime of the black phase, many of which rely on inducing strains in the material that hinder the black-to-yellow phase transition. Physical insight into how these strategies exactly induce strain as well as knowledge of the spatial extent over which these strains impact the material is crucial to optimize these approaches but is still lacking. Herein, we combine machine learning potential-based molecular dynamics simulations with our in silico strain engineering approach to accurately quantify strained large-scale atomic structures on a nanosecond time scale. To this end, we first model the strain fields introduced by atomic substitutions as they form the most elementary strain sources. We demonstrate that the magnitude of the induced strain fields decays exponentially with the distance from the strain source, following a decay rate that is largely independent of the specific substitution. Second, we show that the total strain field induced by multiple strain sources can be predicted to an excellent approximation by summing the strain fields of each individual source. Finally, through a case study, we illustrate how this additive character allows us to explain how complex strain fields, induced by spatially extended strain sources, can be predicted by adequately combining the strain fields caused by local strain sources. Hence, the strain additivity proposed here can be adopted to further our insight into the complex strain behavior in perovskites and to design strain from the atomic level onward to enhance their sought-after phase stability.

Abstract: The present study focuses on the role of mass removal mechanisms in ns-laser ablation. A copper sample is placed in argon, initially set at standard pressure and temperature. Calculations are performed for a 6 ns laser pulse with a wavelength of 532 nm and laser fluences up to 10 J/cm2. The transient behavior in and above the copper target is described by a hydrodynamic model. Transmission profiles and ablation depths are compared with experimental results and similar trends are found. Our calculations reveal an interesting self-inhibiting mechanism: volumetric mass removal in the supercritical region triggers plasma shielding and therefore stops proceeding. This self-limiting process indicates that volumetric mass removal does not necessarily result in large ablation depths.

Abstract: In this paper, the splitting of CO2 in a pulsed plasma system, such as a dielectric barrier discharge (DBD), is evaluated from a chemical point of view by means of numerical modeling. For this purpose, a chemical reaction set of CO2 in an atmospheric pressure plasma is developed, including the vibrational states of CO2, O2, and CO. The simulated pulses are matched to the conditions of a filament (or microdischarge) and repeated with intervals of 1 μs. The influence of vibrationally excited CO2 as well as other neutral species, ions, and electrons on the CO2 splitting is discussed. Our calculations predict that the electrons have the largest contribution to the CO2 splitting at the conditions under study, by electron impact dissociation. The contribution of vibrationally excited CO2 levels in the splitting of CO2 is found be 6.4%, when only considering one microdischarge pulse and its afterglow, but it can be much higher for consecutive discharge pulses, as is typical for a filamentary DBD, when the interpulse time is short enough and accumulation effects in the vibrationally excited CO2 densities can occur.

Abstract: van der Waals (vdW) forces play an important role in the adsorption of molecules on the surface of solids. However, the choice of the most suitable vdW functional for different systems is an essential problem which must be addressed for different systems. The lack of a systematic study on the performance of the vdW functionals in the bulk and adsorption properties of metal-oxides motivated us to examine different vdW approaches and compute the bulk and molecular adsorption properties of alpha-Cr2O3, alpha-Fe2O3, and alpha-Al2O3. For the bulk properties, we compared our results for the heat of formation, cohesive energy, lattice parameters and bond distances between the different vdW functionals and available experimental data. Next we studied the adsorption of benzene and CH3 molecules on top of different oxide surfaces. We employed different approximations to exchange and correlation within DFT, namely, the Perdew-Burke-Ernzerhof (PBE) GGA, (PBE)+U, and vdW density functionals [ DFT(vdW-DF/DF2/optPBE/optB86b/optB88)+U] as well as DFT-D2/D3(+U) methods of Grimme for the bulk calculations and optB86b-vdW(+U) and DFT-D2(+U) for the adsorption energy calculations. Our results highlight the importance of vdW interactions not only in the adsorption of molecules, but importantly also for the bulk properties. Although the vdW contribution in the adsorption of CH3 (as a chemisorption interaction) is less important compared to the adsorption of benzene (as a physisorption interaction), this contribution is not negligible. Also adsorption of benzene on ferryl/chromyl terminated surfaces shows an important chemisorption contribution in which the vdW interactions become less significant.

Abstract: The microstructural complexity of Li-rich cathode materials has so far hampered understanding the critical link between size, morphology and structural defects with both capacity and voltage fadings that this family of materials exhibits. Li2MnO3 is used here as a model material to extract reliable structure–property
relationships that can be further exploited for the development of high-performing and long-lasting Li-rich oxides. A series of samples with microstructural variability have been prepared and thoroughly characterized using the FAULTS software, which allows quantification of planar defects and extraction of
average crystallite sizes. Together with transmission electron microscopy (TEM) and density functional theory (DFT) results, the successful application of FAULTS analysis to Li2MnO3 has allowed rationalizing the synthesis conditions and identifying the individual impact of concurrent microstructural features on
both voltage and capacity fadings, a necessary step for the development of high-capacity Li-ion cathode materials with enhanced cycle life.

Abstract: The physical processes in an Ar/N2 magnetron discharge used for the reactive sputter deposition of TiNx thin films were simulated with a 2d3v particle-in-cell/Monte Carlo collisions (PIC/MCC) model. Cathode currents and voltages were calculated self-consistently and compared with experiments. Also, ion fractions were calculated and validated with mass spectrometric measurements. With this PIC/MCC model, the influence of N2/Ar gas ratio on the particle densities and fluxes was investigated, taking into account the effect of the poisoned state of the target.

Abstract: We propose a simple way of probing the number of modes contributing to the channeling in graphene waveguides which are formed by a gauge potential produced by mechanical strain. The energy mode structure for both homogeneous and non-homogeneous strain regimes is carefully studied using the continuum description of the Dirac equation. We found that high strain values privilege negative (instead of positive) group velocities throughout the guidance, sorting the types of modes flowing through it. We also show how the effect of a substrate-induced gap competes against the strain.

Abstract: This work explains the need for plasma models, introduces arguments for choosing the type of model that better fits the purpose of each study, and presents the basics of the most common nonequilibrium low-temperature plasma models and the information available from each one, along with an extensive list of references for complementary in-depth reading. The paper presents the following models, organised according to the level of multi-dimensional description of the plasma: kinetic models, based on either a statistical particle-in-cell/Monte-Carlo approach or the solution to the Boltzmann equation (in the latter case, special focus is given to the description of the electron kinetics); multi-fluid models, based on the solution to the hydrodynamic equations; global (spatially-average) models, based on the solution to the particle and energy rate-balance equations for the main plasma species, usually including a very complete reaction chemistry; mesoscopic models for plasma–surface interaction, adopting either a deterministic approach or a stochastic dynamical Monte-Carlo approach. For each plasma model, the paper puts forward the physics context, introduces the fundamental equations, presents advantages and limitations, also from a numerical perspective, and illustrates its application with some examples. Whenever pertinent, the interconnection between models is also discussed, in view of multi-scale hybrid approaches.

Abstract: Plasma-catalytic ammonia synthesis is receiving ever increasing attention, especially in packed bed dielectric barrier discharge (DBD) reactors. The latter typically operate in the filamentary regime when used for gas conversion applications. While DBDs are in principle well understood and already applied in the industry, the incorporation of packing materials and catalytic surfaces considerably adds to the complexity of the plasma physics and chemistry governing the ammonia formation. We employ a plasma kinetics model to gain insights into the ammonia formation mechanisms, paying special attention to the role of filamentary microdischarges and their afterglows. During the microdischarges, the synthesized ammonia is actually decomposed, but the radicals created upon electron impact dissociation of N2 and H2 and the subsequent catalytic reactions cause a net ammonia gain in the afterglows of the microdischarges. Under our plasma conditions, electron impact dissociation of N2 in the gas phase followed by the adsorption of N atoms is identified as a rate-limiting step, instead of dissociative adsorption of N2 on the catalyst surface. Both elementary Eley−Rideal and Langmuir−Hinshelwood reaction steps can be found important in plasma-catalytic NH3 synthesis.

Abstract: The dielectric properties of confined water is of fundamental interest and is still controversial. For water confined in channels with height smaller than h = 8 angstrom, we found a commensurability effect and an extraordinary decrease in the out-of-plane dielectric constant down to the limit of the dielectric constant of optical water. Spatial resolved polarization density data obtained from molecular dynamics simulations are found to be antisymmetric across the channel and are used as input in a mean-field model for the dielectric constant as a function of the height of the channel for h > 15 angstrom. Our results are in excellent agreement with a recent experiment [L. Fumagalli et al., Science 360, 1339 (2018)].

Abstract: Optical-photothermal infrared (O-PTIR) spectroscopy is a recently developed technique that provides spectra comparable to traditional transmission FTIR spectroscopy with nanometric spatial resolution. Hence, O-PTIR is a promising candidate for the analysis of historical paintings, as well as other cultural heritage objects, but its potential has not yet been evaluated.

Abstract: Ultrathin metal nanoparticles coatings, synthesized by gas-phase deposition, are emerging as go-to materials in a variety of fields ranging from pathogens control, sensing to energy storage. Predicting their morphology and mechanical properties beyond a trial-and-error approach is a crucial issue limiting their exploitation in real-life applications. The morphology and mechanical properties of Ag nanoparticles ultrathin films, synthesized by supersonic cluster beam deposition, are here assessed adopting a bottom-up, multi-technique approach. A virtual film model is proposed merging high resolution scanning transmission electron microscopy, supersonic cluster beam dynamics and molecular dynamics simulations. The model is validated against mechanical nanometrology measurements and is readily extendable to metals other than Ag. The virtual film is shown to be a flexible and reliable predictive tool to access morphology-dependent properties such as mesoscale gas-dynamics and elasticity of ultrathin films synthesized by gas-phase deposition.

Abstract: A one-dimensional fluid model for a dielectric barrier discharge in CH4/O2 and CH4/CO2 gas mixtures is developed. The model describes the gas-phase chemistry for partial oxidation and for dry reforming of methane. The spatially averaged densities of the various plasma species are presented as a function of time and initial gas mixing ratio. Besides, the conversion of the inlet gases and the selectivities of the reaction products are calculated. Syngas, higher hydrocarbons, and higher oxygenates are typically found to be important reaction products. Furthermore, the main underlying reaction pathways for the formation of syngas, methanol, formaldehyde, and other higher oxygenates are determined.

Abstract: We have performed a first-principles study of the p- and n-type conductivity in CuIn1−xGaxSe2 due to native point defects, based on the HSE06 hybrid functional. Band alignment shows that the band gap becomes larger with x due to the increasing conduction band minimum, rendering it hard to establish n-type conductivity in CuGaSe2. From the defect formation energies, we find that In/GaCu is a shallow donor, while VCu, VIn/Ga and CuIn/Ga act as shallow acceptors. Using the total charge neutrality of ionized defects and intrinsic charge carriers to determine the Fermi level, we show that under In-rich growth conditions InCu causes strongly n-type conductivity in CuInSe2. Under increasingly In-poor growth conditions, the conductivity type in CuInSe2 alters to p-type and compensation of the acceptors by InCu reduces, as also observed in photoluminescence experiments. In CuGaSe2, the native acceptors pin the Fermi level far away from the conduction band minimum, thus inhibiting n-type conductivity. On the other hand, CuGaSe2 shows strong p-type conductivity under a wide range of Ga-poor growth conditions. Maximal p-type conductivity in CuIn1−xGaxSe2 is reached under In/Ga-poor growth conditions, in agreement with charge concentration measurements on samples with In/Ga-poor stoichiometry, and is primarily due to the dominant acceptor CuIn/Ga.

Abstract: Co-rich ZnCoO nanoparticles embedded in wurtzite Zn0.7Co0.3O thin films are grown by pulsed laser deposition on a Si substrate. Local superconductivity with an onset Tc at 5.9 K is demonstrated in the hybrid system. The unexpected superconductivity probably results from Co(3+) in the Co-rich ZnCoO nanoparticles or from the interface between the Co-rich nanoparticles and the Zn0.7Co0.3O matrix.