Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Order parameters”

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 73 records · Page 4

Group-equivariant autoencoder for identifying spontaneously broken symmetries

We introduce the group-equivariant autoencoder (GE autoencoder), a deep neural network (DNN) method that locates phase boundaries by determining which symmetries of the Hamiltonian have spontaneously broken at each temperature. We use group theory to deduce which symmetries of the system remain intact in all phases, and then use this information to constrain the parameters of the GE autoencoder such that the encoder learns an order parameter invariant to these “never-broken” symmetries. This procedure produces a dramatic reduction in the number of free parameters such that the GE-autoencoder size is independent of the system size. We include symmetry regularization terms in the loss function of the GE autoencoder so that the learned order parameter is also equivariant to the remaining symmetries of the system. By examining the group representation by which the learned order parameter transforms, we are then able to extract information about the associated spontaneous symmetry breaking. We test the GE autoencoder on the 2D classical ferromagnetic and antiferromagnetic Ising models, finding that the GE autoencoder (1) accurately determines which symmetries have spontaneously broken at each temperature; (2) estimates the critical temperature in the thermodynamic limit with greater accuracy, robustness, and time efficiency than a symmetry-agnostic baseline autoencoder; and (3) detects the presence of an external symmetry-breaking magnetic field with greater sensitivity than the baseline method. Lastly, we describe various key implementation details, including a quadratic-programming-based method for extracting the critical temperature estimate from trained autoencoders and calculations of the DNN initialization and learning rate settings required for fair model comparisons.

42 ENGINEERING↗

Role of electromagnetic gauge-field fluctuations in the selection between chiral and nematic superconductivity

Motivated by the observation of nematic superconductivity in several systems, we revisit the problem of the leading pairing instability of two-component unconventional superconductors on the triangular lattice—such as (p x , p y )-wave and (d x 2 –y 2 , d xy )-wave superconductors. Such a system has two possible superconducting states: the chiral state (e.g., p + ip or d + id), which breaks time-reversal symmetry, and the nematic state (e.g., p + p or d + d), which breaks the threefold rotational symmetry of the lattice. Weak-coupling calculations generally favor the chiral over the nematic superconducting state, raising the question of what mechanism can stabilize the latter. Here, we show that the electromagnetic field fluctuations can play a crucial role in selecting between these two states. Specifically, we derive and analyze the effective free energy for the two-component superconducting order parameter after integrating out the gauge-field fluctuations, which is formally justified if the spatial order parameter fluctuations can be neglected. A nonanalytic cubic term arises, as in the case of a conventional s -wave superconductor. However, unlike the latter, the cubic term depends on the relative phase and on the relative amplitudes between the two order parameter components, in such a way that it generally favors the nematic state. This result is a direct consequence of the fact that the stiffness of the superconducting order parameter is not isotropic. Competition with the quartic term, which favors the chiral state, leads to a renormalized phase diagram in which the nematic state displaces the chiral state over a wide region in the parameter space. Here, we analyze the stability of the fluctuation-induced nematic phase, generalize our results to tetragonal lattices, and discuss their applicability to candidate nematic superconductors, including twisted bilayer graphene.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Pairing symmetry of twisted bilayer graphene: A phenomenological synthesis

One of the outstanding questions in the study of twisted bilayer graphene—from both experimental and theoretical points of view—is the nature of its superconducting phase. In this paper we perform a comprehensive synthesis of existing experiments, and argue that experimental constraints are strong enough to allow the structure of the superconducting order parameter to be nearly uniquely determined. In particular, we argue that the order parameter is nodal, and is formed from an admixture of spin-singlet and spin-triplet Cooper pairs. This argument is made on phenomenological grounds, without committing to any particular microscopic model of the superconductor. Existing data is insufficient to determine the orbital parity of the order parameter, which could be either p wave or d wave. We propose a way in which the measurement of Andreev edge states can be used to distinguish between the two.

36 MATERIALS SCIENCE↗

Effect of hydride vapor phase epitaxy growth conditions on the degree of atomic ordering in GaInP

