Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Open quantum system”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,261 records · Page 70
All sourcesNASA NTRSDOE OSTI
DOE OSTI2019

Langevin simulations of a long-range electron-phonon model

We present a quantum Monte Carlo (QMC) study, based on the Langevin equation, of a Hamiltonian describing electrons coupled to phonon degrees of freedom. The bosonic part of the action helps control the variation of the field in imaginary time. As a consequence, the iterative conjugate gradient solution of the fermionic action, which depends on the boson coordinates, converges more rapidly than in the case of electron-electron interactions, such as the Hubbard Hamiltonian. Fourier acceleration is shown to be a crucial ingredient in reducing the equilibration and autocorrelation times. After describing and benchmarking the method, we present results for the phase diagram focusing on the range of the electron-phonon interaction. Furthermore, we delineate the regions of charge density wave formation from those in which the fermion density is inhomogeneous, caused by phase separation. We show that the Langevin approach is more efficient than the determinant QMC method for lattice sizes N ≳ 8 × 8 and that it therefore opens a potential path to problems including, for example, charge order in the three-dimensional Holstein model.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2013

A combined DFT and restricted open-shell configuration interaction method including spin-orbit coupling: Application to transition metal L-edge X-ray absorption spectroscopy

A novel restricted-open-shell configuration interaction with singles (ROCIS) approach for the calculation of transition metal L-edge X-ray absorption spectra is introduced. In this method, one first calculates the ground state and a number of excited states of the non-relativistic Hamiltonian. By construction, the total spin is a good quantum number in each of these states. For a ground state with total spin S excited states with spin S Prime = S, S - 1, and S + 1 are constructed. Using Wigner-Eckart algebra, all magnetic sublevels with M{sub S}= S, Horizontal-Ellipsis , -S for each multiplet of spin S are obtained. The spin-orbit operator is represented by a mean-field approximation to the full Breit-Pauli spin-orbit operator and is diagonalized over this N-particle basis. This is equivalent to a quasi-degenerate treatment of the spin-orbit interaction to all orders. Importantly, the excitation space spans all of the molecular multiplets that arise from the atomic Russell-Saunders terms. Hence, the method represents a rigorous first-principles approach to the complicated low-symmetry molecular multiplet problem met in L-edge X-ray absorption spectroscopy. In order to gain computational efficiency, as well as additional accuracy, the excitation space is restricted to single excitations and the configuration interaction matrix is slightly parameterized in order to account for dynamic correlation effects in an average way. To this end, it is advantageous to employ Kohn-Sham rather than Hartree-Fock orbitals thus defining the density functional theory/ROCIS method. However, the method can also be used in an entirely non-empirical fashion. Only three global empirical parameters are introduced and have been determined here for future application of the method to any system containing any transition metal. The three parameters were carefully calibrated using the L-edge X-ray absorption spectroscopy spectra of a test set of coordination complexes containing first row transition metals. These parameters are universal and transferable. Hence, there are no adjustable parameters that are used to fit experimental X-ray absorption spectra. Thus, the new approach classifies as a predictive first-principles method rather than an analysis tool. A series of calculations on transition metal compounds containing Cu, Ti, Fe, and Ni in various oxidation and spin states is investigated and a detailed comparison to experimental data is reported. In most cases, the approach yields good to excellent agreement with experiment. In addition, the origin of the observed spectral features is discussed in terms of the electronic structure of the investigated compounds.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Electronic Origin of T c in Bulk and Monolayer FeSe

