2026-09-09 | | Total: 118
We provide the first fully rigorous justification of time-dependent density functional theory in the continuum: Given a time-dependent density, we prove that there exists an external potential, unique up to a time-dependent constant, such that the corresponding Schrödinger equation reproduces the given density. This potential can be obtained with an iteration scheme. Our argument covers Coulomb interactions and does not need any Taylor expansion in time. It instead requires analyticity in space for the density and potential, which is compatible with extended nuclei.
We report a neutron diffraction study of the structural and magnetic properties of a family of ferro- and ferrimagnetic Prussian blue analogues (PBAs), $A_4$[Fe(CN)$_6$]$_{2.7}$ ($A$ = Co, Mn, Ni), Rb$_2$Ni$_4$[Fe(CN)$_6$]$_{3.3}$, and Mn$_4$[Cr(CN)$_6$]$_{2.7}$, as a function of temperature (2--450~K) and applied magnetic field. All bimetallic compounds of the Fm$\overline{3}$m family exhibit a broad diffuse feature at low scattering angle, which we identify, through comparison with the cation-stabilized Rb$_2$NiFe framework, as an intrinsic signature of correlated vacancies and their associated interstitial water. High-temperature diffraction reveals a continuous crossover from positive to negative thermal expansion in CoFe and MnFe upon dehydration, while NiFe remains structurally robust up to 450~K. At low temperature, all compounds order in a collinear ferrimagnetic state with propagation vector $\mathbf{k}=(0,0,0)$, except MnFe, which adopts a partially frustrated magnetic structure with $\mathbf{k}=(1,0,0)$. A moderate magnetic field of $B_c = 1.3$~T drives a spin reorientation in MnFe toward the collinear $\mathbf{k}=(0,0,0)$ ferrimagnetic state common to the other compounds; the same transition is independently induced by dehydration. A minimal Heisenberg model shows that this transition results from a near-compensation between antiferromagnetic Mn--Fe coupling and a geometrically frustrated antiferromagnetic Mn--Mn interaction on the face-centered-cubic Mn sublattice, placing MnFe in the vicinity of a magnetic compensation point. These results resolve a longstanding ambiguity in the interpretation of the Fe $K$-edge XMCD response of MnFe-based PBAs, and establish water content as a key parameter controlling both the structural and magnetic stability of this family of materials, with direct relevance to their use as battery electrodes.
Despite the pivotal role of twist in modulating physical properties at van der Waals (vdW) interfaces, the mechanics governing torsional response remain poorly understood. Here, we probe twist mechanics at homo- and heterogeneous vdW interfaces, together with their sliding behaviors within a unified experimental framework. For both systems, the peak torque scales nearly linearly with contact area, in contrast to predictions from linear elastic and rigid models. Remarkably, while the sliding friction of the two interfaces diverges by over three orders of magnitude owing to different scaling laws, the corresponding torque follows the same linear scaling and differs by only about twenty-fold. Large-scale atomistic simulations reveal an edge-dominated yielding mechanism for torsional motion, wherein elastic reconstruction shifts the effective load-bearing region toward the edges, eliminating torque from the contact interior. This mechanism contrasts with the bulk-mediated stress transmission governing translational sliding, a distinction rooted in the different loading geometries inherent to the two motion modes, where torsional loading necessitates perimeter actuation, whereas sliding enables center-driven loading. This symmetry-imposed divergence demonstrates that translational and torsional properties cannot be predicted from one another at vdW interfaces, providing critical insights for the design of dynamically reconfigurable micro- and nanoelectromechanical devices.
Molar heat capacity at constant pressure, C_p (T), is a central thermodynamic quantity in materials science, but it remains difficult to calculate accurately from first principles when anharmonic vibrational effects become important. In this work, rocksalt ZrC is used as a high-temperature benchmark for the Beyond Quasi-Harmonic (BQH) method, a prepared-supercell first-principles approach that extracts anharmonic vibrational energy directly from density-functional-theory energy evaluations. The calculated BQH molar heat capacity is compared with CALPHAD reference values, a fully anharmonic thermodynamic-integration calculation, a calculated C_v (T) curve from the present work, and a quasi-harmonic calculation. The BQH curve gives a substantial improvement beyond C_v (T) and beyond the quasi-harmonic approximation. Because ZrC is electronically conductive, an electronic heat-capacity correction was also estimated from the density of states at the Fermi level and added to the vibrational BQH result. This correction brings the calculated molar heat capacity into close agreement with the fully anharmonic theoretical benchmark through approximately 1200 K. These results show that the BQH method fully captures important anharmonic phonon contributions missing from ordinary quasi-harmonic calculations.
Hydrogen isotope retention in plasma facing and structural alloys is a central materials challenge for deuterium-tritium fusion. Motivated by ion beam irradiation experiments and first principles calculations, we identify an irradiation-enabled mechanism whereby solute segregation to defect sinks enhances hydrogen isotope trapping. Triple-ion irradiation of reduced activation ferritic-martensitic steel F82H reveals pronounced segregation of Cr and Ta to cavity surfaces. Density functional theory shows that these segregants markedly increase H stability at cavities by strengthening binding energies and increasing migration barriers, leading to substantially higher attention. The effect originates from solute-tuned electronic structure: Ta promotes strong H 1s-metal d orbital hybridization, whereas Cr shifts the surface d-band toward a more favorable bonding configuration. These findings provide an atomistic link between irradiation-induced segregation and elevated-temperature hydrogen isotope retention and alloy chemistry routes to control tritium inventory in fusion environments.
Within the rapidly growing family of two-dimensional (2D) materials, 2D ferromagnetic (FM) materials have attracted a great deal of recent attention. The ferromagnetism in 2D FM materials originates from the underlying exchange interactions. In particular, metal-insulator transitions (MITs) in 2D FM materials induced via various modulation strategies is critical to the design of 2D electronic devices with reduced sizes (e.g. transistor). However, the critical behaviors of exchange interactions and their effects on the Curie temperatures near MITs remains to be clarified. By taking the 2D FM 1T-FeX2 (X=Cl, Br, I) monolayers as the testbeds, we computationally study the critical behaviors of the exchange interactions in the vicinity of the MITs driven by the on-site Coulomb electron-electon interactions. The FM superexchange interactions and the direct exchange interactions are found to energetically favor the halfmetallic and semiconducting states, respectively. In addition, the halfmetallic states have the stronger FM couplings and thus the higher Curie temperatures. This study provides a comprehensive understanding of the interrelations between MITs and exchange interactions and could be useful for the device design requiring both MITs and exchange-interaction related magnetic properties, particularly robust magnetic ordering and higher Curie temperatures, in 2D FM materials.
Metal phthalocyanines (MPc) are a versatile molecular platform for applications ranging from organic optoelectronic devices to single-atom catalysis (SAC). Their adsorption and layer formation on the prototypical transparent electrode substrates of organic optoelectronic devices is directly relevant for charge injection and transport across the organic-oxide interface. Moreover, the well-defined M-N$_4$ coordination of the metal cation defines their activity as SACs for (electro-)catalysis. Here, the adsorption of the metal-free phthalocyanine (H$_2$Pc) is characterized on In$_2$O$_3$(111) using low-temperature STM, nc-AFM, and STS. In$_2$O$_3$ is not only a model system of indium tin oxide (ITO) but also an active catalytic material for CO$_2$ reduction. H$_2$Pc adsorbs in the same site and configuration reported for copper phthalocyanine [J. Mater. Chem. C 13, 17650-17661 (2025)] and for the majority of cobalt phthalocyanine [Surf. Sci. 722, 122065 (2022)]. Despite this shared preference in adsorption site, H$_2$Pc cannot be organized into extended and ordered monolayer structures by gentle annealing: the molecule starts to decompose at $\approx$50$^\circ$C, well below the temperature used to grow monolayers of CoPc and CuPc. Self-metalation is not observed on stoichiometric In$_2$O$_3$(111). On the reduced surface where In$^0$ adatoms are present, new H$_2$Pc-related features appear but cannot be identified by imaging only. The comparison of H$_2$Pc with CoPc and CuPc identifies distinct roles of the central metal ion in the metal-Pc/In$_2$O$_3$(111) systems: it acts as a structural anchor that stabilizes the macrocycle against decomposition on the surface, and it modifies the frontier-orbital character in ways that determine whether a second adsorption configuration is populated.
SnS as an environmentally friendly, cost-effective, and earth-abundant narrow-bandgap semiconductor material, has demonstrated significant application potential in the field of medium-temperature thermoelectric conversion. However, the thermoelectric performance of its bulk counterpart is inherently constrained by intrinsic point defects (e.g., vacancies) and the material's specific band structure. Low-dimensional engineering has emerged as a pivotal strategy for overcoming these limitations and enhancing thermoelectric performance. In this work, we systematically investigate the thermoelectric properties of SnS nanofilms with distinct thicknesses (82 nm, 199 nm, 616 nm, and 813 nm) across a temperature range of 300-600 K. Measurements were conducted using time-domain thermoreflectance (TDTR) and a dedicated thin-film thermoelectric parameter test system (ZEM-3). Our results confirm that low-dimensionalization effectively boosts the thermoelectric performance of SnS, with the thermoelectric figure of merit (ZT) displaying a pronounced dependence on both film thickness and temperature. All four SnS thin films exhibit thermoelectric performance that is markedly superior to that of bulk SnS. This enhancement is primarily attributed to the quantum confinement effect, energy filtering effect, and intensified phonon scattering, all of which are induced by the low-dimensional structural characteristics. This work provides not only experimental evidence and theoretical insights for the performance optimization of SnS nanofilms but also establishes a foundational framework for the development of high-efficiency, eco-friendly medium-temperature thermoelectric materials, thereby holding significant scientific value and practical implications.
The origin of ferromagnetism in the stable $β$-UH$_3$ phase is still debated. With a high Curie temperature and a short uranium-uranium interatomic distance, $β$-UH$_3$ is placed outside the known systematics of uranium compounds. Pressure provides a powerful means of tuning magnetism toward its suppression, thereby offering valuable insights into its underlying mechanisms, yet the high-pressure behavior of $β$-UH$_3$ has remained largely unexplored. Here, we combine in the Diamond Anvil Cell (DAC) the development of pure $β$-UH$_3$ synthesis and widefield nitrogen-vacancy (NV) center magnetometry to measure the pressure dependence of the Curie temperature in $β$-UH$_3$ up to about 20~GPa. We present two data analysis methods, vector magnetic field reconstruction and statistical analysis of the optically detected magnetic resonance (ODMR) response, which enable us to directly image the magnetic dipole of $β$-UH$_3$ during isobaric warming. We observe a linear decrease in the Curie temperature with pressure, yielding a slope of $\text{d}T_C/\text{d}P=-3.86\,(8)$~K/GPa. Ferromagnetism in $β$-UH$_3$ is predicted to reach 0~K at approximately 44~GPa, where a quantum critical point may emerge.
Time-dose reciprocity, commonly associated with the Bunsen-Roscoe law, states that the response of a photosensitive system depends only on the total exposure dose, regardless of how that energy is delivered over time. Light-sensitive processes in mixed-halide perovskites, such as photoinduced halide segregation, often exhibit threshold-like behavior that may violate this principle and enable material-state control by photon timing. We test this using pulse-burst excitation, which introduces an additional temporal control dimension beyond conventional parameters such as pulse fluence, repetition rate, and average power. By redistributing the same photon dose over microsecond-to-millisecond timescales, we create distinct nonequilibrium excitation conditions and show that mixed-halide perovskites can evolve into different metastable states, revealing a breakdown of time-dose reciprocity in the combined processes of halide segregation and remixing. This additional temporal degree of freedom not only enables control of the material state but also provides a new experimental framework for disentangling the competing processes underlying photoinduced halide redistribution. Our findings establish photon timing as a control parameter for perovskite photochemistry and open additional opportunities for optical memory and neuromorphic photonic applications.
The recent advances in research on layered magnets have established them as invaluable materials for spintronic applications and precursors to 2D magnets. The problem is that the number of such materials available for experimentation is still low. To tackle the problem, one may employ phase engineering; in particular, the interplay between the phase stability and dimensionality provides great opportunities for design of layered magnets. Here, this interrelation is harnessed to produce epitaxial films of an overlooked layered polymorph of GdAlGe. The compound is formed by alternating triangular Gd and honeycomb AlGe lattices. It is stable in ultrathin films, up to 5 monolayers. In thicker films, an additional unwanted phase emerges. The structural quality of the ultrathin films is witnessed by electron and X-ray diffraction as well as electron microscopy. Magnetization measurements detect an anisotropic ferromagnetic state. The GdAlGe magnetism is corroborated by a study of electron transport revealing the anomalous Hall effect and negative magnetoresistance. The material is naturally integrated with Ge, a technological semiconductor, which should facilitate applications. The present study can serve as a blueprint for design of layered magnets via polymorph engineering.
Alloying in two dimensional semiconductors is widely used to tune bandgaps, yet its implications for charge and energy transport remain poorly understood. Here, we investigate MoS2xSe2(1-x) alloys as a model system to study the interplay between composition, thickness, and disorder. Optical transitions and valence band dispersions evolve continuously with both stoichiometry and number of layers, with negligible bandgap bowing and a composition dependent attenuation of thickness driven renormalization. In contrast, time resolved spatial mapping of photoexcited carriers reveals a pronounced and asymmetric collapse of carrier diffusivity at intermediate compositions, which cannot be accounted for by changes in effective mass or band alignment, and instead emerges from strong real space fluctuations in the local energetic landscape generated by random chalcogen substitution. Microscopic simulations reproduce the experimental trends and show that the character of disorder depends critically on the direction of alloying, producing either scattering barriers or deep trapping sites. Together, these results demonstrate that transport in TMD alloys is governed by disorder physics, overlooked by conventional optical and photoemission probes at equilibrium. Our findings establish transport as a stringent metric of electronic quality and high-light intrinsic limitations in the usage of TMD alloys for layered semiconductor devices.
Spin-degenerate antiferromagnets can host hidden unconventional magnetism in their lattice degrees of freedom. We show that coherent phonons activate this magnetism by removing the spin-layer operations that enforce equilibrium band degeneracy without changing the collinear Néel order. The frequency and polarization of a pump electric field select a resonant \(Γ\)-point optical phonon and, within a doublet, its coordinate direction. This choice fixes the residual spin-layer symmetry. In monolayer MnPSe$_3$, an \(A_{2u}\) mode produces \(i\)-wave splitting odd in the mass-weighted phonon coordinate $Q$, reversing sign under $Q\to-Q$, whereas two orthogonal directions of the same doubly degenerate \(E_u\) doublet produce \(d\)- and \(s\)-wave splitting at a common resonance. Rotating the in-plane pump field \(\bm E_\parallel\) therefore programs both the spin-splitting texture and the thermoelectric spin current, continuously tuning the response between transverse pure-spin and longitudinal spin-polarized currents. A complete classification of two-dimensional collinear spin layer groups identifies the \(Γ\)-point coordinates that remove the degeneracy-enforcing operations and activate such unconventional magnetism.
In-situ rapid solidification studies demand measurements of thermal histories with high temporal resolution. We present a simple, effective setup to quantify the cooling response of an uncoated commercial Protochips Fusion MEMS chip in an in-situ scanning transmission electron microscopy (STEM) context. We drive user-defined temperature programs via an arbitrary waveform generator (AWG), while recording the voltage drops across a series shunt to reconstruct chip resistance and temperature at sub-millisecond resolution. We confirm the response times inferred from the current; however, the temperature obtained from the physically linked resistance, $T(R)$, evolves more slowly. Analysis of the maximum cooling step reveals an exponential-like relaxation with time constant $τ=1.80$ ms, consistent with reported thermal lag constants for fast scanning calorimetry. From the time to reach $95\%$ of the temperature difference $ΔT$, we measure an average cooling rate of $\approx 7.9\times 10^{4}$ K/s. Robustness checks include repeated $R(T)$ measurements (revealing a modest downward drift approaching an asymptote), a 10 k$Ω$ test load, and characterization of small off-duty arbitrary waveform generator leakage/offsets. These findings define practical bounds on achievable thermal-path rates when planning in-situ electron microscopy experiments with this chip platform.
Hydrogen embrittlement limits the deployment of high-strength steels and advanced alloys in hydrogen infrastructure. Laser shock peening (LSP) is increasingly considered as a mitigation route because it combines deep compressive residual stresses with near-surface microstructural modification. This review critically assesses LSP not as an isolated strengthening treatment, but as a surface layer design strategy governed by coupled stress, hydrogen transport and trapping mechanisms. Evidence from steels, nickel based alloys and additively manufactured materials shows that compressive residual stresses may suppress stress assisted hydrogen transport and delay crack initiation, while LSP-induced nanostructuring, dislocations, twins and interfaces can either redistribute hydrogen beneficially or promote localized plasticity and damage. Reported trends are frequently confounded by hydrogen charging mode, surface roughness, contamination, residual stress depth profiling and limited structure performance correlations. Two design critical gaps are identified: the lack of quantitative links between post-LSP stress/defect architectures and hydrogen assisted fatigue crack growth, and the near absence of impact toughness data after LSP under hydrogen exposure. The review proposes mechanism informed qualification routes combining residual stress mapping, hydrogen characterization, service representative mechanical testing and modelling.
Combinatorial processing platforms (CPPs), integrating Si microtip arrays with combinatorial thin film synthesis and atom probe tomography (APT), enable near-atomic-scale characterization of compositionally complex solid solutions (CCSSs) under diverse processing and reaction conditions, including oxidation, thermal phase stability and electrocatalytic reactions. Their application at elevated temperatures, however, can be limited when CCSS constituents such as Pd and Pt react with the Si support to form silicides. Although thermally grown SiO2 has proven effective as a diffusion barrier between pure Pt and Si, its performance for multicomponent CCSS thin films is unclear. Here, using Ag-Au-Pd-Pt as a model system, we compare a 25 nm thermally grown SiO2 barrier with native Si oxide during annealing using APT and transmission electron microscopy. Native Si oxide prevents detectable interfacial reactions up to 300°C, but at 400°C Pd and Pt react with Si, causing silicide formation and substantial redistribution of the film constituents. At 600°C, extensive substrate reactions disrupt the CCSS film and produce a pronounced needle-shaped silicide morphology. In contrast, thermally grown SiO2 suppresses CCSS thin film-substrate reactions up to 600°C and retains the CCSS composition. The thermally grown SiO2 thus extends the applicable temperature range of Si-based CPPs to at least 600°C for near-atomic-scale characterization of CCSS thin films.
We investigate the excitonic properties of perylene red microcrystals, whose unit cell contain eight molecules, using both experimental and theoretical methods. Only two of the nominal eight Davydov transitions are experimentally observed, with an apparent splitting of 610 cm$^{-1}$ that is reasonably reproduced by a Frenkel-Holstein exciton model parametrized from density functional theory. A symmetry-based analysis of the eigenstates reveals that the two observed transitions do not belong to the same Davydov pair, so the splitting does not directly report on the Coulomb coupling within the unit cell. The complex mixture of local excitations produces a J-like band for the dominant transitions, favouring exciton transfer through increased spectral overlap. Consistent with this picture, time-correlated single-photon counting of the emission and ultrafast transient absorption spectroscopy show that the dynamics is dominated by highly mobile singlet Frenkel excitons. The exciton mobility extracted from experiment agrees very well with that obtained from Kinetic Monte Carlo simulations, supporting a picture of incoherent hopping transport.
Dynamical electron backscatter diffraction (EBSD) simulations increasingly support systematic studies of how material and experimental parameters affect diffraction patterns. Investigations that combine composition and crystal structure, including changes in atomic order, with accelerating voltage, temperature, and many orientations or detector conditions can require large sets of related calculations. Several packages provide dynamical EBSD simulation; we use the open-source, scriptable EMsoft programs because they are well suited to large studies on high-performance computing systems. An EMsoft investigation nevertheless spans separate configuration, structure, scattering, and orientation files and several dependent calculation stages. As variations multiply, manual preparation becomes difficult to check, interruptions obscure which results remain usable, and the origin of individual patterns becomes laborious to reconstruct. With EBSDmagus, researchers define fixed and varying parameters once. The software prepares the required EMsoft calculations, reuses compatible intermediate results, and checks that the generated files represent the requested investigation before execution. After execution, it checks the expected outputs and, following an interruption, resumes only unresolved calculations while preserving completed work. A portable run record links each output to its inputs and calculation history, providing the provenance needed when the results are deposited as findable, accessible, interoperable, and reusable (FAIR) data. We assess this approach using a completed representative calculation, a larger study prepared and checked before submission, two observed cluster interruptions, and a controlled job cancellation. During recovery, EBSDmagus retained completed calculations and resubmitted only incomplete work.
We present a combined density functional theory (DFT) and density functional tight binding (DFTB) study of zirconia (ZrO2) nanoparticles of experimentally relevant sizes of several nanometers and their interactions with the CO molecule. A hybrid DFTB - Force Field (DFTB-FF) framework is developed, whereby band structure calculations rely on an existing Slater-Koster framework, while the accuracy of structural optimization and adsorption properties is controlled by the introduction of classical long-range interatomic potentials into DFTB instead of the traditional repulsive potentials. Additionally, coordination-dependent Zr-C potentials are introduced to account for the distinct local chemical environments of bulk-like facet sites and under-coordinated tip and edge sites, thereby improving the description of CO adsorption. This hybrid DFTB-FF approach substantially improves the robustness of geometry optimization and provides a practical strategy for extending the applicability of DFTB to complex oxide nanostructures. The calculations reveal that termination stoichiometry can be used to engineer intrinsic, p-type, or n-type electronic structures and thereby tune the adsorption activity of zirconia nanoparticles. While stoichiometric nanoparticles do not activate the C-O bond, low-coordinated sites in Zr-rich (n-type) nanoparticles exhibit chemisorption accompanied by charge donation into a CO antibonding LUMO-derived orbital, resulting in C-O bond activation. These results demonstrate that stoichiometry-controlled electronic structure and under-coordinated surface sites introduced by nanostructuring play a key role in governing the adsorption strength and reactivity of zirconia nanoparticles.
Strain-gradient engineering via out-of-plane wrinkling offers a powerful route to tune electronic and electromechanical properties in two-dimensional (2D) materials. Here, we systematically investigate the electronic and flexoelectric response of wrinkled Janus MoSSe/MoSeS monolayers using density functional theory (DFT) calculations coupled with continuum elastica modeling. Exploring varying wrinkle sizes and compressive strain levels ($5\%-20\%$), we show that the global out-of-plane polarization follows a linear behavior when parameterized by the projected aspect ratio of the nanowrinkles. On this basis, we develop a physically grounded geometric model incorporating an effective 2D flexoelectric coefficient, which accurately predicts DFT polarizations without requiring higher-order nonlinear parameters. Atom- and orbital-resolved charge density analyses reveal the microscopic origin of this behavior: while the central Mo $4d$-orbital manifold acts as a robust, linear flexoelectric core, local curvature drives continuous, chemically asymmetric charge transfer between the S $3p$ and Se $4p$ sublayer manifolds. Our findings establish clear geometric design rules for harnessing Janus-based flexoelectricity for flexible nanoelectronics and optoelectronics.
Wurtzite ferroelectrics combine large remanent polarization with full CMOS compatibility, positioning them as a leading platform for next-generation non-volatile memory. Their practical deployment, however, is hindered by an intrinsically large coercive field, rooted in the high energy barrier separating the polar wurtzite phase from the nonpolar hexagonal phase that mediates polarization switching. Here, using first-principles calculations, we propose an alternative, field-free strategy for lowering this barrier: above-bandgap electronic photoexcitation. Taking LaN as a representative wurtzite ferroelectric, we show that light induced carriers dramatically reduce the energy difference between the hexagonal intermediate phase and the wurtzite ground state, sharply reducing the energy barrier to ferroelectric switching. This effect originates from a photoinduced partial metallization of the polar phase, which screens the dipole-dipole interactions that stabilize ferroelectric order and thereby favors the competing nonpolar structure. The robustness of this mechanism is further confirmed for the rocksalt polymorph. Our results establish light as a powerful, non-invasive route to controlling ferroelectric switching in wurtzites, opening a path toward faster, lower-voltage, and more energy-efficient non-volatile memory technologies.
We investigated the evolution of the well-known parabolic dependence of atomic volume on atomic number in transition metals under extreme compression at pressures up to 400 GPa using density functional theory calculations. Our results reveal that the ambient-pressure parabolic trend transforms into a characteristic cubic-like behavior at high pressures. This evolution is attributed to the higher compressibility of bcc transition metals associated with comparatively large increases in the total energy. The present findings are discussed in relation to previous experimental observations and first-principles calculations.
Molecular beam epitaxy (MBE) is renowned for its potential for atomic layer control, but unexpected growth mechanisms can potentially compromise this level of precision. In this study, we employ first-principles calculations to investigate the atomistic processes governing the MBE growth of perovskite SrTiO3 and Ruddlesden-Popper Sr2TiO4 films on a SrTiO3 substrate. We systematically explore the potential molecular species in the gas phase and their reactions and diffusion dynamics on the film surfaces and layer edges. Our analyses uncover three mechanisms of importance for understanding this type of growth. First, oxygen vacancies can be dynamically induced during surface diffusion on defect-free substrate and noticeably accelerate the diffusion processes. Second, while the SrO layer is expected to grow in a single-layer growth mode, the presence of potential TiSr defects may promote the formation of SrO islands. Lastly, adsorbed Ti atoms and TiO2 molecules on the SrO bilayer can insert into SrO bilayers, resulting in an unexpected growth sequence. These findings may have broader implications for the MBE growth of metal oxide films and provide guidance for achieving improved control over their growth processes.
We present a first-principles real-time approach to study non-linear response of solids at finite temperature. Finite temperature effects are included as renormalization of the quasiparticle energies and a dephasing term proportional to the quasiparticle lifetimes. We evaluate electron-phonon matrix elements from Density-Functional Perturbation Theory for lattice dynamics and then calculate quasiparticles renormalization and lifetime from the Fan and Debye-Waller terms of the electron self-energy. Electron excitations are treated at the independent particle level of approximation. We apply the approach to the second-harmonic generation (SHG) in monolayer and trilayer MoS$_2$. We observe a nontrivial temperature-dependence of the SHG due to a strong crystal-momentum dependent quasiparticle renormalization. From the phonon-mode analysis we find that the coupling with acoustic and shear modes determines the overall crystal-momentum dependence respectively in monolayer and trilayer MoS$_2$. The nontrivial temperature-dependence of the SHG can help rationalize the increase of SHG intensity with increasing temperature observed in monolayer MoS$_2$ at a given laser energy [Adv. Optical Mater. 8, 2000441 (2020)].
AlGaN alloys are promising wide and ultra-wide-bandgap semiconductors for next-generation power and RF electronics applications. To realize the full potential of AlGaN based devices, it is important to understand the electron transport to accurately predict device performance and identify material limits under various operating conditions. In this work, the high-field electron transport properties of AlxGa1-xN are investigated using a supercell based full band Monte Carlo method. The supercell approach is employed to explicitly capture the true disorder of the alloy system, enabling a more realistic description of carrier transport. The velocity field characteristics are calculated across a range of Al compositions to evaluate key transport metrics, including peak velocity, saturation velocity, and critical electric field. The role of different scattering mechanisms is studied in detail to understand the high field transport mechanism in the AlGaN alloy system. In addition to steady state transport, transient electron dynamics are examined for various Al fractions to study velocity-overshoot behavior, which is especially important for improving the performance of scaled RF devices. Finally, the temperature dependence of the velocity field characteristics in ultra-wide-bandgap Al0.75Ga0.25N is investigated to assess its transport performance under high temperature conditions. These results provide a detailed understanding of high-field transport in AlGaN alloys and offer guidance for the design of AlGaN-based RF and power electronic devices.