Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “symmetry projection”

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 199 records · Page 11

Coulomb Crystallization of Charged Microspheres Levitated in a Gas Discharge Plasma

The technical topic of the project was the experimental observation of Coulomb crystallization of charged microspheres levitated in a gas discharge plasma. This suspension, sometimes termed a dusty plasma, is closely analogous to a colloidal suspension, except that it has a much faster time response, is more optically thin, and has no buoyancy forces to suspend the particles. The particles are levitated by electric fields. Through their collective Coulomb repulsions, the particles arrange themselves in a lattice with a crystalline symmetry, which undergoes an order-disorder phase transition analogous to melting when the effective temperature of the system is increased. Due to gravitational sedimentation, the particles form a thin layer in the laboratory, so that the experimental system is nearly 2D, whereas in future microgravity experiments they are expected to fill a larger volume and behave like a 3D solid or liquid. The particles are imaged using a video camera by illuminating them with a sheet of laser light. Because the suspension is optically thin, this imaging method will work as well in a 3D microgravity experiment as it does in a 2D laboratory system.

Goree, John↗

Theoretical Investigation of Thermo-Mechanical Behavior of Carbon Nanotube-Based Composites Using the Integral Transform Method

This project uses the integral transform technique to model the problem of nanotube behavior as an axially symmetric system of shells. Assuming that the nanotube behavior can be described by the equations of elasticity, we seek a stress function x which satisfies the biharmonic equation: del(exp 4) chi = [partial deriv(r(exp 2)) + partial deriv(r) + partial deriv(z(exp 2))] chi = 0. The method of integral transformations is used to transform the differential equation. The symmetry with respect to the z-axis indicates that we only need to consider the sine transform of the stress function: X(bar)(r,zeta) = integral(from 0 to infinity) chi(r,z)sin(zeta,z) dz.

Pawloski, Janice S.↗

One Dimensional Topological Nanomaterials and Superconductivity (Final Technical Report)

The project synthesized and made transport properties of topological crystalline insulator tin telluride (SnTe) nanowires and other topological materials, such as MoTe 2 and WTe 2 . The transport properties of the nanowires were correlated with the diameter of the nanowire, in terms of the ferroelectric phase transition temperature, magnetoresistance response, and Josephson junction behaviors as a function of the nanowire diameter. The most important observation was the proximity-induced superconductivity in SnTe nanowires, which suggests time reversal symmetry breaking.

36 MATERIALS SCIENCE↗

A search for Heavy Neutral Leptons with masses up to 150 MeV and studies on scintillation light detection in SBND

Introduction Neutrino physics is one of today’s most promising fields in particle physics. The discovery of neutrino oscillations is an evidence of their mass not being null. However, the Standard Model of elementary particles (SM) does not contain a mechanism that generates them, so we have to search beyond the SM to find answers. Other unknowns remain open, such as the amount of CP symmetry violated in the lepton sector that could help explain the differences between matter and antimatter that formed the present universe. Nor is it understood why these particles are so light, more than six orders of magnitude below the next fermion, the electron. A possible explanation is the existence of heavier neutrino mass states usually referred as heavy neutral leptons (HNLs). One of the most advanced technologies for detecting neutrino interactions is the liquid argon time projection chamber (LArTPC). LArTPC detectors form a 3D image of the ionized electrons in the interaction, producing a detailed description of each neutrino event. Argon is also a prolific scintillator (40k photons/MeV) with light signals indicating the interaction time. The Short-Baseline Near Detector (SBND), a 112-active mass LArTPC at Fermilab, will measure millions of neutrino interactions from the Booster Neutrino Beam (BNB). As the near (L=110 m) detector of the Short-Baseline neutrino program, SBND will search for a fourth light neutrino (∼ eV) in the coming years. This hypothesis could explain the anomalous data observed in the LSND and MiniBooNE experiments.