FeSe is classed as a Hund’s metal, with a multiplicity of d bands near the Fermi level. Correlations in Hund’s metals mostly originate from the exchange parameter J, which can drive a strong orbital selectivity in the correlations. The Fe-chalcogens are the most strongly correlated of the Fe-based superconductors, with d xy the most correlated orbital. Yet little is understood whether and how such correlations directly affect the superconducting instability in Hund’s systems. By applying a recently developed ab initio theory, we show explicitly the connections between correlations in d xy and the superconducting critical temperature T c . Starting from the ab initio results as a reference, we consider various kinds of excursions in parameter space around the reference to determine what controls T c . We show small excursions in J can cause colossal changes in T c . Additionally we consider changes in hopping by varying the Fe-Se bond length in bulk, in the free standing monolayer M-FeSe, and M-FeSe on a SrTiO 3 substrate (M-FeSe/STO). The twin conditions of proximity of the dxy state to the Fermi energy, and the strength of J emerge as the primary criteria for incoherent spectral response and enhanced single- and two-particle scattering that in turn controls T c . Using c-RPA, we show further that FeSe in monolayer form (M-FeSe) provides a natural mechanism to enhance J. We explain why M-FeSe/STO has a high T c , whereas M-FeSe in isolation should not. Our study opens a paradigm for a unified understanding what controls Tc in bulk, layers, and interfaces of Hund’s metals by hole pocket and electron screening cloud engineering.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2017

Modification of Transition-Metal Redox by Interstitial Water in Hexacyanometalate Electrodes for Sodium-Ion Batteries

A sodium-ion battery (SIB) solution is attractive for grid-scale electrical energy storage. Low-cost hexacyanometalate is a promising electrode material for SIBs because of its easy synthesis and open framework. Most hexacyanometalate-based SIBs work with aqueous electrolyte, and interstitial water in the material has been found to strongly affect the electrochemical profile, but the mechanism remains elusive. Here we provide a comparative study of the transition-metal redox in hexacyanometalate electrodes with and without interstitial water based on soft X-ray absorption spectroscopy and theoretical calculations. We found distinct transition-metal redox sequences in hydrated and anhydrated Na x MnFe(CN) 6 ·zH 2 O. The Fe and Mn redox in hydrated electrodes are separated and are at different potentials, leading to two voltage plateaus. On the contrary, mixed Fe and Mn redox in the same potential range is found in the anhydrated system. This work reveals for the first time how transition-metal redox in batteries is strongly affected by interstitial molecules that are seemingly spectators. The results suggest a fundamental mechanism based on three competing factors that determine the transition-metal redox potentials. Finally, because most hexacyanometalate electrodes contain water, this work directly reveals the mechanism of how interstitial molecules could define the electrochemical profile, especially for electrodes based on transition-metal redox with well-defined spin states.

25 ENERGY STORAGE↗
DOE OSTI2022

Magnon heat transport in a two-dimensional Mott insulator

Whether or not anomalies in the thermal conductivity in insulating cuprates can be attributed to antiferromagnetic order and magnons in a 2D Mott insulator remains an intriguing open question. To shed light on this issue, we investigate the thermal conductivity $\kappa$ and its relationship with the specific heat $c_v$ in the half-filled 2D single-band Hubbard model, using the numerically exact determinant quantum Monte Carlo algorithm and maximum entropy analytic continuation. At low temperatures where the charge degrees of freedom are gapped-out and $c_v$ exhibits a clear magnon peak, we observe that thermal conductivity $\kappa$ also tends to form a peak at similar temperatures. Reducing temperature further produces a sharp upturn in $\kappa$, associated with an increasing mean-free path. Here we identify this as the high-temperature side of the anomalous peak in insulating cuprates, where the mean-free path eventually is cut-off by other scattering effects, including phonons, disorder, and physical size. Different scattering effects in our model are identified and analyzed in the thermal diffusivity.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2014

The role of non-equilibrium fluxes in the relaxation processes of the linear chemical master equation

