Abstract: Two-dimensional materials (2DMs) continue to attract a lot of attention, particularly for their extreme flexibility and superior thermal properties. Molecular dynamics simulations are among the most powerful methods for computing these properties, but their reliability depends on the accuracy of interatomic interactions. While first principles approaches provide the most accurate description of interatomic forces, they are computationally expensive. In contrast, classical force fields are computationally efficient, but have limited accuracy in interatomic force description. Machine learning interatomic potentials, such as Gaussian Approximation Potentials, trained on density functional theory (DFT) calculations offer a compromise by providing both accurate estimation and computational efficiency. In this work, we present a systematic procedure to develop Gaussian approximation potentials for selected 2DMs, graphene, buckled silicene, and h-XN (X = B, Al, and Ga, as binary compounds) structures. We validate our approach through calculations that require various levels of accuracy in interatomic interactions. The calculated phonon dispersion curves and lattice thermal conductivity, obtained through harmonic and anharmonic force constants (including fourth order) are in excellent agreement with DFT results. HIPHIVE calculations, in which the generated GAP potentials were used to compute higher-order force constants instead of DFT, demonstrated the first-principles level accuracy of the potentials for interatomic force description. Molecular dynamics simulations based on phonon density of states calculations, which agree closely with DFT-based calculations, also show the success of the generated potentials in high-temperature simulations.

Abstract: X-ray diffraction and energy dispersive x-ray spectroscopic analyses showed a natural galena (PbS) crystal from Freiberg in Saxony (Germany) to be a single phase specimen [rock salt (NaCl) structure type, space group Fm3m, a = 5.932(1) angstrom] with stoichiometric composition and an enhanced dislocation density (8 approximate to 1011 cm-2). The latter parameter leads to an increase of the electrical resistivity in the high-temperature regime, as well as to the appearance of phonon resonance with a characteristic frequency coPR = 3.8(1) THz. Being in the same range (i.e., 3-5.5 THz) with the sulfur optical modes of highest group velocities, it results in a drastic reduction (by similar to 75%) of thermal conductivity (K) at lower temperatures (i.e., < 100 K), as well as in the appearance of a characteristic minimum in K at T approximate to 30 K. Furthermore, the studied galena is characterized by phonon-drag behavior and by temperature dependent switch of the charge carrier scattering mechanism regime (i.e., scattering on dislocations for T < 100 K, on acoustic phonons for 100 K < T < 170 K and on both acoustic and optical phonons for 170 K < T < 300 K). The combined theoretical calculation and optical spectroscopic study confirm this mineral to be a direct gap degenerate semiconductor. The possible origins of the second-order Raman spectrum are discussed.

Abstract: We investigate theoretically the ballistic transport in a two-dimensional electron gas, which is rolled up as a tube and is micro-structured into a Hall bar. A uniform magnetic field applied to such a curved surface results in a non-uniform perpendicular magnetic field. The bend resistances become asymmetric with respect to the orientation of the magnetic field due to the varying magnetic field along the junction. The resistance asymmetry is strongly affected by corrugation due to the varying mobility along different crystallographic directions. We compare our results with a recent transport measurement.

Abstract: Environmental magnetism, and the magnetic leaf signal in particular, is amply investigated and applied as proxy for atmospheric particulate matter pollution. In this study, we investigated the magnetic signal of annual segments of tree branches, and the composition of particles deposited hereon. Branches are, contrary to leaves of deciduous trees, available during leaf-off seasons and exposed to air pollution year-round. We examined the intra- and inter-tree variation in saturation isothermal remanent magnetization (SIRM) of branch internodes of London plane (Platanus x aeerifolia Willd.) trees in an urban environment. The branch SIRM, normalized by surface area, ranged from 18 to 650 x 10(-6) angstrom; the median amounted to 106 x 10(-6) angstrom. Most of the branch magnetic signal was attributed to the epidermis or bark, and the presence of metal-containing particles on the branch surfaces was confirmed by SEM-EDX. The location of the trees and the height, the depth in the crown and the age of the branches significantly influenced the branch SIRM. The median branch SIRM was up to 135% higher near a busy ring road than in quiet environments (city park and quiet street canyon), and was linked to the presence of Fe-rich particles with co-occurrence of trace metals such as Cr, Cu, Zn and Mn on the branch surface. Within the tree crowns, the branch SIRM generally decreased with increasing height, and was 22% higher in the interior than at the periphery of the crowns. Within the branches, the SIRM increased with each year of exposure, but did not relate to year-to-year variation in particle concentrations due to branch surface changes (epidermis shedding). Our results provide indications that branches can be a valuable alternative for biomagnetic monitoring of particulate pollution, but intra-tree variability in branch SIRM can be substantial due to the branch's location in the tree and branch age.

