Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “continuum modeling”

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 325 records · Page 18
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Sound speed measurements in shock-compressed cemented Tungsten carbides at pressures up to 100 GPa

To gain insights into thermodynamic states attained during shock compression of cemented tungsten carbide with 3.7 wt.% cobalt binder, we present results of longitudinal sound (release wave) speed measurements and their analysis at peak stresses up to 100 GPa (volumetric compression ratio ~ 15%). The sound speeds are determined using front-surface impact and release-wave overtake plate impact experimental configurations using laser interferometry. The measured sound speed data along with estimates for bulk sound speeds obtained using the fourth-order Birch-Murnaghan EoS and thermodynamics are used to determine the longitudinal moduli and shear moduli of shocked tungsten carbide at the various peak compression states attained in the experiments. Here, the longitudinal sound speeds were found to increase linearly with volume compression ratio from 6.97 ± 0.010 km/s at ambient conditions to 8.26 ± 0.156 km/s at a volume compression ratio of ~ 15%. The corresponding longitudinal elastic moduli also increase nearly linearly with the volume compression ratio but remain consistently lower than their theoretical predictions based on continuum models with no damage. Also, the sensitivity of shear moduli to pressure, as predicted by the Steinberg-Guinan model, is reduced substantially and the shear moduli of cemented WC with 3.7 wt.% Co remains nearly constant at ~ 310 GPa at the various peak compression stress states investigated in the present study.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Low-frequency Raman spectrum of 2D layered perovskites: Local atomistic motion or superlattice modes?

We report the low-frequency Raman spectrum (ω = 10 cm -1 –150 cm -1 ) of a wide variety of alkylammonium iodide based 2D lead halide perovskites (2D LHPs) as a function of A-site cation (MA = methylammonium and FA = formamidinium), octahedral layer thickness (n = 2–4), organic spacer chain length (butyl-, pentyl-, hexyl-), and sample temperature (T = 77 K–293 K). Using density functional theory calculations under the harmonic approximation for n = 2 BA:MAPbI, we assign several longitudinal/transverse optical phonon modes between 30 cm -1 and 100 cm -1 , the eigen displacements of which are analogous to that observed previously for octahedral twists/distortions in bulk MAPbI. Additionally, we propose an alternative assignment for low-frequency modes below this band (<30 cm -1 ) as zone-folded longitudinal acoustic phonons corresponding to the periodicity of the entire layered structure. Furthermore, we compare measured spectra to predictions of the Rytov elastic continuum model for zone-folded dispersion in layered structures. Our results are consistent across the various 2D LHPs studied herein, with energetic shifts of optical phonons corresponding to microscopic structural differences between materials and energetic shifts of acoustic phonons according to changes in the periodicity and elastic properties of the perovskite/organic subphases. This study highlights the importance of both the local atomic order and the superlattice structure on the vibrational properties of layered 2D materials.

2D materials↗
DOE OSTI2021

Quadrupole-mediated dielectric response and the charge-asymmetric solvation of ions in water

Treating water as a linearly responding dielectric continuum on molecular length scales allows very simple estimates of the solvation structure and thermodynamics for charged and polar solutes. While this approach can successfully account for basic length and energy scales of ion solvation, computer simulations indicate not only its quantitative inaccuracies but also its inability to capture some basic and important aspects of microscopic polarization response. Here, we consider one such shortcoming, a failure to distinguish the solvation thermodynamics of cations from that of otherwise-identical anions, and we pursue a simple, physically inspired modification of the dielectric continuum model to address it. The adaptation is motivated by analyzing the orientational response of an isolated water molecule whose dipole is rigidly constrained. Its free energy suggests a Hamiltonian for dipole fluctuations that accounts implicitly for the influence of higher-order multipole moments while respecting constraints of molecular geometry. Finally, we propose a field theory with the suggested form, whose nonlinear response breaks the charge symmetry of ion solvation. An approximate variational solution of this theory, with a single adjustable parameter, yields solvation free energies that agree closely with simulation results over a considerable range of solute size and charge.

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

Thermal conductivity tensor of β -1,3,5,7-tetranitro-1,3,5,7-tetrazoctane ( β -HMX) as a function of pressure and temperature