We propose a non-equilibrium thermodynamical description in terms of the Chemical Master Equation (CME) to characterize the dynamics of a chemical cycle chain reaction among m different species. These systems can be closed or open for energy and molecules exchange with the environment, which determines how they relax to the stationary state. Closed systems reach an equilibrium state (characterized by the detailed balance condition (D.B.)), while open systems will reach a non-equilibrium steady state (NESS). The principal difference between D.B. and NESS is due to the presence of chemical fluxes. In the D.B. condition the fluxes are absent while for the NESS case, the chemical fluxes are necessary for the state maintaining. All the biological systems are characterized by their “far from equilibrium behavior,” hence the NESS is a good candidate for a realistic description of the dynamical and thermodynamical properties of living organisms. In this work we consider a CME written in terms of a discrete Kolmogorov forward equation, which lead us to write explicitly the non-equilibrium chemical fluxes. For systems in NESS, we show that there is a non-conservative “external vector field” whose is linearly proportional to the chemical fluxes. We also demonstrate that the modulation of these external fields does not change their stationary distributions, which ensure us to study the same system and outline the differences in the system's behavior when it switches from the D.B. regime to NESS. We were interested to see how the non-equilibrium fluxes influence the relaxation process during the reaching of the stationary distribution. By performing analytical and numerical analysis, our central result is that the presence of the non-equilibrium chemical fluxes reduces the characteristic relaxation time with respect to the D.B. condition. Within a biochemical and biological perspective, this result can be related to the “plasticity property” of biological systems and to their capabilities to switch from one state to another as is observed during synaptic plasticity, cell fate determination, and differentiation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2018

Self-organized Segregation on the Grid

We consider an agent-based model with exponentially distributed waiting times in which two types of agents interact locally over a graph, and based on this interaction and on the value of a common intolerance threshold τ, decide whether to change their types. This is equivalent to a zero-temperature ising model with Glauber dynamics, an asynchronous cellular automaton with extended Moore neighborhoods, or a Schelling model of self-organized segregation in an open system, and has applications in the analysis of social and biological networks, and spin glasses systems. Some rigorous results were recently obtained in the theoretical computer science literature, and this work provides several extensions. We enlarge the intolerance interval leading to the expected formation of large segregated regions of agents of a single type from the known size ϵ>0 to size ≈0.134. Namely, we show that for 0.433<τ<1/2 (and by symmetry 1/2<τ<0.567), the expected size of the largest segregated region containing an arbitrary agent is exponential in the size of the neighborhood. We further extend the interval leading to expected large segregated regions to size ≈0.312 considering “almost segregated” regions, namely regions where the ratio of the number of agents of one type and the number of agents of the other type vanishes quickly as the size of the neighborhood grows. In this case, we show that for 0.344<τ≤0.433 (and by symmetry for 0.567≤τ<0.656) the expected size of the largest almost segregated region containing an arbitrary agent is exponential in the size of the neighborhood. This behavior is reminiscent of supercritical percolation, where small clusters of empty sites can be observed within any sufficiently large region of the occupied percolation cluster. The exponential bounds that we provide also imply that complete segregation, where agents of a single type cover the whole grid, does not occur with high probability for p=1/2 and the range of intolerance considered.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2018

Two-step perfectly matched layer for arbitrary-order pseudo-spectral analytical time-domain methods

Numerical simulation of an electrodynamic system in empty space requires the implementation of open boundary conditions (BC) to terminate the solution of Maxwell’s equations on the boundaries of the computational domain. The Perfectly Matched Layer (PML) has become the method of choice for open BC with wave equations, as it is straightforward and relatively easy to implement, and offers very efficient and user-adjustable absorption rates. PMLs are most often employed with the Finite-Difference Time-Domain (FDTD) algorithm, which in its most common implementation offers second-order accuracy in space and time on Cartesian grids. Yet, simulations (including some class of electromagnetic Particle-In-Cell simulations) that require higher precision may resort to higher-order Maxwell solvers employing extended finite-difference stencils, or even to pseudo-spectral Maxwell solvers, for which a general, versatile and efficient formulation of the PML has been missing so far. Here we propose a novel “two-step” formulation of the PML that is simple, very versatile and can be used as is with any Maxwell solver. In particular, it is applicable to a large class of Maxwell solvers including the arbitrary-order Pseudo-Spectral Analytical Time-Domain (PSATD) solver, which offers arbitrarily low numerical dispersion when increasing solver order and becomes dispersion-free at infinite order. Analysis and numerical simulations demonstrate that the new formulation is as efficient as the standard PML formulation, both for the FDTD and the PSATD implementations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2013

