Abstract: We calculate bubble nucleation rates in a Lennard-Jones fluid through explicit molecular dynamics simulations. Our approach-based on a recent free energy method (dubbed reweighted Jarzynski sampling), transition state theory, and a simple recrossing correction-allows us to probe a fairly wide range of rates in several superheated and cavitation regimes in a consistent manner. Rate predictions from this approach bridge disparate independent literature studies on the same model system. As such, we find that rate predictions based on classical nucleation theory, direct brute force molecular dynamics simulations, and seeding are consistent with our approach and one another. Published rates derived from forward flux sampling simulations are, however, found to be outliers. This study serves two purposes: First, we validate the reliability of common modeling techniques and extrapolation approaches on a paradigmatic problem in materials science and chemical physics. Second, we further test our highly generic recipe for rate calculations, and establish its applicability to nucleation processes.

Abstract: Ptychography provides highly efficient imaging in scanning transmission electron microscopy (STEM), but questions have remained over its applicability to strongly scattering samples such as those most commonly seen in materials science. Although contrast reversals can appear in ptychographic phase images as the projected potentials of the sample increase, we show here how these can be easily overcome by a small amount of defocus. The amount of defocus is small enough that it not only can exist naturally when focusing using the annular dark field (ADF) signal but can also be adjusted post acquisition. The ptychographic images of strongly scattering materials are clearer at finite doses than other STEM techniques and can better reveal light atomic columns within heavy lattices. In addition, data for ptychography can now be collected simultaneously with the fastest of ADF scans. This combination of sensitivity and interpretability presents an ideal workflow for materials science.

Abstract: The evaluation of nucleation rates from molecular dynamics trajectories is hampered by the slow nucleation time scale and impact of finite size effects. Here, we show that accurate nucleation rates can be obtained in a very general fashion relying only on the free energy barrier, transition state theory, and a simple dynamical correction for diffusive recrossing. In this setup, the time scale problem is overcome by using enhanced sampling methods, in casu metadynamics, whereas the impact of finite size effects can be naturally circumvented by reconstructing the free energy surface from an appropriate ensemble. Approximations from classical nucleation theory are avoided. We demonstrate the accuracy of the approach by calculating macroscopic rates of droplet nucleation from argon vapor, spanning 16 orders of magnitude and in excellent agreement with literature results, all from simulations of very small (512 atom) systems.

Abstract: We present here a new interatomic potential parameter set to predict the thermal conductivity of zirconium trisulfide monolayers. The generated Tersoff-type force field is parameterized using data collected with first-principles calculations. We use non-equilibrium molecular dynamics simulations to predict the thermal conductivity. The generated parameters result in very good agreement in structural, mechanical, and dynamical parameters. The room temperature lattice thermal conductivity ( kappa) of the considered crystal is predicted to be kappa x x = 25.69Wm – 1K – 1 and kappa y y = 42.38Wm – 1K – 1, which both agree well with their corresponding first-principles values with a discrepancy of less than 5%. Moreover, the calculated kappa variation with temperature (200 and 400 K) are comparable within the framework of the accuracy of both first-principles and molecular dynamics simulations.

Abstract: Recently, a two-dimensional (2D) MoSi 2N 4 (MSN) structure has been successfully synthesized [Hong et al., Science 369(6504), 670-674 (2020)]. Motivated by this result, we investigate the structural, electronic, and optical properties of MSN monolayer (MSN-1L) and bilayer (MSN-2L) under the applied electric field (E-field) and strain using density functional theory calculations. We find that the MSN-2L is a semiconductor with an indirect bandgap of 1.60 (1.80)eV using Perdew-Burke-Ernzerhof (HSE06). The bandgap of MSN-2L decreases as the E-field increases from 0.1 to 0.6V/angstrom and for larger E-field up to 1.0V/angstrom the bilayer becomes metallic. As the vertical strain increases, the bandgap decreases; more interestingly, a semiconductor to a metal phase transition is observed at a strain of 12 %. Furthermore, the optical response of the MSN-2L is in the ultraviolet (UV) region of the electromagnetic spectrum. The absorption edge exhibits a blue shift by applying an E-field or a vertical compressive strain. The obtained interesting properties suggest MSN-2L as a promising material in electro-mechanical and UV opto-mechanical devices.