Here, we have used reverse non-equilibrium molecular dynamics (RNEMD) simulations to determine the full thermal conductivity tensor for the monoclinic high explosive crystal β-1,3,5,7-tetranitro-1,3,5,7-tetrazoctane (β-HMX). In order to do so for the monoclinic crystal, four directions for heat propagation are used. Effects of the temperature and pressure are investigated between 200 and 500 K and 0 and 5 GPa, respectively, which approximately covers the range where the β polymorph is stable. Simulations are carried out with the Smith–Bharadwaj non-reactive empirical potential [Smith and Bharadwaj, J. Phys. Chem. B 103, 3570 (1999)], which is known to reproduce well the thermo-elastic properties of HMX. Our results indicate that the thermal conductivity, κ, is highly anisotropic, with 36% difference between the two extreme values at 300 K and 0 GPa. A simple function is used to interpolate κ in the pressure-temperature regime considered in this study, which can be used in continuum models. The results from RNEMD simulations compare well with available experimental results from the literature and allow the determination of κ for any direction and temperature and pressure within or around the fitting interval.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Mechanically activated and deactivated ion transport across nanopores with heterogeneous surface charge distributions

To mimic the intricate and adaptive functionalities of biological ion channels, electrohydrodynamic ion transport has been studied extensively, albeit mostly, across uniformly charged nanochannels. Here, we analyze the ion transport under coupled electric field and pressure across heterogeneously charged nanopores with oppositely charged sections on their lateral surface. We only consider such pores with symmetric hourglass-like and cylindrical shapes to focus on the effects of the non-uniform surface charge distribution. Finite-element simulations of a continuum model demonstrate that a pressure applied in either direction of the pore-axis equally suppresses or amplifies the ionic conductance, depending on the electric field polarity, by distorting the quasi-static distribution of ions in the pore. The resulting anomalous mechanical deactivation and activation of ionic current under opposite voltage biases exhibit the functional modularity of our setup, while their intensities are highly tunable, substantially greater than those of analogous behaviors in other nanochannels, and fundamentally correlated to ionic current rectification (ICR) in our pores. A detailed study of ICR subsequently reveals its counterintuitive non-monotonous variations, in the pores, with the magnitude of applied voltage and the pore length, that can help optimize their diode-like behavior. We further illustrate that while the hourglass-shaped nanopores yield the more efficient mechanical suppressors of ion transport, their cylindrical analogs are the superior rectifiers and mechanical amplifiers of ion conduction. Therefore, this article provides a blueprint for the strategic design of nanofluidic circuits to attain a robust, modular, and tunable control of ion transport under external electrical and mechanical stimuli.

Physics↗
DOE OSTI2026

Finite-element-based simulations of electrodes for CO 2 cascade reduction reactions

The multielectron reduction of CO 2 to liquid fuels could be a path to scalable energy storage, but reaching this goal requires major advances in catalysis and systems engineering. Cascade catalysis, which couples sequential reactions without isolating intermediates, has emerged as a promising route to enhance selectivity and efficiency in CO 2 reduction (CO 2 R). In this review, we examine how finite-element-based simulations of continuum model [finite element method (FEM)] approaches are being used to analyze and guide CO 2 R cascade systems. We first outline the fundamentals of cascade catalysis and recent advances in catalytic materials (metallic, molecular, and hybrid architectures). We then focus on FEM developments at the electrode and device scales, emphasizing how these models capture transport phenomena, local microenvironments, and geometry-dependent effects. To clarify design principles, we present case studies of cascade electrodes organized in systems without and with integrated semiconductors. We further emphasize the integration of FEM with multiscale frameworks (density functional theory, molecular dynamics, kinetic Monte Carlo) and its role in bridging atomic-level insights with device-level performance. Finally, we identify current limitations and future prospects, including improved boundary conditions, coupling with operando experiments, and machine learning-accelerated model development. Together, these insights provide design principles for next-generation CO 2 R cascade systems for efficient solar fuel production.

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

Autonomous phototaxis of hydrogel swimmers