Power law inflation with electromagnetism

We generalize Ringström’s global future causal stability results (Ringström 2009) [11] for certain expanding cosmological solutions of the Einstein-scalar field equations to solutions of the Einstein–Maxwell-scalar field system. In particular, after noting that the power law inflationary spacetimes (M{sup n+1},g{sup -hat}, ϕ{sup -hat}) considered by Ringström (2009) in [11] are solutions of the Einstein–Maxwell-scalar field system (with exponential potential) as well as of the Einstein-scalar field system (with the same exponential potential), we consider (nonlinear) perturbations of initial data sets of these spacetimes which include electromagnetic perturbations as well as gravitational and scalar perturbations. We show that if (as in Ringström (2009) [11]) we focus on pairs of relatively scaled open sets U{sub R{sub 0}}⊂U{sub 4R{sub 0}} on an initial slice of (M{sup n+1},g{sup -hat}), and if we choose a set of perturbed data which on U{sub 4R{sub 0}} is sufficiently close to that of (M{sup n+1},g{sup -hat},ϕ{sup -hat}, A{sup -hat} = 0), then in the maximal globally hyperbolic spacetime development (M{sup n+1},g,ϕ,A) of this data via the Einstein–Maxwell-scalar field equations, all causal geodesics emanating from U{sub R{sub 0}} are future complete (just as in (M{sup n+1},g{sup -hat})). We also verify that, in a certain sense, the future asymptotic behavior of the fields in the spacetime developments of the perturbed data sets does not differ significantly from the future asymptotic behavior of (M{sup n+1},g{sup -hat}, ϕ{sup -hat}, A{sup -hat} = 0). -- Highlights: •We prove stability of expanding solutions of the Einstein–Maxwell-scalar field equations. •All nearby solutions are geodesically complete. •The topology of the initial slice is irrelevant to our stability results.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2014

Towards identifying the mechanisms underlying field-aligned edge-loss of HHFW power on NSTX

Fast-wave heating will be a major heating scheme on ITER, as it can heat ions directly and is relatively unaffected by the large machine size unlike neutral beams. However, fast-wave interactions with the plasma edge can lead to deleterious effects such as, in the case of the high-harmonic fast-wave (HHFW) system on NSTX, large losses of fast-wave power in the scrape off layer (SOL) under certain conditions. In such scenarios, a large fraction of the lost HHFW power is deposited on the upper and lower divertors in bright spiral shapes. The responsible mechanism(s) has not yet been identified but may include fast-wave propagation in the scrape off layer, parametric decay instability, and RF currents driven by the antenna reactive fields. Understanding and mitigating these losses is important not only for improving the heating and current-drive on NSTX-Upgrade but also for understanding fast-wave propagation across the SOL in any fast-wave system. This talk summarizes experimental results demonstrating that the flow of lost HHFW power to the divertor regions largely follows the open SOL magnetic field lines. This lost power flux is relatively large close to both the antenna and the last closed flux surface with a reduced level in between, so the loss mechanism cannot be localized to the antenna. At the same time, significant losses also occur along field lines connected to the inboard edge of the bottom antenna plate. The power lost within the spirals is roughly estimated, showing that these field-aligned losses to the divertor are significant but may not account for the total HHFW loss. To elucidate the role of the onset layer for perpendicular fast-wave propagation with regards to fast-wave propagation in the SOL, a cylindrical cold-plasma model is being developed. This model, in addition to advanced RF codes such as TORIC and AORSA, is aimed at identifying the underlying mechanism(s) behind these SOL losses, to minimize their effects in NSTX-U, and to predict their importance in ITER.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Single-Crystalline InGaAs/InP Dense Micro-Pillar Forest on Poly-Silicon Substrates for Low-Cost High-Efficiency Solar Cells