We report a systematic study of CuPt-type ordering in hydride vapor phase epitaxy (HVPE)-grown Ga 0.5 In 0.5 P. Selected-area electron diffraction reveals ordering in samples grown on (001) GaAs substrates offcut toward (1$\overline{1}$1)B. The ordering is single-variant, occurring only on (1$\overline{1}$1)B planes and not on ($\overline{1}$11)B. Quantitative analysis of the order parameter by high-resolution x-ray diffraction (HRXRD) indicates that ordering increases with deposition temperature in samples grown at 600–700 °C with a constant gas-phase V/III ratio ~3. Ordering increases with V/III ratio in the range of 1.3–6.7 at a constant deposition temperature of 650°C. Photoluminescence measurements correlate the order parameter with Ga0.5In0.5P bandgap contraction, though the contraction is larger than expected based on the magnitude of order parameters measured by HRXRD. A possible reason for this discrepancy is that the photoluminescence emission occurs in the lower bandgap ordered domains, which are small and evenly dispersed throughout the material. We also show that the degree of ordering decreases with growth rate, disappearing at ~60 µm/h. The HVPE-grown material exhibits a generally weaker ordering than organometallic vapor phase epitaxy (OMVPE)-grown material, likely due to the moderate V/III ratios employed, in contrast to the V/III ratios in the 100s typical of OMVPE. However, the tendency for ordering to increase with V/III ratio suggests that the same dimer-induced stress mechanism used to explain the occurrence of ordering in OMVPE-grown material also applies to HVPE. The tendencies for ordering to increase with deposition temperature and decrease with growth rate show that kinetics limit the degree of ordering, also in agreement with OMVPE trends.

74 ATOMIC AND MOLECULAR PHYSICS↗

Impurity-Induced Anomalous Thermal Hall Effect in Chiral Superconductors

Chiral superconductors exhibit novel transport properties that rely on the topology of the order parameter, topology of the Fermi surface, the spectrum of bulk and edge Fermionic excitations, and the structure of the impurity potential. In the case of electronic heat transport, impurities induce an anomalous (zero-field) thermal Hall conductivity that is easily orders of magnitude larger than the quantized edge contribution. The effect originates from branch-conversion scattering of Bogoliubov quasiparticles by the chiral order parameter, induced by potential scattering. The former transfers angular momentum between the condensate and the excitations that transport heat. The anomalous thermal Hall conductivity is shown to depend to the structure of the electron-impurity potential, as well as the winding number ν of the chiral order parameter $Δ(p)=|Δ(p)|e^{iν\phi_p}$. The results provide quantitative formulas for interpreting heat transport experiments seeking to identify broken T and P symmetries, as well as the topology of the order parameter for chiral superconductors.

36 MATERIALS SCIENCE↗

Probing configurational disorder in ZnGeN 2 using cluster-based Monte Carlo

ZnGeN 2 is sought as a semiconductor with comparable lattice constant to GaN and tunable band gap for integration in optoelectronic devices. Configurational disorder on the cation sublattice of ZnGeN 2 can strongly modify the electronic structure compared to the ordered material, and both ordered and disordered forms of ZnGeN 2 are candidates for light-emitting diodes and other emitters. The nonisovalent character of the disordered species (Zn 2+ and Ge 4+ ) subjects the cation ordering to strong short-range order effects. To model these effects, we use Monte Carlo (MC) simulations utilizing a cluster expansion to approximate formation enthalpy. Representative disordered configurations in 1024-atom supercells are relaxed in density functional theory calculations. From the MC structures, we extract a short-range order parameter (the N-cation coordination motif), and two long-range order parameters (Bragg-Williams and stretching parameters), and examine their correlations. We perform a thermodynamic integration to determine the mixing entropy and free energy. ZnGeN 2 exhibits a first-order phase transition with pronounced discontinuities in enthalpy and entropy, as well as in the structural order parameters. We discuss the relationship between the effective temperature used in the MC simulation and the growth temperatures in experiment in relation to the crossover from the nonequilibrium to the equilibrium growth regime. This work expands on current models of site disorder in ZnGeN 2 and provides atomic structure models with a systematic variation of the degree of cation disorder.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Charge density wave fluctuation driven composite order in layered kagome metals