The design of synthetic soft matter capable of emulating the complex behaviors of living organisms, such as sensing and adapting to their environment, remains an important challenge in developing biomimetic materials. Functionalized hydrogels are ideal candidates for such materials since they are highly responsive to their environment and can be operated in water. In this work, we investigate a hybrid bonding hydrogel composed of peptide amphiphile supramolecular nanofibers covalently attached to a photoresponsive network, in which high-aspect-ratio ferromagnetic nanowires are aligned along the length of the sample, designed to swim under oscillating magnetic fields. This hybrid hydrogel swimmer can autonomously swim toward a light source by utilizing photoinduced interactions between supramolecular and covalent networks reminiscent of phototactic swimming in living systems. Using a combination of experimental techniques and a continuum model incorporating photochemistry, magnetoelasticity, and hydrodynamics, we explain the swimming mechanism and predict phototactic behavior. Our work highlights the potential role of hybrid bonding polymers, which leverage the interplay between supramolecular assemblies and covalent networks. We demonstrate how these polymers can be tailored to react dynamically to their environment, paving the way for developing intelligent and autonomous robotic systems.

Science & Technology - Other Topics↗
DOE OSTI2021

Electric-field-tunable electronic nematic order in twisted double-bilayer graphene

Graphene-based moiré systems have attracted considerable interest in recent years as they display a remarkable variety of correlated phenomena. Besides insulating and superconducting phases in the vicinity of integer fillings of the moiré unit cell, there is growing evidence for electronic nematic order both in twisted bilayer graphene and twisted double-bilayer graphene (tDBG), as signaled by the spontaneous breaking of the threefold rotational symmetry of the moiré superlattices. Here, we combine symmetry-based analysis with a microscopic continuum model to investigate the structure of the nematic phase of tDBG and its experimental manifestations. First, we perform a detailed comparison between the theoretically calculated local density of states and recent scanning tunneling microscopy data (arXiv:2009.11645) to resolve the internal structure of the nematic order parameter in terms of the layer, sublattice, spin, and valley degrees of freedom. We find strong evidence that the dominant contribution to the nematic order parameter comes from states at the moiré scale rather than at the microscopic scale of the individual graphene layers, which demonstrates the key role played by the moiré degrees of freedom and confirms the correlated nature of the nematic phase in tDBG. Secondly, our analysis reveals an unprecedented tunability of the orientation of the nematic director in tDBG by an externally applied electric field, allowing the director to rotate away from high-symmetry crystalline directions. Here, we compute the expected fingerprints of this rotation in both STM and transport experiments, providing feasible ways to probe it. Rooted in the strong sensitivity of the flat bands of tDBG to the displacement field, this effect opens an interesting route to the electrostatic control of electronic nematicity in moiré systems.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Electromagnetic response and emergent topological orders in transition metal dichalcogenide MoTe 2 bilayers

Twisted bilayer transition metal dichalcogenides, such as MoTe 2 , provide a versatile platform for exploring correlated topological phases. This work investigates the interplay between perpendicular magnetic and electric fields in tuning the electronic structure and emergent topological orders of twisted bilayer MoTe 2 (t-MoTe 2 ) across two distinct regimes: a low-twist-angle phase (θ ≈ 2.1°) hosting multiple Chern bands of identical Chern numbers per valley, and a higher-angle phase (θ ≈ 3.89°) featuring Haldane-like bands with opposite Chern numbers. Using a continuum model incorporating moiré potentials up to second harmonics, we compute the Hofstadter fractal spectra under applied fields, revealing Landau fan structures and magnetic-flux-dependent band topology. These fractal spectra are useful in studying emergent topological orders in terms of the composite fermion picture, where the statistical Chern-Simons flux is approximated as a uniform gauge field. We demonstrate that the system hosts both Jain-sequence fractional Chern insulators (FCIs) and non-Jain fractal FCIs with higher Chern numbers. As a result, the electric field suppresses composite fermion gaps and induces topological quantum phase transitions. Furthermore, our analysis extends to valley-contrasting flux attachment, proposing pathways to describe fractional quantum spin Hall states.

Anyons↗
DOE OSTI2020

Correlated insulating phases of twisted bilayer graphene at commensurate filling fractions: A Hartree-Fock study