Abstract: Polarization fatigue, i.e., the loss of polarization of ferroelectric capacitors upon field cycling, has been widely discussed as an interface related effect. However, mechanism(s) behind the development of fatigue have not been fully identified. Here, we study the fatigue mechanisms in Pt/PbZr0.52Ti0.48O3/SrRuO3 (Pt/PZT/SRO) capacitors in which all layers are fabricated by pulsed laser deposition without breaking the vacuum. With scanning transmission electron microscopy, we observed that in the fatigued capacitor, the Pt/PZT interface becomes structurally degraded, forming a 5 nm-10 nm thick non-ferroelectric layer of crystalline ZrO2 and diffused Pt grains. We then found that the fatigued capacitors can regain the full initial polarization switching if the externally applied field is increased to at least 10 times the switching field of the pristine capacitor. These findings suggest that polarization fatigue is driven by a two-step mechanism. First, the transient depolarization field that repeatedly appears during the domain switching under field cycling causes decomposition of the metal/ferroelectric interface, resulting in a non-ferroelectric degraded layer. Second, this interfacial non-ferroelectric layer screens the external applied field causing an increase in the coercive field beyond the usually applied maximum field and consequently suppresses the polarization switching in the cycled capacitor. Our work clearly confirms the key role of the electrode/ferroelectric interface in the endurance of ferroelectric-based devices.

Abstract: The metal-insulator transition (MIT) in strongly correlated oxides is a topic of great interest for its potential applications, such as Mott field effect transistors and sensors. We report that the MIT in high quality epitaxial SrVO3 (SVO) thin films is present as the film thickness is reduced, lowering the dimensionality of the system, and electron-electron correlations start to become the dominant interactions. The critical thickness of 3 u.c is achieved by avoiding effects due to off-stoichiometry using optimal growth conditions and excluding any surface effects by a STO capping layer. Compared to the single SVO thin films, conductivity enhancement in SVO/STO superlattices is observed. This can be explained by the interlayer coupling effect between SVO sublayers in the superlattices. Magnetoresistance and Hall measurements indicate that the dominant driving force of MIT is the electron–electron interaction.

Abstract: Positron annihilation spectroscopy in the Doppler and coincidence Doppler mode was applied on Ge1 xSnx epitaxial layers, grown by chemical vapor deposition with different total As concentrations (1019-1021 cm3), high active As concentrations (1019 cm3), and similar Sn concentrations (5.9%-6.4%). Positron traps are identified as mono-vacancy complexes. Vacancy-As complexes, V-Asi, formed during the growth were studied to deepen the understanding of the electrical passivation of the Ge1 xSnx:As epilayers. Larger monovacancy complexes, V-Asi (i 2), are formed as the As doping increases. The total As concentration shows a significant impact on the saturation of the number of As atoms (i 1/4 4) around the vacancies in the sample epilayers. The presence of V-Asi complexes decreases the dopant activation in the Ge1 xSnx:As epilayers. Furthermore, the presence of Sn failed to hinder the formation of larger V-Asi complexes and thus failed to reduce the donor-deactivation.

Abstract: We examined the vacancy trapping proficiency of Sn and P atoms in germanium using positron annihilation spectroscopy measurements, sensitive to the open-volume defects. Epitaxial Ge1 xSnx films were grown by chemical vapor deposition with different P concentrations in the 3: 0 1019-1: 5 1020 cm 3 range. We corroborate our findings with first principles simulations. Codoping of Ge with a Sn concentration of up to 9% is not an efficient method to suppress the free vacancy concentration and the formation of larger phosphorus-vacancy complexes. Experimental results confirm an increase in the number of P atoms around the monovacancy with P-doping, leading to dopant deactivation in epitaxial germanium-tin layers with similar Sn content. Vice versa, no impact on the improvement of maximum achieved P activation in Ge with increasing Sn-doping has been observed. Theoretical calculations also confirm that Pn-V (vacancy) complexes are energetically more stable than the corresponding SnmPn-V and Snm-V defect structures with the same number of alien atoms (Sn or P) around the monovacancy. he strong attraction of vacancies to the phosphorus atoms remains the dominant dopant deactivation mechanism in Ge as well as in Ge1 xSnx. Published under license by AIP Publishing.