Álvarez Garrote, Rodrigo [Madrid U.] (ORCID:000000↗

Dynamics and control of coherent structures in the turbulent wall layer: An overview

We expand the velocity field in the vicinity of the wall in empirical eigenfunctions obtained from experiment. Truncating our system and using Galerkin projection, we obtain a closed set of non-linear ordinary differential equations with ten degrees of freedom. We find a rich dynamical behavior, including in particular a heteroclinic attracting orbit giving rise to intermittency. The intermittent jump from one attracting point to the other resembles in many respects the bursts observed in experiments. Specifically, the time between jumps and the duration of the jumps, is approximately that observed in a burst; the jump begins with the formation of a narrowed and intensified updraft, like the ejection phase of a burst, and is followed by a gentle, diffuse downdraft like the sweep phase of a burst. The magnitude of the Reynolds stress spike produced during a burst is limited by our truncation. The behavior is quite robust, much of it being due to the symmetries present (Aubry's group has examined dimensions up to 128 with persistence of the global behavior). We have examined eigenvalues and coefficients obtained from experiment, and from exact simulation, which differ in magnitude. Similar behavior is obtained in both cases; in the latter case, the heteroclinic orbits connect limit cycles instead of fixed points, corresponding to cross-stream waving of the streamwise rolls. The bifurcation diagram remains structurally similar, but somewhat distorted. The role of the pressure term is made clear - it triggers the intermittent jumps, which otherwise would occur at longer and longer intervals, as the system trajectory is attracted closer and closer to the heteroclinic cycle. The pressure term results in the jumps occurring at essentially random times, and the magnitude of the signal determines the average timing. Stretching of the wall region shows that the model is consistent with observations of polymer drag reduction. Change of the third order coefficients, corresponding to acceleration or deceleration of the mean flow, changes the heteroclinic cycles from attracting to repelling, increasing or decreasing the stability, in agreement with observations. The existence of fixed points is an artifact introduced by the projection; however, a decoupled model still displays the rich dynamics. Numerous assumptions made in Aubry et al. (1988) can now be proved exactly. Feeding back eigenfuncitons with the proper phase can delay the bursting, (the heteroclinic jump to the other fixed point), decreasing the drag. It is also possible to speed up the bursting, increasing mixing to control separation. Our approach is optimal for short time tracking in control.

Berkooz, Gal↗

Perturbative Quantum Monte Carlo Method for Nuclear Physics

While first order perturbation theory is routinely used in quantum Monte Carlo (QMC) calculations, higher-order terms present significant numerical challenges. We present a new approach for computing perturbative corrections in projection QMC calculations. Here we demonstrate the method by computing nuclear ground state energies up to second order for a realistic chiral interaction. We calculate the binding energies of several light nuclei up to 16 O by expanding the Hamiltonian around the Wigner SU(4) limit and find good agreement with data. In contrast to the natural ordering of the perturbative series, we find remarkably large second-order energy corrections. This occurs because the perturbing interactions break the symmetries of the unperturbed Hamiltonian. Our method is free from the sign problem and can be applied to QMC calculations for many-body systems in nuclear physics, condensed matter physics, ultracold atoms, and quantum chemistry.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Virtual and Experimental Visualization of Flows in Packed Beds of Spheres Simulating Porous Media Flows

A videotape presentation of flow in a packed bed of spheres is provided. The flow experiment consisted of three principal elements: (1) an oil tunnel 76.2 mm by 76.2 mm in cross section, (2) a packed bed of spheres in regular and irregular arrays, and (3) a flow characterization methodology, either (a) full flow field tracking (FFFT) or (b) computational fluid dynamic (CFD) simulation. The refraction indices of the oil and the test array of spheres were closely matched, and the flow was seeded with aluminum oxide particles. Planar laser light provided a two-dimensional projection of the flow field, and a traverse simulated a three-dimensional image of the entire flow field. Light focusing and reflection rendered the spheres black, permitting visualization of the planar circular interfaces in both the axial and transverse directions. Flows were observed near the wall-sphere interface and within the set of spheres. The CFD model required that a representative section of a packed bed be formed and gridded, enclosing and cutting six spheres so that symmetry conditions could be imposed at all cross-boundaries. Simulations had to be made with the flow direction at right angles to that used in the experiments, however, to take advantage of flow symmetry. Careful attention to detail was required for proper gridding. The flow field was three-dimensional and complex to describe, yet the most prominent finding was flow threads, as computed in the representative 'cube' of spheres with face symmetry and conclusively demonstrated experimentally herein. Random packing and bed voids tended to disrupt the laminar flow, creating vortices.

Hendricks, R. C.↗

Toward Realistic Dynamics of Rotating Orbital Debris, and Implications for Lightcurve Interpretation

Optical observations of rotating space debris near GEO contain important information on size, shape, composition, and rotational states, but these aspects are difficult to extract due to data limitations and the high number of degrees of freedom in the modeling process. For tri-axial rigid debris objects created by satellite fragmentations, the most likely initial rotation state has a large component of initial angular velocity directed along the intermediate axis of inertia, leading to large angular reorientations of the body on the timescale of the rotation period. This lends some support to the simplest possible interpretation of light curves -- that they represent sets of random orientations of the objects of study. However, effects of internal friction and solar radiation are likely to cause significant modification of rotation states within a time as short as a few orbital periods. In order to examine the rotational dynamics of debris objects under the influences of these effects, a set of seven first-order coupled equations of motion were assembled in state form: three are Euler equations describing the rates of change of the components of angular velocity in the body frame, and four describe the rates of change of the components of the unit quaternion. Quaternions are a four-dimensional extension of complex numbers that form a seamless, singularity-free representation of body orientation on S3. The Euler equations contain explicit terms describing torque from solar radiation in terms of spherical harmonics, and terms representing effects of a prescribed rate of internal friction. Numerical integrations of these equations of motion are being performed, and results will be presented. Initial tests show that internal friction without solar radiation torque leads to rotation about the maximum principal axis of inertia, as required, and solar radiation torque is expected to lead to spin-up of objects. Because the axis of maximum rotational inertia tends to be roughly coincident with the normal to the largest projected cross-sectional area, internal friction is expected to lead to reduced variation of light curve amplitudes at a given phase angle, but a large dependence of the same on phase angle. At a given phase angle, databases are generated which contain reflected intensities for comprehensive sets of equally-likely orientations, represented as unit quaternions. When projected onto three dimensions (S2) and color-coded by intensity, the set is depicted as points within a solid, semi-transparent unit sphere, within which all possible reflected intensities for an object at a given phase angle may be inspected simultaneously. Rotational sequences are represented by trajectories through the sphere. Databases are generated for each of a set of phase angles separately, forming a comprehensive dataset of reflected intensities spanning all object orientations and solar phase angles. Symmetries in the problem suggest that preferred rotation states are likely, defined relative to the object-sun direction in inertial space and relative to the maximum principal axis of inertia in the body coordinate system. Such rotation states may greatly simplify the problem of light curve interpretation by reducing the number of degrees of freedom in the problem.

Ojakangas, Gregory W.↗

KOTO vs. NA62 dark scalar searches

The two kaon factories, KOTO and NA62, are at the cutting edge of the intensity frontier, with an unprecedented numbers of long lived and charged Kaons,~10 13 , being measured and analyzed. These experiments have currently a unique opportunity to search for dark sectors. In this paper, we demonstrate that searches done at KOTO and NA62 are complementary, both probing uncharted territories. We consider two qualitatively different physics cases. In the first, we analyze models of axion-like-particles (ALP) which couple to gluons or electroweak gauge bosons. In the second, we introduce a model based on an approximate strange flavor symmetry that leads to a strong violation of the Grossman-Nir bound. For the first scenario, we design a new search strategy for the KOTO experiment,K L → π 0 a →4γ. Its expected sensitivity on the branching ratio is at the level of 10 -9 . This demonstrates the great potential of KOTO as a discovery machine. In addition, we revisit other bounds on ALPs from Kaon factories, highlighting the main sources of theoretical uncertainty, and collider experiments, and show new projections. For the second scenario, we show that the model may be compatible with the preliminary analysis of the KOTO-data that shows a hint for New Physics.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Interacting fermions in narrow-gap semiconductors with band inversion

Abstract Highly unconventional behavior of the thermodynamic response functions has been experimentally observed in a narrow gap semiconductor samarium hexaboride. Motivated by these observations, we use renormalization group technique to investigate many-body instabilities in the f -orbital narrow gap semiconductors with band inversion in the limit of weak coupling. By projecting out the double occupancy of the f -states we formulate a low-energy theory describing the interacting particles in two hybridized electron- and hole-like bands. The interactions are assumed to be weak and short-ranged. We take into account the difference between the effective masses of the quasiparticles in each band. Upon carrying out the renormalization group analysis we find that there is only one stable fixed point corresponding to the excitonic instability with time-reversal symmetry breaking for small enough mismatch between the effective masses.

Physics↗

Gauge invariants of linearized gravity with a general background metric

In linearized gravity with distributed matter, the background metric has no generic symmetries, and decomposition of the metric perturbation into global normal modes is generally impractical. This complicates the identification of the gauge-invariant part of the perturbation, which is a concern, for example, in the theory of dispersive gravitational waves (GWs) whose energy–momentum must be gauge-invariant. Here, we propose how to identify the gauge-invariant part of the metric perturbation and the six independent gauge invariants per se for an arbitrary background metric. For the Minkowski background, the operator that projects the metric perturbation on the invariant subspace is proportional to the well-known dispersion operator of linear GWs in vacuum. For a general background, this operator is expressed in terms of the Green’s operator of the vacuum wave equation. Further, if the background is smooth, it can be found asymptotically using the inverse scale of the background metric as a small parameter.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Degenerate coupled-cluster theory

A size-extensive, converging, black-box, ab initio coupled-cluster (ΔCC) ansatz is introduced that computes the energies and wave functions of states from any degenerate or nondegenerate Slater-determinant references with any numbers of α- and β-spin electrons, any patterns of orbital occupancy, any spin multiplicities, and any spatial symmetries. For a nondegenerate reference, it reduces to the single-reference coupled-cluster ansatz. For a degenerate multireference, it is a natural coupled-cluster extension of degenerate Møller–Plesset perturbation (ΔMP) theory. For ionized and electron-attached references, it is a coupled-cluster Green’s function, although the present theory is convergent toward the full-configuration-interaction limits, while the Feynman–Dyson many-body Green’s function (MBGF) theory generally is not. Its single-excitation instance is a projection Hartree–Fock theory as per the Thouless theorem, which may be useful for core ionizations, high-spin states, and possibly electron affinities. Additionally, a new multireference coupled-cluster theory for a general model space is developed. This quasidegenerate coupled-cluster (QCC) theory is exactly converging, but not black-box, and intended for strong correlation. Determinant-based, general-order algorithms of ΔCC and QCC theories are implemented and compared with configuration-interaction (CI) and equation-of-motion coupled-cluster (EOM-CC) theories through octuple excitations and with ΔMP and MBGF theories up to the nineteenth order. An algebraic, optimal-scaling algorithm of the ΔCC theory is computer-synthesized at the levels of single excitations (ΔCCS) and of single and double excitations (ΔCCSD). As a result, the order of performance is QCC ≈ ΔCC > EOM-CC > CI at the same order or QCC ≈ ΔCC > ΔMP > MBGF at the same cost scaling.

Hirata, So [University of Illinois at Urbana-Champ↗

New Theories for Neutrino Masses and Dark Matter

The main goal of the project is to investigate new theories for neutrino masses and dark matter, and understand their testability in different experiments. The PI proposed several mechanisms for neutrino masses in the context of B-L theories, left-right symmetric theories, supersymmetric theories and grand unified theories. The PI would like to investigate new ways to understand the origin of neutrino masses where the seesaw scale is in the multi-TeV region. One can have a multi-TeV seesaw scale in scenarios where the same symmetry relevant for neutrino masses defines the dark matter relic density, or in the minimal supersymmetric theory based on local B-L. In these theories one can have new ways to understand the testability of the origin of neutrino masses at colliders and low energy experiments. In these theories one predicts the possibility to observe lepton number violating signatures at the LHC, one can predict large contributions to lepton number violating processes such as mu to e conversion, neutrinoless double beta decay experiments and others. The nature of the dark matter in the Universe is one of the most important problems in cosmology. Two of the most popular candidates are the Axions and the Weakly Interacting Massive Particles (WIMPs). The PI would like to investigate the implications of a new electroweak theta term similar to the QCD vacuum angle but in the SU(2) gauge sector of the SM which is physical if the baryon and lepton numbers are broken symmetries. We will investigate the implications of having a new electroweak theta term dark matter portal, the implications for baryogenesis and understand the constraints coming from different experiments such as the searches for electric dipole moments and axion experiments. The implications of having the electroweak theta term in theories with different sources of baryon and lepton number violation will be investigated. WIMPs are perhaps the most appealing candidates for many reasons. The PI would like to investigate the properties of WIMPs in different gauge theories, study the implications of the cosmological bounds on the WIMPs relic density for the symmetry breaking scale in models for new physics. The PI would like to investigate in great detail the predictions for gamma lines, neutrino lines and other striking signatures which can be present in these theories. These studies can have a profound impact in the testability of different theories for physics beyond the Standard Model.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Quantum materials for nanosensing and fault-tolerant quantum computing

New concepts of symmetry related to topological order emerged from the discovery of the fractional quantum Hall effect and high-temperature superconductivity in strongly correlated electron systems. This led to the study of quantum materials-- materials exhibiting emergent quantum phenomena with no classical analogues. While these materials have engendered exciting basic materials science and physics, realizing novel devices is a key challenge in the field. The goal of this proposal is to harness the unique properties of topological materials for quantum computing and quantum sensing applications. In this project, we investigated a variety of topological superconducting platforms and identified three technologies that can benefit from their quantum properties: quantum memory, single-photon detection, and non reciprocal electronics. The platforms developed in this work will be broadly useful to National Security and Basic Science.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Twisted bilayer graphene. I. Matrix elements, approximations, perturbation theory, and a k · p two-band model

We investigate the twisted bilayer graphene (TBG) model of Bistritzer and MacDonald (BM) [Bistritzer and MacDonald, Proc. Natl. Acad. Sci. 108, 12233 (2011)] to obtain an analytic understanding of its energetics and wave functions needed for many-body calculations. We provide an approximation scheme for the wave functions of the BM model, which first elucidates why the BM K M -point centered original calculation containing only four plane waves provides a good analytical value for the first magic angle (θ M ≈ 1°). The approximation scheme also elucidates why most of the many-body matrix elements in the Coulomb Hamiltonian projected to the active bands can be neglected. By applying our approximation scheme at the first magic angle to a Γ M -point centered model of six plane waves, we analytically understand the reason for the small Γ M -point gap between the active and passive bands in the isotropic limit w 0 = w 1 . Furthermore, we analytically calculate the group velocities of the passive bands in the isotropic limit, and show that they are almost doubly degenerate, even away from the Γ M point, where no symmetry forces them to be. Furthermore, moving away from the Γ M and K M points, we provide an explicit analytical perturbative understanding as to why the TBG bands are flat at the first magic angle, despite the first magic angle is defined by only requiring a vanishing K M -point Dirac velocity. We derive analytically a connected “magic manifold” w 1 = $2\sqrt{1 + w^{2}_{0}}$ $-\sqrt{2 + 3w^2_0}$, on which the bands remain extremely flat as w 0 is tuned between the isotropic (w 0 = w 1 ) and chiral (w 0 = 0) limits. We analytically show why going away from the isotropic limit by making w 0 less (but not larger) than w 1 increases the Γ M -point gap between the active and the passive bands. Finally, by perturbation theory, we provide an analytic Γ M point k ∙ p two-band model that reproduces the TBG band structure and eigenstates within a certain w 0 , w 1 parameter range. Further refinement of this model are discussed, which suggest a possible faithful representation of the TBG bands by a two-band Γ M point k ∙ p model in the full w 0 , w 1 parameter range.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Finite and infinite matrix product states for Gutzwiller projected mean-field wave functions

Matrix product states (MPS) and “dressed” ground states of quadratic mean fields (e.g., Gutzwiller projected Slater determinants) are both important classes of variational wave functions. This latter class has played important roles in understanding superconductivity and quantum spin liquids. We present a method to obtain both the finite and infinite MPS (iMPS) representation of the ground state of an arbitrary fermionic quadratic mean-field Hamiltonian (which in the simplest case is a Slater determinant and in the most general case is a Pfaffian). We also show how to represent products of such states (e.g., determinants times Pfaffians). From this representation one can project to single occupancy and evaluate the entanglement spectra after Gutzwiller projection. We then obtain the MPS and iMPS representation of Gutzwiller projected mean-field states that arise from the variational slave-fermion approach to the S = 1 bilinear-biquadratic quantum spin chain. To accomplish this, we develop an approach to orthogonalize degenerate iMPS to find all the states in the degenerate ground-state manifold. We find the energies of the MPS and iMPS states match the variational energies closely, indicating the method is accurate and there is minimal loss due to truncation error. We then present an exploration of the entanglement spectra of projected slave-fermion states, exploring their qualitative features and finding good qualitative agreement with the respective exact ground-state spectra found from density matrix renormalization group.

1-dimensional spin chains↗

Numerically Optimized Coronagraph Designs for the Habitable Exoplanet Imaging Mission (HabEx) Concept

The primary science goal of the Habitable Exoplanet Imaging Mission (HabEx), one of four candidate flagship missions under investigation, is to image and spectrally characterize Earth-like exoplanets. It is well known that pupil obscurations degrade coronagraphic performance and complicate coronagraph design, so HabEx is planned to have an off-axis, unobscured primary mirror. We utilize the circular symmetry of the aperture to investigate 1D-radial coronagraph optimization methods that are prohibitively time-consuming or intractable in 2D, such as diffractive pupil remapping and concurrent, multi-plane optimization. We also directly constrain sensitivities to dynamic, low-order Zernike aberrations, which are separable in polar coordinates and can thus be propagated as 1D-radial integrals. The mask technologies in our designs claim heritage from the extensive modeling and testbed experiments performed by the Wide-Field Infrared Survey Telescope (WFIRST) Coronagraph Instrument (CGI) project. In this paper, we detail our optimization methods and outline future work to complete our design survey.

Balasubramaian, Kunjithapatham↗

A hybrid 3D/2D field response calculation for liquid argon detectors with PCB based anode plane

Liquid Argon Time Projection Chamber (LArTPC) technology is commonly utilized in neutrino detector designs. It enables detailed reconstruction of neutrino events with high spatial precision and low energy threshold. Its field response (FR) model describes the time-dependent electric currents induced in the anode-plane electrodes when ionization electrons drift nearby. An accurate and precise FR is a crucial input to LArTPC detector simulations and charge reconstruction. Established LArTPC designs have been based on parallel wire planes. It allows accurate and computationally economic two-dimensional (2D) FR models utilizing the translational symmetry along the direction of the wires. Recently, novel LArTPC designs utilize electrodes formed on printed circuit board (PCB) in the shape of strips with through holes. The translational symmetry is no longer a good approximation near the electrodes and a new FR calculation that employs regions with three dimensions (3D) has been developed. Extending the 2D models to 3D would be computationally expensive. Fortuitously, the nature of strips with through holes allows for a computationally economic approach based on the finite-difference method (FDM). In this paper, we present a new software package pochoir that calculates LArTPC field response for these new strip-based anode designs. This package combines 3D calculations in the volume near the electrodes with 2D far-field solutions to achieve fast and precise field response computation. We apply the resulting FR to simulate and reconstruct samples of cosmic-ray muons and 39 Ar decays from a Vertical Drift (VD) detector prototype operated at CERN. We find the difference between real and simulated data within 5%. Current state-of-the-art LArTPC software requires a 2D FR which we provide by averaging over one dimension and estimate that variations lost in this average are smaller than 7%.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