The newly discovered kagome metals AV 3 Sb 5 (A = K, Rb, Cs) offer an exciting route to study exotic phases arising due to interplay between electronic correlations and topology. Besides superconductivity, these materials exhibit a charge-density wave (CDW) phase occurring at ~ 100 K, whose origin still remains elusive. The robust multi-component 2 × 2 CDW phase in these systems is of great interest due to the presence of an unusually large anomalous Hall effect. In quasi-2D systems with weak inter-layer coupling fluctuation driven exotic phases may appear. In particular in systems with multi-component order parameters fluctuations may lead to establishment of composite order when only products of individual order parameters condense while the individual ones themselves remain disordered. Here we argue that such fluctuation-driven regime of composite CDW order may exist in thin films of kagome metals above the CDW transition temperature. It is suggested that the melting of the Trihexagonal state in the material doped way from the van Hove singularities gives rise to a pseudogap regime where the spectral weight is concentrated in small pockets and most of the original Fermi surface is gapped. Our findings suggest possible presence of exotic phases in the weakly coupled layered kagome metals, more so in the newly synthesized thin films of kagome metals.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Molecular Static Third-Order Polarizabilities of Carbon-Cage Fullerenes and their Correlation with Three Geometric Parameters: Group Order, Aromaticity, and Size

Calculations were performed on the valence contribution to the static molecular third-order polarizabilities (gamma) of thirty carbon-cage fullerenes (C60, C70, five isomers of C78, and twenty-three isomers of C84). The molecular structures were obtained from B3LYP/STO-3G calculations. The values of the tensor elements and an associated numerical uncertainty were obtained using the finite-field approach and polynomial expansions of orders four to eighteen of polarization versus static electric field data. The latter information was obtained from semiempirical calculations using the AM1 hamiltonian.

Moore, Craig E.↗

Phase diagrams of single-layer two-dimensional transition metal dichalcogenides: Landau theory

We report single layer (SL) two-dimensional transition metal dichalcogenides (TMDs), such as MoS 2 , ReS 2 , WSe 2 , and MoTe 2 have now become the focus of intensive fundamental and applied research due to their intriguing and tunable physical properties. These materials exhibit a broad range of structural phases that can be induced via elastic strain, chemical doping, and electrostatic field effect. These transformations, in turn, can open and close the band gap of SL-TMDs, leading to metal-insulator transitions, and lead to emergence of more complex quantum phenomena. These considerations necessitate detailed understanding of the mesoscopic mechanisms of these structural phase transitions. Here we develop a Landau-type thermodynamic description of SL-TMDs on a sample of the SL-(MoS 2 ) 1-x -(ReS 2 ) x system and analyze the free energy surfaces, phase diagrams, and order parameter behavior. Our results predict the existence of multiple structural phases with two-, six- and 12-fold degenerated energy minima for in-plane and out-of-plane order parameters. This analysis suggests that out-of-plane ferroelectricity can exist in many of these phases, with the switchable polarization being proportional to the out-of-plane order parameter. We further predict that the domain walls in SL-(MoS 2 ) 1-x -(ReS 2 ) x should become conductive above a certain strain threshold.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Command of active and responsive elastomers by topological defects and patterns