Abstract: The electrochemical properties and phase transformations during (de)insertion of Li+ in LiFePO4, LiFe0.9Mn0.1PO4 and LiFe0.5Mn0.5PO4 are studied by means of galvanostatic cycling, potential intermittent titration technique (PITT) and in situ X-ray powder diffraction. Different modes of switching between the solid solution and two-phase regimes are revealed which are influenced by the Mn content in Li1-xFe1-yMnyPO4. Additionally, an increase in electrochemical capacity with the Mn content is observed at high rates of galvanostatic cycling (10C, 20C), which is in good agreement with the numerically estimated contribution of the solid solution mechanism determined from PITT data. The observed asymmetric behavior of the phase transformations in Li1-xFe0.5Mn0.5PO4 during charge and discharge is discussed. For the first time, the crystal structures of electrochemically deintercalated Li1-xFe0.5Mn0.5PO4 with different Li content – LiFe0.5Mn0.5PO4, Li0.5Fe0.5Mn0.5PO4 and Li0.1Fe0.5Mn0.5PO4 – are refined, including the occupancy factors of the Li position. This refinement is done using electron diffraction tomography data. The crystallographic analyses of Li1-xFe0.5Mn0.5PO4 reveal that at x = 0.5 and 0.9 the structure retains the Pnma symmetry and the main motif of the pristine x = 0 structure without noticeable short range order effects.

Abstract: We review the energy spectrum and transport properties of several types of one-dimensional superlattices (SLs) on single-layer and bilayer graphene. In single-layer graphene, for certain SL parameters an electron beam incident on an SL is highly collimated. On the other hand, there are extra Dirac points generated for other SL parameters. Using rectangular barriers allows us to find analytical expressions for the location of new Dirac points in the spectrum and for the renormalization of the electron velocities. The influence of these extra Dirac points on the conductivity is investigated. In the limit of δ-function barriers, the transmission T through and conductance G of a finite number of barriers as well as the energy spectra of SLs are periodic functions of the dimensionless strength P of the barriers, Graphic, with vF the Fermi velocity. For a KronigPenney SL with alternating sign of the height of the barriers, the Dirac point becomes a Dirac line for P = π/2+nπ with n an integer. In bilayer graphene, with an appropriate bias applied to the barriers and wells, we show that several new types of SLs are produced and two of them are similar to type I and type II semiconductor SLs. Similar to single-layer graphene SLs, extra Dirac points are found in bilayer graphene SLs. Non-ballistic transport is also considered.

Abstract: Radiation damage places a fundamental limitation on the ability of microscopy to resolve many types of materials at high resolution. Here we evaluate the dose efficiency of phase contrast imaging with electron ptychography. The method is found to be far more resilient to temporal incoherence than conventional and spherical aberration optimized phase contrast imaging, resulting in significantly greater clarity at a given dose. This robustness is explained by the presence of achromatic lines in the four dimensional ptychographic dataset.

Abstract: Electron microscopy is a powerful tool for studying the properties of materials down to their atomic structure. In many cases, the quantitative interpretation of images requires simulations based on atomistic structure models. These typically use the independent atom approximation that neglects bonding effects, which may, however, be measurable and of physical interest. Since all electrons and the nuclear cores contribute to the scattering potential, simulations that go beyond this approximation have relied on computationally highly demanding all-electron calculations. Here, we describe a new method to generate ab initio electrostatic potentials when describing the core electrons by projector functions. Combined with an interface to quantitative image simulations, this implementation enables an easy and fast means to model electron scattering. We compare simulated transmission electron microscopy images and diffraction patterns to experimental data, showing an accuracy equivalent to earlier all-electron calculations at a much lower computational cost.

Abstract: We have studied the temperature and field dependencies of the critical current I(C) in the Nb-Fe(0.1)Si(0.9)-Nb Josephson junction with a tunneling barrier formed by a paramagnetic insulator. We demonstrate that in these junctions coexistence of both the 0 and the pi states within one tunnel junction occurs, and leads to the appearance of a sharp cusp in the temperature dependence I(C)(T), similar to the I(C)(T) cusp found for the 0-pi transition in metallic pi junctions. This cusp is not related to the 0-pi temperature-induced transition itself, but is caused by the different temperature dependencies of the opposing 0 and pi supercurrents through the barrier.