Abstract: Strain has a strong effect on the properties of materials and the performance of electronic devices. Their ever shrinking size translates into a constant demand for accurate and precise measurement methods with a very high spatial resolution. In this regard, transmission electron microscopes are key instruments thanks to their ability to map strain with a subnanometer resolution. Here, we present a method to measure strain at the nanometer scale based on the diffraction of electron Bessel beams. We demonstrate that our method offers a strain sensitivity better than 2.5 × 10−4 and an accuracy of 1.5 × 10−3, competing with, or outperforming, the best existing methods with a simple and easy to use experimental setup.

Abstract: Particle-in-cell/Monte Carlo collision (PIC/MCC) simulations are performed to investigate the asymmetric secondary electron emission (SEE) effects when electrons strike two different material electrodes in low pressure capacitively coupled plasmas (CCPs). To describe the electron-surface interactions, a realistic model, considering the primary electron impact energy and angle, as well as the corresponding surface property-dependent secondary electron yields, is employed in PIC/MCC simulations. In this model, three kinds of electrons emitted from the surface are considered: (i) elastically reflected electrons, (ii) inelastically backscattered electrons, and (iii) electron induced secondary electrons (SEs, i.e., delta-electrons). Here, we examined the effects of electron-surface interactions on the ionization dynamics and plasma characteristics of an argon discharge. The discharge is driven by a voltage source of 13.56MHz with amplitudes in the range of 200-2000V. The grounded electrode material is copper (Cu) for all cases, while the powered electrode material is either Cu or silicon dioxide (SiO2). The simulations reveal that the electron impact-induced SEE is an essential process at low pressures, especially at high voltages. Different electrode materials result in an asymmetric response of SEE. Depending on the instantaneous local sheath potential and the phase of the SEE, these SEs either are reflected by the opposite sheath or strike the electrode surface, where they can induce delta-electrons upon their residual energies. It is shown that highly energetic delta-electrons contribute significantly to the ionization rate and a self-bias forms when the powered electrode material is assumed to be made of SiO2. Complex dynamics is observed due to the multiple electron-surface interaction processes and asymmetric yields of SEs in CCPs.

Abstract: The study presents a 2D fluid plasma model of a gliding arc discharge for dissociation of CO2 which allows its subsequent conversion into value-added chemicals. The model is based on the balance equations of charged and neutral particles, the electron energy balance equation, the gas thermal balance equation and the current continuity equation. By choosing the modeling domain to be the plane perpendicular to the arc current, the numerical calculations are significantly simplified. Thus, the model allows us to explore the influence of the gas instabilities (turbulences) on the energy efficiency of CO2 conversion. This paper presents results for plasma parameters at different values of the effective turbulent thermal conductivity leading to enhanced energy transport.

Abstract: Nanoindentation tests performed on nanocrystalline palladium films subjected to hydriding/dehydriding cycles demonstrate a significant softening when compared to the as-received material. The origin of this softening is unraveled by combining in situ TEM nanomechanical testing with automated crystal orientation mapping in TEM and high resolution TEM. The softening is attributed to the presence of a high density of stacking faults and of Shockley partial dislocations after hydrogen loading. The hydrogen induced defects affect the elementary plasticity mechanisms and the mechanical response by acting as preferential sites for twinning/detwinning during deformation. These results are analyzed and compared to previous experimental and simulation works in the literature. This study provides new insights into the effect of hydrogen on the atomistic deformation and cracking mechanisms as well as on the mechanical properties of nanocrystalline thin films and membranes.

Abstract: Due to Klein tunneling, electrostatic confinement of electrons in graphene is not possible. This hinders the use of graphene for quantum dot applications. Only through quasi-bound states with finite lifetime has one achieved to confine charge carriers. Here, we propose that bilayer graphene with a local region of decoupled graphene layers is able to generate bound states under the application of an electrostatic gate. The discrete energy levels in such a quantum blister correspond to localized electron and hole states in the top and bottom layers. We find that this layer localization and the energy spectrum itself are tunable by a global electrostatic gate and that the latter also coincides with the electronic modes in a graphene disk. Curiously, states with energy close to the continuum exist primarily in the classically forbidden region outside the domain defining the blister. The results are robust against variations in size and shape of the blister which shows that it is a versatile system to achieve tunable electrostatic confinement in graphene. Published by AIP Publishing.