The project resulted in the development of stimuli-responsive elastomer coatings formed by photopolymerized liquid crystal molecules. The molecular orientation of the liquid crystal elastomers is coupled to rubber-like elasticity. The project developed an approach to produce complex patterns of molecular orientation in elastomers by employing photopatterning technique based on plasmonic metamasks that convert unpolarized incoming light into a transmitted light beam with spatially-varying linear polarization. The polarization-modulated light beam aligns the film of azodye molecules, which serves as the substrate for liquid crystal elastomer coating. The alignment of the azodye molecules follows the polarization pattern with a phase shift of the alignment direction by 90 degrees. The azodye substrate imposes orientation of molecules in the adjacent liquid crystal layer which is preserved after the monomers are polymerized into a liquid crystal elastomer. By using differently aligned substrates one could produce various patterns of molecular orientation in the elastomer, which dictate the mechanical properties of the elastomers and their response to external stimuli such as temperature, humidity, and ultraviolet irradiation. In particular, the predesigned molecular orientation produces deterministically defined topography response of the liquid crystal coatings, in which their free surface changes from flat to either locally elevated or locally compressed, depending on whether the pre-inscribed molecular orientation is predesigned with a splay or bend deformation. The research established that the mechanism of this deterministic relation between the dynamic topography and pre-inscribed molecular orientation is rooted in the change of the tensorial order parameter that characterizes the degree of orientational order, and in the contraction/expansion of polymer molecules in response to the order parameter changes. In particular, lowering the order causes contraction of polymer globules. The work demonstrated that the dynamic topography of elastomer coatings with three-dimensional pattern of orientational order includes up and down motion of the material (along the normal to the coating) and also shifts in lateral directions. We demonstrate that the deterministic relationship between the orientational patterns and variations of coatings’ profile is caused by forces that are mathematically equivalent to active forces in the system of “pullers” or “pushers” in active matter, demonstrating universality of the out-of-equilibrium behavior. The project resulted in theoretical models that predictively describe the dynamic response of coatings. The topography of liquid crystal elastomers is sensitive to stimuli such as temperature, ultraviolet irradiation, and humidity. The applicability of the patterned approach to produce dynamic coatings was expanded from nematic elastomers to their smectic counterparts. The project combined a battery of imaging techniques to unveil the dynamic properties of the coatings, ranging from the state-of-the-art dynamic holographic microscopy to fluorescent confocal polarizing microscopy and polscope microscopy. The research uncovered potential applications of dynamic elastomer coatings, such as control placement of colloidal particles, engineering of living tissues with alignment patterns of cells that follow the orientational order of elastomer coatings. The approach developed in the project can be used to design programmable dynamic coatings with functionalities that mimic biological tissues such as skin.

36 MATERIALS SCIENCE↗

Simulation and characterization of cation disorder in ZnGeP 2

New optoelectronic materials are needed for improving the efficiency and reliability of devices such as solar cells. Cation ordering presents one means of controlling optoelectronic properties while introducing potential to also diversify the mineral constituents of electronic devices; however, the mechanisms of ordering are not yet well understood. To better understand cation ordering in a system integratable with current devices, we assess short- and long-range order parameters of ZnGeP 2 , a material closely lattice matched to Si. Structures are simulated using cluster-based Monte Carlo and first-principles calculations to compare structural distortions, periodicity, and local coordination environments in ZnGeP 2 to experimental data both from the literature and presented here. Comparing order parameters, we relate the transition in order parameters of ZnGeP 2 to that of ZnGeP 2 , discuss the reduction of band gaps with disorder, and show that traditional structural characterization alone is insufficient for understanding order in ZnGeP 2 .

14 SOLAR ENERGY↗

Interplay between magnetism and superconductivity in UTe 2

Time-reversal symmetry breaking (TRSB) in UTe 2 was inferred from observations of a spontaneous Kerr response in the superconducting state after cooling in zero magnetic field, while a finite c-axis magnetic field training was further used to determine the nature of the nonunitary composite order parameter of this material. Here, we present an extensive study of the magnetic-field-trained Kerr effect, which unveils a unique critical state of pinned “ferromagnetic vortices.” We show that a remanent Kerr signal appears following the removal of a training magnetic field, which reflects the response to the TRSB order parameter and the external magnetic field through the paramagnetic susceptibility. Furthermore, this unambiguously demonstrates the importance of the magnetic fluctuations and their intimate relation to the composite order parameter. Focusing the beam on the center of the sample, we are able to accurately determine the maximum field that is screened by the critical state and the respective critical current. Measurements in the presence of magnetic field show the tendency of the superconductor to produce shielding currents that oppose the increase in vortex-induced magnetization due to the diverging paramagnetic susceptibility.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

A consistent and conservative volume distribution algorithm and its applications to multiphase flows using Phase-Field models