Motivated by the recently observed insulating states in twisted bilayer graphene, we study the nature of the correlated insulating phases of the twisted bilayer graphene at commensurate filling fractions. Iin this work, we use the continuum model and project the Coulomb interaction onto the flat bands to study the ground states by using a Hartree-Fock approximation. In the absence of the hexagonal boron nitride substrate, the ground states are the intervalley coherence states at charge neutrality (filling $\textit{ν}$ = 0 , or four electrons per moiré cell) and at $\textit{ν}$ = –1/4 and –1/2 (three and two electrons per cell, respectively) and the $C_2\mathcal{T}$ symmetry-broken state at $\textit{ν}$ = –3/4 (one electron per cell). The hexagonal boron nitride substrate drives the ground states at all $\textit{ν}$ into $C_2\mathcal{T}$ symmetry broken-states. Our results provide good reference points for further study of the rich correlated physics in the twisted bilayer graphene.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Moiré quantum chemistry: Charge transfer in transition metal dichalcogenide superlattices

Transition metal dichalcogenide (TMD) bilayers have recently emerged as a robust and tunable moiré system for studying and designing correlated electron physics. In this Rapid Communication, by combining a large-scale first-principles calculation and continuum model approach, we provide an electronic structure theory that maps long-period TMD heterobilayer superlattices onto diatomic crystals with cations and anions. We find that the interplay between the moiré potential and Coulomb interaction leads to filling-dependent charge transfer between different moiré superlattice regions. Further, we show that the insulating state at half filling found in recent experiments on WSe 2 /WS 2 is a charge-transfer insulator rather than a Mott-Hubbard insulator. Our work reveals the richness of simplicity in moiré quantum chemistry.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Efficient simulation of moiré materials using the density matrix renormalization group

Here we present an infinite density-matrix renormalization group (DMRG) study of an interacting continuum model of twisted bilayer graphene (tBLG) near the magic angle. Because of the long-range Coulomb interaction and the large number of orbital degrees of freedom, tBLG is difficult to study with standard DMRG techniques - even constructing and storing the Hamiltonian already poses a major challenge. To overcome these difficulties, we use a recently developed compression procedure to obtain a matrix product operator representation of the interacting tBLG Hamiltonian which we show is both efficient and accurate even when including the spin, valley, and orbital degrees of freedom. To benchmark our approach, we focus mainly on the spinless, single-valley version of the problem where, at half filling, we find that the ground state is a nematic semimetal. Remarkably, we find that the ground state is essentially a k-space Slater determinant, so that Hartree-Fock and DMRG give virtually identical results for this problem. Our results show that the effects of long-range interactions in magic angle graphene can be efficiently simulated with DMRG and open up a new route for numerically studying strong correlation physics in spinful, two-valley tBLG, and other moiré materials in future work.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Electronic structures, charge transfer, and charge order in twisted transition metal dichalcogenide bilayers

Moiré superlattices of transition metal dichalcogenide (TMD) bilayers have been shown to host correlated electronic states, which arise from the interplay of long wavelength moiré potential and long-range Coulomb interaction. Here, we theoretically investigate structural relaxation and single-particle electronic structure of twisted TMD homobilayer. From the large-scale density functional theory calculation and continuum model with layer degrees of freedom, we find that the out-of-plane gating field creates a tunable charge transfer gap at the Dirac point between the first and second moiré valence bands. We further study the charge orders at the fractional band fillings. In the flat band limit, we find from Monte Carlo simulations a series of charge-ordered insulating states at various fillings $n = \frac{1}{4}, \frac{1}{3}, \frac{1}{2}, \frac{2}{3}, 1$. Further, we predict that gating field induces a phase transition between different electron crystals at fixed filling $n = \frac{1}{2}$ or $\frac{2}{3}$. At half-filling $\textit{n} = 1$, the ground state is a Mott insulator with electronically driven ferroelectricity. This work demonstrates that TMD homobilayer provides a powerful platform for the investigation of tunable charge transfer insulator and charge orders.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Mirror symmetry breaking and lateral stacking shifts in twisted trilayer graphene

Here, we construct a continuum model of twisted trilayer graphene using ab initio density-functional-theory calculations, and we apply it to address twisted trilayer electronic structure. Our model accounts for moiré variation in site energies, hopping between outside layers, and hopping within layers. We focus on the role of a mirror symmetry present in ABA graphene trilayers with a middle layer twist. The mirror symmetry is lost intentionally when a displacement field is applied between layers, and unintentionally when the top layer is shifted laterally relative to the bottom layer. We calculate two band-structure characteristics that are directly relevant to transport measurements, namely the Drude weight and the weak-field Hall conductivity. We relate them via the Hall density to assess the influence of the accidental lateral stacking shifts currently present in all experimental devices on electronic properties. Finally, we comment on the role of the possible importance of accidental lateral stacking shifts for superconductivity in twisted trilayers.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Electric field tunable layer polarization in graphene/boron-nitride twisted quadrilayer superlattices