The ultimate goal of this project is to develop a photovoltaic system high conversion efficiency (>20%) using high quality III-V compound-based three-dimensional micro-structures on silicon and poly-silicon. Such a PV-system could be of very low cost due to minimum usages of III-V materials. This project will address the barriers that currently hamper the performance of solar cells based on three-dimensional micro-structures. To accomplish this goal the project is divided into 4 tasks, each dealing with a different aspect of the project: materials quality, micropillar growth control, light management, and pillar based solar cells. Materials Quality: the internal quantum efficiency (IQE) - by which is meant here the internal fluorescence yield - of the micro-pillars has to be increased. We aim at achieving an IQE of 45% by the end of the first year. By the end of the second year there will be a go-no-go milestone of 65% IQE. By the end of year 3 and 4 we aim to achieve 75% and 90% IQE, respectively. Micropillar growth control: dense forests of micropillars with high fill ratios need to be grown. Pillars within forests should show minimum variations in size. We aim at achieving fill ratios of 2%, 10%, >15%, >20% in years 1, 2, 3, and 4, respectively. Variations in dimension should be minimized by site-controlled growth of pillars. By the end of year 1 we will aim at achieving site-controlled growth with > 15% yield. By end of year 2 the variation of critical pillar dimensions should be less than 25%. Light management: high light absorption in the spectral range of the sun has been to be demonstrated for the micropillar forests. By the end of year 1 we will employ FDTD simulation techniques to demonstrate that pillar forests with fill ratios <20% can achieve 99% light absorption. By end of year 2 our original goal was to demonstrate >85% absorption. By end of year 3 > 90% absorption should be demonstrated. Pillar based solar cells: devices will be studied to explore ways to achieve high open-circuit voltages which will lead to high efficiency micropillar-based solar cells. We will start on single pillar devices and the findings in these studies should pave the way for devices based on forests/ arrays of pillars. By the end of the second year we aim to demonstrate a single pillar device with an open-circuit voltage of 0.7 V, as well as a pillar-forest based device with 8% conversion efficiency. By the end of year 3 these numbers should be improved to 0.9 V open-circuit voltage for single pillar devices and >15% efficiency for forest/array-based devices. We will aim to realize a device with 20% efficiency by the end of the project period.

Chang-Hasnain, Constance↗
DOE OSTI2017

Effects of heat exchanger tubes on hydrodynamics and CO 2 capture of a sorbent-based fluidized bed reactor

In virtual design and scale up of pilot-scale carbon capture systems, the coupled reactive multiphase flow problem must be solved to predict the adsorber’s performance and capture efficiency under various operation conditions. This paper focuses on the detailed computational fluid dynamics (CFD) modeling of a pilot-scale fluidized bed adsorber equipped with vertical cooling tubes. Multiphase Flow with Interphase eXchanges (MFiX), an open-source multiphase flow CFD solver, is used for the simulations with custom code to simulate the chemical reactions and filtered models to capture the effect of the unresolved details in the coarser mesh for simulations with reasonable simulations and manageable computational effort. Previously developed two filtered models for horizontal cylinder drag, heat transfer, and reaction kinetics have been modified to derive the 2D filtered models representing vertical cylinders in the coarse-grid CFD simulations. The effects of the heat exchanger configurations (i.e., horizontal or vertical) on the adsorber’s hydrodynamics and CO2 capture performance are then examined. The simulation result subsequently is compared and contrasted with another predicted by a one-dimensional three-region process model.

bubbling bed↗
DOE OSTI2024

Dynamic embedding of effective harmonic normal mode vibrations in all-atomistic energy gap fluctuations: Case study of light harvesting 2 complex