Abstract: We demonstrate for the first time a possibility to vary the anion content in perovskites over a wide range through a long-range-ordered arrangement of crystallographic shear (CS) planes. Anion-deficient perovskites (Pb,Bi)1−xFe1+xO3−y with incommensurately modulated structures were prepared as single phases in the compositional range from Pb0.857Bi0.094Fe1.049O2.572 to Pb0.409Bi0.567Fe1.025O2.796. Using a combination of electron diffraction and high-resolution scanning transmission electron microscopy, we constructed a superspace model describing a periodic arrangement of the CS planes. The model was verified by refinement of the Pb0.64Bi0.32Fe1.04O2.675 crystal structure from neutron powder diffraction data ((3 + 1)D S.G. X2/m(α0γ), X = [1/2,1/2,1/2,1/2], a = 3.9082(1) Å, b = 3.90333(8) Å, c = 4.0900(1) Å, β = 91.936(2)°, q = 0.05013(4)a* + 0.09170(3)c* at T = 700 K, RP = 0.036, RwP = 0.048). The (Pb,Bi)1−xFe1+xO3−y structures consist of perovskite blocks separated by CS planes confined to nearly the (509)p perovskite plane. Along the CS planes, the perovskite blocks are shifted with respect to each other over the 1/2[110]p vector that transforms the corner-sharing connectivity of the FeO6 octahedra in the perovskite framework to an edge-sharing connectivity of the FeO5 pyramids at the CS plane, thus reducing the oxygen content. Variation of the chemical composition in the (Pb,Bi)1−xFe1+xO3−y series occurs mainly because of a changing thickness of the perovskite block between the interfaces, that can be expressed through the components of the q vector as Pb6γ+2αBi1−7γ−αFe1+γ−αO3−3γ−α. The Pb, Bi, and Fe atoms are subjected to strong displacements occurring in antiparallel directions on both sides of the perovskite blocks, resulting in an antiferroelectric-type structure. This is corroborated by the temperature-, frequency-, and field-dependent complex permittivity measurements. Pb0.64Bi0.32Fe1.04O2.675 demonstrates a remarkably high resistivity >0.1 T Ω cm at room temperature and orders antiferromagnetically below TN = 608(10) K.

Abstract: Context: Asthma affects 20 million Americans resulting in an economic burden of approximately $18 billion in the US alone (Allergies and Asthma Foundation 2000; National Center for Environmental Health (NCEH) 1999). Research studies based on differences in patient-specific airway morphology for asthma and the associated effect on deposition of inhaled aerosols are currently not available in the literature. Therefore, the role of morphological variations such as upper airway (extrathoracic) occlusion is not well documented. Objective: Functional imaging based computational fluid dynamics (CFD) of the respiratory airways for five asthmatic subjects is performed in this study using computed tomography (CT) based patient-specific airway models and boundary conditions. Methods: CT scans for 5 asthma patients were used to reconstruct 3D lung models using segmentation software. An averaged inhalation profile and patient-specific lobar flow distribution were used to perform the simulation. The simulations were used to obtain deposition for BDP/Formoterol (R) HFA pMDI in the patient-specific airway models. Results: The lung deposition obtained using CFD was in excellent agreement with available in vivo data using the same product. Specifically, CFD resulted in 30% lung deposition, whereas in vivo lung deposition was reported to be approximately 31%. Conclusion: It was concluded that a combination of patient-specific airway models and lobar boundary conditions can be used to obtain accurate lung deposition estimates. Lower lung deposition can be expected for patients with higher extrathoracic resistance. Novel respiratory drug delivery devices need to accommodate population subgroups based on these morphological and anatomical differences in addition to subject age.

Abstract: The (AO)(ABO3)n Ruddlesden-Popper structure is an archetypal complex oxide consisting of two distinct structural units, an (AO) rock salt layer separating an n-octahedra thick perovskite block. Conventional high-temperature oxide synthesis methods cannot access members with n > 3, but low-temperature layer-by-layer thin film methods allow the preparation of materials with thicker perovskite blocks, exploiting high surface mobility and lattice matching with the substrate. This paper describes the growth of an n = 6 member CaO[(CSMO)2(LCMO)2 (CSMO)2] in which the six unit cell perovskite block is sub-divided into two central La0.67Ca0.33MnO3 (LCMO) and two terminal Ca0.85Sm0.15MnO3 (CSMO) layers to allow stabilization of the rock salt layer and variation of the transition metal charge.

Abstract: We observe magnetic-field-induced level mixing and quantum superposition phenomena between three approaching single-particle states in a quantum dot probed via the ground state of an adjacent quantum dot by single-electron resonant tunneling. The mixing is attributed to anisotropy and anharmonicity in realistic dot confining potentials. The pronounced anticrossing and transfer of strengths (both enhancement and suppression) between resonances can be understood with a simple coherent level mixing model. Superposition can lead to the formation of a dark state by complete cancellation of an otherwise strong resonance, an effect resembling coherent population trapping in a three-level-system of quantum and atom optics.

Abstract: In this contribution, we present an overview of the various techniques for combining atomistic molecular dynamics with Monte Carlo simulations, mainly in the context of condensed matter systems, as well as a brief summary of the main accelerated dynamics techniques. Special attention is given to the force bias Monte Carlo technique and its combination with molecular dynamics, in view of promising recent developments, including a definable timescale. Various examples of the application of combined molecular dynamics / Monte Carlo simulations are given, in order to demonstrate the enhanced simulation efficiency with respect to either pure molecular dynamics or Monte Carlo.