In the present study, the multiphase volume distribution problem, where there can be an arbitrary number of phases, is addressed using a consistent and conservative volume distribution algorithm. The proposed algorithm satisfies the summation constraint, the conservation constraint, and the consistency of reduction. The first application of the volume distribution algorithm is to determine the Lagrange multipliers in multiphase Phase-Field models that enforce the mass conservation, and a multiphase conservative Allen-Cahn model that satisfies the consistency of reduction is developed. A corresponding consistent and conservative numerical scheme is developed for the model. The multiphase conservative Allen-Cahn model has a better ability than the multiphase Cahn-Hilliard model to preserve under-resolved structures. The second application is to develop a numerical procedure, called the boundedness mapping, to map the order parameters, obtained numerically from a multiphase model, into their physical interval, and at the same time to preserve the physical properties of the order parameters. Along with the consistent and conservative schemes for the multiphase Phase-Field models, the numerical solutions of the order parameters are reduction consistent, conservative, and bounded, which are theoretically analyzed and numerically validated. Then, the multiphase Phase-Field models are coupled with the momentum equation by satisfying the consistency of mass conservation and the consistency of mass and momentum transport, thanks to the consistent formulation. Finally, it is demonstrated that the proposed model and scheme converge to the sharp-interface solution and are capable of capturing the complicated multiphase dynamics even when there is a large density and/or viscosity ratio.

42 ENGINEERING↗

Revealing the unusual rate-dependent mechanical behaviors of nematic liquid crystal elastomers

Liquid crystal elastomers (LCEs) exhibit unique mechanical properties of soft elasticity and enhanced energy dissipation with rate dependency. They are potentially transformative materials for applications in mechanical impact mitigation and vibration isolation. However, previous studies have primarily focused on the mechanics of LCEs under equilibrium and quasistatic loading conditions. Critical knowledge gaps exist in understanding their rate-dependent behaviors, which are a complex mixture of traditional network viscoelasticity and the soft elastic behaviors with changes in the mesogen orientation and order parameter. Together, these inelastic mechanisms lead to unusual rate-dependent energy absorption responses of LCEs. In this work, we developed a viscoelastic constitutive theory for monodomain nematic LCEs to investigate how multiple underlying sources of inelasticity manifest in the rate-dependent and dissipative behaviors of monodomain LCEs. Here, the theoretical modeling framework combines the neo-classical network theory with evolution rules for the mesogen orientation and order parameter with conventional viscoelasticity. The model is calibrated with uniaxial tension and compression data spanning six decades of strain rates. The established 3D constitutive model enables general loading predictions taking the initial mesogen orientation and order parameter as inputs. Additionally, parametric studies were performed to further understand the rate dependence of monodomain LCEs in relation to their energy absorption characteristics. Based on the parametric studies, particularly loading scenarios are identified as conditions where LCEs outperform conventional elastomers regarding energy absorption.

42 ENGINEERING↗

Electron– K -phonon interaction in twisted bilayer graphene

We develop an analytic theory to describe the interaction between electrons and K phonons and study its influence on superconductivity in the bare bands of twisted bilayer graphene (TBG). We find that, due to symmetry and the two-center approximation, only one optical K phonon (160meV) of graphene is responsible for the intervalley electron-phonon interaction. Here, this phonon has recently been found in angular-resolved photoemission spectroscopy to be responsible for replicas of the TBG flat bands. By projecting the interaction to the TBG flat bands, we perform the full symmetry analysis of the phonon-mediated attractive interaction and pairing channels in the Chern basis, and show that several channels are guaranteed to have gapless order parameters. From the linearized gap equations, we find that the highest T c pairing induced by this phonon is a singlet gapped s-wave inter-Chern-band order parameter, followed closely by a gapless nematic d-wave intra-Chern-band order parameter. We justify these results analytically, using the topological heavy-fermion mapping of TBG which has allowed us to obtain an analytic form of a phonon-mediated attractive interaction and to analytically solve the linearized and T = 0 gap equations. For the intra-Chern-band channel, the nematic state with nodes is shown to be stabilized in the chiral flat-band limit. While the flat-band Coulomb interaction can be screened sufficiently enough—around the Van Hove singularities—to allow for electron-phonon based superconductivity, it is unlikely that this effect can be maintained in the lower density of states excitation bands around the correlated insulator states.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Superconducting gap symmetry from Bogoliubov quasiparticle interference analysis on Sr 2 RuO 4