Environmental effects in excitation energy transfer have mostly been modeled by baths of harmonic oscillators, but to what extent such modeling provides a reliable description of actual interactions between molecular systems and environments remains an open issue. Here, we address this issue by investigating fluctuations in the excitation energies of the light harvesting 2 complex using a realistic all-atomistic simulation of the potential energy surface. Our analyses reveal that molecular motions exhibit significant anharmonic features, even for underdamped intramolecular vibrations. In particular, we find that the anharmonicity contributes to the broadening of spectral densities and substantial overlaps between neighboring peaks, which complicates the meaning of mode frequencies constituting a bath model. Thus, we develop a strategy to construct a minimally underdamped harmonic bath that has a clear connection to all-atomistic dynamics by utilizing actual normal modes of molecules but optimizing their frequencies such that the resulting bath model can best reproduce the all-atomistic simulation results. By subtracting the underdamped contribution from the entire fluctuations, we also show that identifying a residual spectral density representing all other contributions with overdamped behavior is possible. We find that this can be fitted well with a well-established analytic form of a spectral density function or, alternatively, modeled as explicit time dependent fluctuations with muti-exponential or power law type correlation functions. We provide an assessment and the implications of these possibilities. The approach presented here can also serve as a general strategy to construct a simplified bath model that can effectively represent the underlying all-atomistic bath dynamics.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Spacecraft Geometry Effects on Kinetic Impactor Missions

Abstract The DART (Double Asteroid Redirection Test) mission will impact a spacecraft on the secondary (Dimorphos) of the binary asteroid system Didymos in 2022 September, with the goal of altering the orbital period of Dimorphos about Didymos sufficiently to be observed from ground-based observations. Numerical impact modeling is a crucial component in understanding the outcome of the DART experiment, and while many have investigated the effects of target properties, such as material strength and porosity (which remain unknown), an often overlooked factor is the importance of accurately representing the spacecraft itself in such models. Most impact modeling to date has considered simple impactor geometries such as a solid uniform sphere, but in reality the spacecraft is a complex shape full of different components, open spaces, and thin walled structures. At a minimum, a simple solid representation underestimates the surface area of the impact: for a small body such as Dimorphos (approximately 160 m in diameter), the difference between a spacecraft spanning 20 m (including solar arrays) impacting and a sub-1 m idealized shape may be important. In this paper, we compare models impacting high-fidelity models of the spacecraft based on the CAD geometry with various simplified impactors, in order to assess the potential importance of this effect. We find that the difference between the simplest impactor geometries (such as a uniform sphere) and the real spacecraft is measurable, and has an interesting dependence on the material properties of the asteroid itself.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Magic-Angle Twisted Bilayer Graphene as a Topological Heavy Fermion Problem

Magic-angle (θ=1.05°) twisted bilayer graphene (MATBG) has shown two seemingly contradictory characters: the localization and quantum-dot-like behavior in STM experiments, and delocalization in transport experiments. We construct a model, which naturally captures the two aspects, from the Bistritzer-MacDonald (BM) model in a first principle spirit. A set of local flat-band orbitals (f) centered at the AA-stacking regions are responsible to the localization. A set of extended topological semimetallic conduction bands (c), which are at small energetic separation from the local orbitals, are responsible to the delocalization and transport. The topological flat bands of the BM model appear as a result of the hybridization of f and c electrons. This model then provides a new perspective for the strong correlation physics, which is now described as strongly correlated f electrons coupled to nearly free c electrons—we hence name our model as the topological heavy fermion model. Using this model, we obtain the U(4) and U(4)×U(4) symmetries of Refs. [1–5] as well as the correlated insulator phases and their energies. Simple rules for the ground states and their Chern numbers are derived. Moreover, features such as the large dispersion of the charge ±1 excitations [2,6,7], and the minima of the charge gap at the Γ M point can now, for the first time, be understood both qualitatively and quantitatively in a simple physical picture. Our mapping opens the prospect of using heavy-fermion physics machinery to the superconducting physics of MATBG.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2026

Design Rules for Open-Shell Molecular Wires: Insights from Correlated Many-Body Transport