The recently observed unconventional ferroelectricity in AB bilayer graphene sandwiched by hexagonal boron nitride (hBN) presents a new platform to manipulate correlated phases in multilayered van der Waals heterostructures [Zheng et al., Nature 588, 71 (2020).]. Here, we present a low-energy continuum model for AB bilayer graphene encapsulated by the top and bottom layers of either hBN or graphene, with two independent twist angles. For the graphene/hBN heterostructures, we show that twist angle asymmetry leads to a layer polarization of the valence and conduction bands. We also show that an out-of-plane displacement field not only tunes the layer polarization but also flattens the low-energy bands. We extend the model to show that the electronic structures of quadrilayer graphene heterostructure consisting of AB bilayer graphene encapsulated by the top and bottom graphene layers can similarly be tuned by an external electric field.

Zhu, Ziyan↗
DOE OSTI2022

Formation of the skyrmionic polaron by Rashba and Dresselhaus spin-orbit coupling

Skyrmions in reduced dimensions such as thin layers and interfaces are of both fundamental and technological importance. In these systems, itinerant electrons are often present together with the Rashba and Dresselhaus spin-orbit coupling (SOC). Here, in this study, we show that an itinerant electron in the presence of these interactions can nucleate the skyrmion state, even when the standard Dzyaloshinskii-Moriya interaction (DMI) is absent, and the electron can become self-trapped in the skyrmion core, forming the “skyrmionic polaron” (SkP). The formation of the SkP is investigated from a continuum model of the electron, exchange coupled to the lattice spins, by solving the appropriate Euler-Lagrange equations. The skyrmion (antiskyrmion) texture is favored by the Rashba (Dresselhaus) SOC, with the binding energy increasing quadratically with the strength of the interaction. In contrast, if the skyrmion is already formed due to a nonzero DMI, the electron is delocalized and avoids the skyrmion core until the strength of the Rashba or Dresselhaus SOC exceeds a critical value. Below this critical value, the electron is not bound to the skyrmion core, the polaron does not form, and the electron has little effect on the skyrmion state. Our work envisions the possibility of manipulating the skyrmion state in device applications by altering the strength of the Rashba or Dresselhaus interactions, e.g., by an external electric field.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Spin-triplet superconductivity from intervalley Goldstone modes in magic-angle graphene

Here, we consider magic-angle graphene in the doping regime around charge neutrality, and we study the connection between a recently proposed intervalley coherent insulator at zero doping and the neighboring superconducting domes. The magic-angle graphene continuum model has an emergent U(1) valley-charge conservation symmetry and an emergent SU(2) symmetry corresponding to opposite spin rotations in the two valleys. The intervalley coherent insulator spontaneously breaks both these emergent symmetries, and as a result has four Goldstone modes that couple to doped charge carriers. We derive the effective interaction mediated by the Goldstone modes, and we study its role in electron pair formation. The SU(2) Goldstone modes generate a ferromagnetic interaction, which is attractive in spin-triplet pairing channels and repulsive in spin-singlet channels. From a weak-coupling BCS calculation, we find the leading superconducting instability in the p-wave channel.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Time-reversal invariant topological moiré flat band: A platform for the fractional quantum spin Hall effect

Motivated by recent observation of the quantum spin Hall effect in monolayer germanene and twisted bilayer transition-metal-dichalcogenides (TMDs), we study the topological phases of moir twisted bilayers with time-reversal symmetry and spin sz conservation. By using a continuum model description which can be applied to both germanene and TMD bilayers, we show that at small twist angles the emergent moir flat bands can be topologically nontrivial due to inversion symmetry breaking. Each of these flat bands admits a lowest-Landau-level description for each spin projection in the chiral limit and at magic twist angle. Furthermore, this allows for the construction of a many-body Laughlin state with time-reversal symmetry which can be stabilized by a short-range pseudopotential, and therefore serves as an ideal platform for realizing the so-far elusive fractional quantum spin Hall effect with emergent spin-1/2 U(1) symmetry.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