Here, the nature of the superconducting order parameter in Sr 2 RuO 4 has generated intense interest in recent years. Since the superconducting gap is very small, high-resolution methods such as scanning tunneling spectroscopy might be the best chance to directly resolve the gap symmetry. Recently, a Bogoliubov quasiparticle interference imaging (BQPI) experiment has suggested that the d x 2 –y 2 gap symmetry is appropriate for Sr 2 RuO 4 . In this work, we use a material-specific theoretical approach based onWannier functions of the surface of Sr 2 RuO 4 to calculate the continuum density of states as detected in scanning tunneling microscopy experiments. We examine several different proposed gap order parameters and calculate the expected BQPI pattern for each case. Comparing to the available experimental data, our results suggest that a s' + id xy gap order parameter is the most probable state, but the measured BQPI patterns still display features unaccounted for by the theory for any of the states currently under discussion.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Nonequilibrium universality of the nonreciprocally coupled 𝑂⁡(𝑛 1 ) × 𝑂⁡(𝑛 2 ) model

Nonequilibrium dynamics play an important role in all contexts of physics, both classical and quantum as well as living and nonliving, so it is crucial to develop a foundational understanding of nonequilibrium phase transitions. In this work we investigate an important class of nonequilibrium dynamics in the form of nonreciprocal interactions. In particular we study how nonreciprocal coupling between two 𝑂⁡(𝑛𝑖) order parameters (with 𝑖 = 1,2) affects the universality at a multicritical point, extending the analysis of J. T. Young et al. [Phys. Rev. X 10, 011039 (2020)], which considered the case 𝑛 1 = 𝑛 2 = 1, i.e., a ℤ 2 × ℤ 2 model. We show that nonequilibrium fixed points (NEFPs) emerge for a broad range of 𝑛 1 ,𝑛 2 and exhibit intrinsically nonequilibrium critical phenomena, namely a violation of fluctuation-dissipation relations at all scales and underdamped oscillations near criticality in contrast to the overdamped relaxational dynamics of the corresponding equilibrium models. Furthermore, the NEFPs exhibit an emergent discrete scale invariance in certain physically relevant regimes of 𝑛 1 ,𝑛 2 , but not others, depending on whether the critical exponent 𝜈 is real or complex. The boundary between these two regions is described by an exceptional point in the renormalization group (RG) flow, leading to distinctive features in correlation functions and the phase diagram. Another contrast with the previous work is the number and stability of the NEFPs as well as the underlying topology of the RG flow. Lastly, we investigate an extreme form of nonreciprocity where one order parameter is independent of the other order parameter but not vice versa. Unlike the ℤ 2 × ℤ 2 model, which becomes nonperturbative in this case, we identify a distinct nonequilibrium universality class whose dependent field similarly violates fluctuation-dissipation relations but does not exhibit discrete scale invariance or underdamped oscillations near criticality.

Critical phenomena↗

Predictive outlook for experiments resolving prompt vs local redeposition of high- Z materials in tokamaks

High-Z plasma facing components redeposit within the sheath through a combination of two distinct mechanisms: prompt (or geometric-driven) and local (or sheath-driven) redeposition. Experimental efforts are needed to determine the leading-order parameters influencing prompt-vs-local trade-off, which sets the fraction of material entering the scrape-off layer. In preparation for such experiments, leading-order parameters are isolated within the PYEAD-RustBCA-GITR coupled net erosion code using Sobol’ sensitivity analysis. Then, experiments resolving prompt-vs-local trade-off under variation of these leading-order parameters are proposed using an isotopic coupon design with multifaceted diagnostic coverage. The measurability of these experiments is evaluated using synthetic diagnostics.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