Open-shell π-conjugated systems have emerged as promising molecular wires that can sustain unusually high low-bias conductance over tens of nanometers. However, predictive design rules remain limited because near-degeneracy and spin polarization give rise to correlated electronic states that are not reliably captured by standard density functional theory and Landauer transport descriptions. Here, we establish physically transparent design rules for long-range, weakly length-dependent charge transport based on the interplay between bond-length alternation and electron–electron interactions. Using a fully correlated transport framework that combines density matrix renormalization group calculations with nonequilibrium Green’s function embedding, we show that weak dimerization maintains spatially extended edge states, while moderate on-site interaction stabilizes open-shell character without excessively separating transport-relevant resonances. This balance yields zero-bias transmission that is largely insensitive to molecular length, providing a route to connect chemically tunable structure to correlated transport beyond mean-field descriptions.

Charge transport↗
DOE OSTI2017

Modeling of wave processes in blocky media with porous and fluid-saturated interlayers

The wave processes in blocky media are analyzed by applying different mathematical models, wherein the elastic blocks interact with each other via pliant interlayers with the complex mechanical properties. Four versions of constitutive equations are considered. In the first version, an elastic interaction between the blocks is simulated within the framework of linear elasticity theory, and the model of elastic–plastic interlayers is constructed to take into account the appearance of irreversible deformation of interlayers at short time intervals. In the second one, the effects of viscoelastic shear in the interblock interlayers are taken into the consideration using the Poynting–Thomson rheological scheme. In the third option, the model of an elastic porous material is used in the interlayers, where the pores collapse if an abrupt compressive stress is applied. In the fourth case, the model of a fluid-saturated material with open pores is examined based on Biot's equations. The collapse of pores is modeled by the generalized rheological approach, wherein the mechanical properties of a material are simulated using four rheological elements. Three of them are the traditional elastic, viscous and plastic elements, the fourth element is the so-called rigid contact, which is used to describe the behavior of materials with the different resistance to tension and compression. It was shown that the thermodynamically consistent model is provided, which means that the energy balance equation is fulfilled for an entire blocky structure, where the kinetic and potential energy of the system is the sum of the kinetic and potential energies of the blocks and interlayers. Under numerical implementation of the interlayers models, the dissipationless finite difference Ivanov's method was used. The splitting method by spatial variables in the combination with the Godunov gap decay scheme was applied in the blocks. As a result, robust and stable computational algorithms are built and tested. Using MPI technology, the parallel software was designed for the modeling of wave processes in 2D setting. The numerical results are presented, discussed and future studies are outlined. - Highlights: • We model the propagation of waves in a blocky medium with compliant interlayers. • Models of porous and fluid saturated interlayers are proposed using Biot's approach. • A parallel computational algorithm is developed for supercomputers of cluster type. • Problems of instant rotation of central and boundary blocks are solved numerically. • Wave motions are compared in the cases of elastic and fluid-saturated interlayers.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Robust Spin-Moiré-Superlattice-Driven Gap Opening in EuAg 4 Sb 2 under in-Plane Magnetic Field

Moiré superlattices introduce new length and energy scales, enabling discoveries, such as unconventional superconductivity, in 2D systems. This concept has recently been extended to bulk materials with multiple-q spin textures, opening exciting opportunities for spin moiré physics. A notable example is EuAg 4 S b2 , where a spin moiré superlattice (SMS), manifested as a double-q spin modulation, induces a superzone gap opening. In this work, we investigate the tunability of this gap by an in-plane magnetic field in EuAg 4 Sb 2 by using neutron scattering, magnetization, and transport measurements. We reveal a cycloidal ground state and multiple spin-reoriented phases induced by the in-plane field, highlighting the critical role of in-plane components in driving magnetic transitions. Moreover, we demonstrate that a robust gap opening persists in the double-q phase, regardless of in-plane field orientation. Model calculations attribute this robustness to the stability of the SMS under tilted fields. Here, our results establish EuAg 4 Sb 2 as a tunable platform for exploring the spin-texture-driven superzone gap opening in electronic states.

36 MATERIALS SCIENCE↗