2025-04-25 | | Total: 32
Characterizing microstructural effects on the dynamical response of materials is challenging due to the extreme conditions and the short timescales involved. For example, little is known about how grain boundary characteristics affect spall strength. This study explores 6H-SiC bicrystals under shock waves via large-scale molecular dynamics simulations. We focused on symmetric tilt grain boundaries with a wide range of misorientations and found that spall strength and dynamical fracture surface energy are strongly affected by the grain boundary microstructure, especially the excess free volume. Grain boundary energy also plays a considerable role. As expected, low-angle grain boundaries tend to have higher spallation strengths. We also extracted cohesive models for the dynamical strength of bulk systems and grain boundaries that can be used in continuum simulations.
Spintronics is a rapidly evolving technology that utilizes the spin of electrons along with their charge to enable high speed, low power and non volatile electronic devices. The development of novel materials with tailored magnetic and electronic properties is critical to exploit the full potential of spintronic applications. Among these, Heusler alloys stand out due to their tunable multifunctional properties. This review presents a comprehensive overview of various Heusler based materials including half metallic ferromagnets, spin gapless semiconductors, magnetic semiconductors, spin semimetals, and nearly zero moment materials focusing on their synthesis, structural and magnetic characterizations, and transport behavior. The role of crystal structure, and structural disorder in governing their magnetic and electronic properties is discussed in detail. Emphasis is placed on experimental results and their implications for spintronic devices. By bringing together recent advancements, the review highlights the critical role of Heusler alloys in advancing the next-generation spintronic technologies and outlines future directions for their integration in practical applications.
We present in depth study of p-type doping concentration of mid-wave infrared (MWIR) and long-wave infrared (LWIR) mercury cadmium Telluride (HgCdTe) thin films. Annealing time was changed under specific conditions to achieve a stable copper (Cu) doping concentration for HgCdTe thin films. Both MWIR and LWIR HgCdTe material were grown by molecular beam epitaxy (MBE), where different trends were observed between LWIR and MWIR HgCdTe thin films by increasing anneal time. We also report the impact of different thickness (4 micron, 6 micron and 9 micron) along with annealing time on doping level of LWIR HgCdTe thin films.
We implement a new experimental approach to directly measure the lifetime of guided polaritons arising from the strong-coupling of GaN excitons and the guided photonic modes of a slab waveguide. Using a Fourier imaging setup, combined with spatial filtering of the emission, the emission associated to polaritonic modes with well-defined propagation constants can be selectively analyzed in the temporal domain. By directing it to the entrance slit of a streak camera, time-resolved photoluminescence (TRPL) measurements along the polariton dispersion branch were performed at 40 K, enabling to assess the time decay of polariton modes. By combining this information with the photonic/excitonic fraction corresponding to each polariton mode, extracted from a coupled-oscillators model that indicate a Rabi splitting of $\Omega$ = 80 meV, we could extract the photon lifetime in the waveguide $\tau_\gamma\, =\, 3\pm 1$ ps. This corresponds to a record $Q$-factor in the UV of 16 000. The excitonic reservoir lifetime, which contributes to polariton formation, was determined through TRPL measurements on excitonic luminescence. Finally, measurements conducted at lower temperature highlight secondary feeding mechanisms for the guided polaritonic mode, either via photon recycling from the AlGaN cladding layer or through resonant injection of photons from transitions below the band gap.
This work presents a fully physical model of the hydrogen diffusion and trapping kinetics in metals, integrating permeation and thermal desorption within a unified framework. Based on the McNabb and Foster approach, it requires only binding energy and number density of trap sites. It correctly reproduces the physics of the system and the results of the analytical solutions of the permeation kinetics. It is also capable of reproducing thermal desorption spectra with considerable accuracy. The sensitivity analysis has elucidated the relationships among the processing conditions and the parameters commonly used to characterize permeation and thermal desorption experiments. An equation empirically derived from the simulation results, expressing the dependence of time lag in desorption on specimen thickness, number density of occupied trap sites, and cathodic concentration, is proposed. In summary, the model represents a valuable tool in supporting the interpretation and rationalization of experiments also from a quantitative viewpoint.
Zeolites, as representative porous materials, possess intricate three-dimensional frameworks that endow them with high surface areas and remarkable catalytic properties. There are a few factors that give a huge influence on the catalytic properties, including the size and connectivity of these three-dimensional channels and atomic level defects. In additional to that, the surface morphology and thickness variation of zeolites particles are essential to their catalytic performances as well. However, it is a significant challenge to characterize these macroscopic properties of zeolites using conventional techniques due to their sensitivity to electron beams. In this study, we introduce surface-adaptive electron ptychography, an advanced approach based on multi-slice electron ptychography, which enables high-precision reconstruction of both local atomic configurations and global structural features in zeolite nanoparticles. By adaptively optimizing probe defocus and slice thickness during the reconstruction process, SAEP successfully resolves surface morphology, thickness variations and atomic structure simultaneously. This integrated framework facilitates a direct and intuitive correlation between zeolite channel structures and particle thickness. Our findings open new pathways for large-scale, comprehensive structure property analysis of beam-sensitive porous materials, deepening the understanding of their catalytic behavior.
We report the design, synthesis, crystal structure, and physical properties of a layered intergrowth compound, Eu$_2$CuMn$_2$P$_3$. The structure of Eu$_2$CuMn$_2$P$_3$ features an alternating arrangement of hexagonal EuCuP block layers and trigonal EuMn$_2$P$_2$ block layers, interconnected through shared Eu planes. This structural hybridization leads to multiple magnetic orderings in Eu$_2$CuMn$_2$P$_3$: weak antiferromagnetic (AFM) ordering of Mn at $T_\mathrm{N}^\mathrm{Mn}$ = 80 K, AFM ordering of Eu at $T_\mathrm{N}^\mathrm{Eu}$ = 29 K, a spin-reorientation transition at $T_\mathrm{SR}$ = 14.5 K, and weak ferromagnetism below $T_\mathrm{N}^\mathrm{Mn}$. The spin configurations at different temperature regions were discussed based on the calculations of magnetic energies for various collinear arrangements. Resistivity measurements reveal a pronounced transition peak at $T_\mathrm{N}^\mathrm{Eu}$, which is suppressed in the presence of a magnetic field, resulting in a significant negative magnetoresistance effect. The computed semimetallic band structure, characterized by a small density of states at the Fermi level, aligns well with experimental observations. The successful synthesis of Eu$_2$CuMn$_2$P$_3$ and its fascinating magnetic properties highlight the effectiveness of our block-layer design strategy. By assembling magnetic block layers of compounds with compatible crystal symmetries and closely matched lattice parameters, this approach opens exciting avenues for discovering layered materials with unique magnetic behaviors.
We present an implementation of spin-orbit coupling (SOC) for the computation of nuclear magnetic resonance (NMR) chemical shielding tensors within linear response theory. Our implementation in the Vienna {\it Ab initio} Simulation Package (VASP) is tailored to solid-state systems by employing periodic boundary conditions and the gauge-including projector augmented waves (GIPAW) approach. Relativistic effects are included on the level of the zeroth-order regular approximation (ZORA). We discuss the challenges posed by the PAW partial wave basis in describing SOC regarding chemical shielding tensors. Our method is in good agreement with existing local-basis ZORA implementations for a series of Sn, Hg, and Pb molecules and cluster approximations for crystalline systems.
Driven quantum materials with on demand properties controlled by external stimuli are critical for emergent quantum technology. In optically tunable superconducting heterostructures, the lattice responses at the buried interface may hold the key to the light susceptibility but is very challenging to detect. In this work, a nondestructive synchrotron-based X-ray scattering phase-retrieval technique is implemented in monolayer-FeSe/SrTiO3 heterostructures to capture the three-dimensional interfacial atomic displacements in-situ as the interface superconductivity is actively manipulated by light. It is found that the interlayer sliding between FeSe and SrTiO3 can drastically alter how the lattice responds to the light. In domains with selected stacking configurations, the interface transforms the very weak photoexcitation in SrTiO3 into significant Fe-atom displacements in FeSe and generate metastable interfacial structures that can lead to a persistent superconductivity enhancement. These findings demonstrate an effective strategy for achieving greatly amplified light-lattice coupling for efficient quantum phase manipulations at designed interfaces.
Thermal quenching has been used to find metastable materials such as hard steels and metallic glasses. More recently, quenching-based phase control has been applied to correlated electron systems that exhibit metal--insulator, magnetic or superconducting transitions. Despite the discovery of metastable electronic phases, however, how metastability is achieved through the degrees of freedom, which can vary even at low temperatures such as those of an electron, is unclear. Here, we show a thermally quenched metastable phase in the Ising model without conservation of magnetization by Monte Carlo simulations. When multiple types of interactions that stabilize different long-range orders are introduced, the ordering kinetics divergently slow toward low temperatures, meaning that the system will reach a low temperature without ordering if the cooling rate is high enough. Quantitative analysis of the divergent behavior suggests that the energy barrier for eliminating the local structure of competing orders is the origin of this metastability. Thus, the present simulations show that competing interactions play a key role in realizing metastability.
The crux of the present study is to analyze the Mode I crack propagation behavior in a pre-stressed monoclinic crystalline strip of finite thickness and infinite extent. The investigation focuses on the effects of collinear Griffith cracks and dynamic punch loading induced by plane wave propagation. The cracks are assumed to be in motion, and a Galilean transformation is employed to formulate the problem within a moving coordinate system. The boundary value problem is transformed into a system of coupled Cauchy-type singular integral equations, which are solved analytically using the Hilbert transform method. This approach yields elegant closed-form solutions for both the stress intensity factor and the crack opening displacement. The study considers two monoclinic crystalline materials, Lithium Niobate and Lithium Tantalate, and compares their behavior with that of an isotropic material to assess the role of material anisotropy. Numerical simulations and graphical analysis are performed for the crystalline materials with monoclinic symmetry to evaluate the influence of crack velocity, punch loading, material anisotropy, initial stress, and crack geometry on the fracture parameters. As a special case, the system is analyzed under the action of point loading from the punch pressure, and a comparative assessment is conducted between point loading and constant normal punch pressure. The results unveil critical insights into the dynamic fracture behavior of anisotropic materials under localized loading. This understanding enhances failure prediction in high-precision fields such as geomechanics, MEMS, surface acoustic devices, and biosensors.
Ultrafast coherent control of electric dipoles using strong terahertz (THz) pulses provides a means to discover hidden phases of materials and potentially leads to applications in high-speed electro-optic devices. The effectiveness of this means, albeit demonstrated in architype (incipient) ferroelectric systems such as SrTiO3, hinges on a spectral overlapping between their soft phonon modes within the excitation bandwidth of THz pulses. Generally this can only induce an appreciable coupling close to the phase transition temperatures, where the lattice phonons substantially soften. Because of their emergent subterahertz collective dynamics, a fundamentally distinct and effective THz coupling mechanism can be envisaged in topological polar structures recently discovered in PbTiO3/SrTiO3 superlattices. Here, we show that polar skyrmions can be coherently driven into a hidden phase with transient macroscopic polarization, as probed based on THz field-induced second harmonic generation and optical Kerr effects. Such an ultrafast THz-driven phase transition is found to sustain across a broad temperature range of 4-470 K, in accordance with the equilibrium stability field of the skyrmions. Supplemented by dynamical phase-field simulations, we identify the spatial symmetries and relaxation behaviors of the excited collective modes, thereby revealing their correlation with the emergence of the polar phases. Our results unveil the exotic dynamical properties of topological polar structures, which could be technologically exploited given their remarkable flexibility in structure design and tunability under external fields.
The phase field approach to modeling fracture uses a diffuse damage field to represent a crack. This addresses the singularities that arise at the crack tip in computations with sharp interface models, mollifying some of the difficulties associated with the mathematical and numerical treatment of fracture. The introduction of the diffuse field helps with crack propagation dynamics, enabling phase-field approaches to model all phases of damage from crack initiation to propagation, branching, and merging. Specific formulations, beginning with brittle fracture, have also been shown to converge to classical solutions. Extensions to cover the range of material failure, including ductile and cohesive fracture, leads to an array of possible models. There exists a large body of studies of these models and their consequences for crack evolution. However, there have not been systematic studies into how optimal models may be chosen. Here we take a first step in this direction by developing formal methods for identification of the best parsimonious model of phase field fracture given full-field data on the damage and deformation fields. We consider some of the main models that have been used to model damage, its degradation of elastic response, and its propagation. Our approach builds upon Variational System Identification (VSI), a weak form variant of the Sparse Identification of Nonlinear Dynamics (SINDy). In this first communication we focus on synthetically generated data but we also consider central issues associated with the use of experimental full-field data, such as data sparsity and noise.
In this thesis, I discuss several basic science studies in the field of energy materials using neutron scattering as a probe for the lattice dynamics. To enable understanding of neutron scattering spectra, I also use computational and theoretical methods. These methods and neutron scattering in general are discussed in detail in Chapter 2. It is assumed that the reader is familiar with basic quantum mechanics as well as with solid state physics topics including the band theory of electrons, harmonic lattice dynamics, and molecular dynamics. For the unfamiliar reader, the details of electronic structure theory and lattice dynamics that are needed to understand the methods in Chapter 2 are provided in Chapters 3 and 4. In the remaining chapters, these methods are applied to the study of several energy materials: cuprate La2CuO4,(hybrid) solar perovskite CH3NH3PbI3, and thermoelectric clathrate Ba8Ga16Ge30.
Magnetic bilayers consisting of an epitaxially grown ferrimagnetic insulator and a heavy metal layer are attractive for spintronic application because of the opportunity for electric control and read-out of spin textures via spin orbit torque. Here, we investigate ultrathin thulium iron garnet (TmIG)/Pt bilayers when the TmIG layer thickness is 3 nm and below using a sensitive Sagnac magneto-optical Kerr effect technique. We compare the hysteresis loops from out-of-plane and in-plane applied magnetic fields. The preferred magnetization orientation evolves with the TmIG thickness and the presence of the Pt overlayer. We quantify the evolution of the magnetic anisotropy in these ultrathin films and find a significant change even when the TmIG thickness is varied by less than 1 nm. In these ultrathin films, the presence of a Pt overlayer changes the effective anisotropy field by more than a factor of 2, suggesting that the interfacial anisotropy at the Pt/TmIG interface plays a critical role in this regime.
In this work, we study the influence of disorder on the electronic structure of YbSb -- a rare-earth monopnictide featuring a simple rock-salt (B1) crystal structure and a well-defined Fermi surface topology -- by employing first-principles density functional theory (DFT). We focus on chemical disorder introduced through Te and Al doping, selected based on their thermodynamic stability in alloyed configurations, to understand how such perturbations modify the electronic states of YbSb. Our results indicate that Te doping predominantly introduces electron-like states at the \textit{X} and \textit{L} points, while Al doping leads to a suppression of hole-like states at $\Gamma$, effectively driving the system from a semimetallic state to one characterized by very narrow-gap behavior at $\Gamma$. This modulation of the Fermi surface, particularly the reduction of central hole pockets at $\Gamma$, plays a central role in altering inter-pocket scattering -- a mechanism critical for tuning quantum transport properties, including superconductivity. This disorder-driven modulation of the Fermi surface, particularly the suppression of central hole pockets at $\Gamma$, controls inter-pocket scattering, which is essential for optimizing quantum transport properties, including superconductivity. Our results show that disorder can be effectively used as a means of engineering band topology, thereby tuning quantum-related responses through tailored electronic structure.
The Helmholtz energy is a key thermodynamic quantity representing available energy to do work at a constant temperature and volume. Despite a well-established methodology from first-principles calculations, a comprehensive tool and database are still lacking. To address this gap, we developed an open-source Density Functional Theory Tool Kit (DFTTK), which automates first-principles thermodynamics using the quasiharmonic approximation (QHA) for Helmholtz energy predictions. This Python-based package provides a solution for automating the calculation and analysis of various contributions to Helmholtz energy, including the static total energy contributions at 0 K in terms of DFT-based energy-volume curves, vibrational contributions from the Debye-Gruneisen model and phonons, and thermal electronic contributions via the electronic density of states. The QHA is also implemented to calculate the Gibbs energy and associated properties at constant temperature and pressure. The present work demonstrates DFTTK's capabilities through case studies on a simple FCC Al and various collinear magnetic configurations of Invar Fe3Pt, where DFTTK enumerates all unique configurations and their associated multiplicities. DFTTK is freely available on GitHub, and its modular design allows for the easy addition of new workflows.
Dielectric spectroscopy is known as one of the most powerful techniques for studying ferroelectric and other polar materials. Since the discovery of ferroelectricity in Liquid Crystals, it has been successfully employed for the characterization of ferro-, antiferro- and ferri-electric liquid crystalline phases. However, recently the Boulder group raised the question of the applicability of dielectric spectroscopy for characterizing ferroelectric nematics due to parasitic effects from the insulating alignment layers. This affects the apparent/measured values of the dielectric permittivity. In this paper, we study this effect in greater detail. The following issues will receive special attention: Are the real values of dielectric permittivity lower or higher than the measured ones? Can the real values of dielectric permittivity be recovered in the cell with alignment layers? We also provide an example of the effect of insulating alignment layers in a non-ferroelectric nematic phase.
Non-trivial geometry of electronic Bloch states gives birth to topological insulators that are robust against sufficiently weak randomness, inevitably present in any quantum material. However, increasing disorder triggers a quantum phase transition into a featureless normal insulator. As the underlying quantum critical point is approached from the topological side, small scattered droplets of normal insulators start to develop in the system and their coherent nucleation causes ultimate condensation of a trivial insulation. Unless disorder is too strong, the normal insulator accommodates disjoint tiny topological puddles. Furthermore, in the close vicinity of such a transition the emergent islands of topological and trivial insulators display spatial fractal structures, a feature that is revealed only by local topological markers. Here we showcase this (possibly) generic phenomenon that should be apposite to dirty topological crystals of any symmetry class in any dimension from the Bott index and local Chern marker for a square lattice-based disordered Chern insulator model.
Graphite occurs in a range of metastable stacking orders characterized by both the number and direction of shifts between adjacent layers by the length of a single carbon-carbon bond. At the extremes are Bernal (or ``ABAB...'') stacking, where the direction of the interlayer shift alternates with each layer, and rhombohedral (or ``ABCABC...'') stacking order where the shifts are always in the same direction. However, for an N-layer system, there are in principle $N-1$ unique metastable stacking orders of this type. Recently, it has become clear that stacking order has a strong effect on the low energy electronic band structure with single-layer shifts completely altering the electronic properties. Most experimental work has focused on the extremal stacking orders in large part due to the difficulty of isolating and identifying intermediate orders. Motivated by this challenge, here we describe two atomic force microscopy (AFM) based techniques to unambiguously distinguish stacking orders and defects in graphite flakes. Photo-thermal infrared atomic force microscope (AFM-IR) is able to distinguish stacking orders across multiple IR wavelengths and readily provides absolute contrast via IR spectral analysis. Scanning microwave impedance microscopy (sMIM) can distinguish the relative contrast between Bernal, intermediate and rhombohedral domains. We show that both techniques are well suited to characterizing graphite van der Waals devices, providing high contrast determination of stacking order, subsurface imaging of graphene flakes buried under a hexagonal boron nitride (hBN) dielectric layer, and identifying nanoscale domain walls. Our results pave the way for the reliable fabrication of graphene multilayer devices of definite interlayer registry.
Pyrochlore oxides with chemical formula of A2B2O7 exhibit a diverse range of electronic properties as a representative family of quantum materials. These properties mostly stem from strong electron correlations at the transition metal B site and typical geometrical frustration effects on the pyrochlore lattice. Furthermore, the coupling between the magnetic moments of the rare-earth A site and the conduction electrons at the B site, along with the relativistic spin-orbit coupling particularly affecting the 4d/5d electrons at the B site, gives rise to the topological characteristics of the correlated electrons. This review paper focuses on the metal-insulator transitions in pyrochlore oxides as evidence of the strong electron correlation, which is highlighted as a rich source of intriguing charge dynamics coupled with frustrated spin-orbital entangled magnetism.
Synthetic polymeric materials underpin fundamental technologies in the energy, electronics, consumer goods, and medical sectors, yet their development still suffers from prolonged design timelines. Although polymer informatics tools have supported speedup, polymer simulation protocols continue to face significant challenges: on-demand generation of realistic 3D atomic structures that respect the conformational diversity of polymer structures. Generative algorithms for 3D structures of inorganic crystals, bio-polymers, and small molecules exist, but have not addressed synthetic polymers. In this work, we introduce polyGen, the first latent diffusion model designed specifically to generate realistic polymer structures from minimal inputs such as the repeat unit chemistry alone, leveraging a molecular encoding that captures polymer connectivity throughout the architecture. Due to a scarce dataset of only 3855 DFT-optimized polymer structures, we augment our training with DFT-optimized molecular structures, showing improvement in joint learning between similar chemical structures. We also establish structure matching criteria to benchmark our approach on this novel problem. polyGen effectively generates diverse conformations of both linear chains and complex branched structures, though its performance decreases when handling repeat units with a high atom count. Given these initial results, polyGen represents a paradigm shift in atomic-level structure generation for polymer science-the first proof-of-concept for predicting realistic atomic-level polymer conformations while accounting for their intrinsic structural flexibility.
Surface roughness plays a critical role in ultrashort pulse laser ablation, particularly for industrial applications using burst mode operations, multi-pulse laser processing, and the generation of laser-induced periodic surface structures. Hence, we address the impact of surface roughness on the resulting laser ablation topography predicted by a simulation model and compared to experimental results. We present a comprehensive multi-scale simulation framework that first employs finite-difference-time-domain simulations for calculating the surface fluence distribution on a rough surface measured by an atomic-force-microscope followed by the two-temperature model coupled with hydrodynamic/solid mechanics simulation for the initial material heating. Lastly, a computational fluid dynamics model for material relaxation and fluid flow is developed and employed. Final state results of aluminum and AISI 304 stainless steel simulations demonstrated alignment with established ablation models and crater dimension prediction. Notably, Al exhibited significant optical scattering effects due to initial surface roughness of 15 nm - being 70 times below the laser wavelength, leading to localized, selective ablation processes and substantially altered crater topography compared to idealized conditions. Contrary, AISI 304 with RMS roughness of 2 nm showed no difference. Hence, we highlight the necessity of incorporating realistic, material-specific surface roughness values into large-scale ablation simulations. Furthermore, the induced local fluence variations demonstrated the inadequacy of neglecting lateral heat transport effects in this context.
While nearly single cycle THz pulse is conventionally accepted as the stimulus for the fastest and the most energy efficient control of spins in magnets, all-dielectric metasurfaces have been recently demonstrated to be the least dissipative mean to enhance and control the coupling of light to spins. All-dielectric metasurfaces for the THz control of spins hold great potential in the field of spintronics and related technologies, pushing the boundaries of speed and energy efficiency in spin-based information processing. Here we demonstrate such a metasurface for an advanced THz control of spins in a ferrimagnetic film of iron garnet. Structuring a nonmagnetic substrate one can force a THz electromagnetic field, otherwise described by plane waves, to acquire an out-of-plane magnetic field and thus enable arbitrary direction of the torque acting on spins in all three dimensions. Hence, metaphotonics opens up a plethora of opportunities for advanced control of spins at THz rates in many hot fields of contemporary science, including spintronics, magnonics and quantum computing.
We study a graphene Josephson junction where the inner graphene layer is subjected to spin-orbit coupling by proximity effect. This could be achieved, for example, by growing the graphene layer on top of a transition metal dichalcogenide, such as WS$_2$. Here, we focus on the ballistic, wide, and short junction limits and study the effects of the spin-orbit interaction on the supercurrent. In particular, we analyze the current phase relation using an analytical approach based on the continuum model. We find combinations of types of spin-orbit coupling that significantly suppress the supercurrent by opening a gap in the graphene band structure. At the same time, other combinations enhance it, acting as an effective spin-valley resolved chemical potential. Moreover, we find that a strong Rashba spin-orbit coupling leads to a junction with a highly voltage tunable harmonic content.