Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Density matrix methods”

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 811 records · Page 45

Microporous structure in silicone polymers with polylactic acid and polyethylene glycol microspheres: Synthesis and characterization

Porosity in polymer parts could improve mechanical properties while reducing weight and density. Individual pores must be small enough that bulk mechanical properties are retained but significant enough to reduce the density of the material. Filler materials can be incorporated into polymer matrices and removed to create voids. Here, this paper focuses on creating spherical pores in silicone matrices by inserting and removing core materials or fillers during post-cure. A method was developed to incorporate pore formers into an uncured polymer resin and remove the pore former materials from the matrix. Polyethylene glycol (PEG) and polylactic acid (PLA) microspheres were created and incorporated into silicone matrices as pore formers. Microspheres with controllable size distributions were formed via an emulsion formation. Microspheres are removed from the matrix by solvent extraction and calcination. Thermogravimetric Analysis (TGA), Dynamic Mechanical Analysis (DMA), and optical imaging characterize and determine the successful creation of micropores in the polymer matrices.

36 MATERIALS SCIENCE↗

Features that Further Performance Limits of Nuclear Fuel Fabrication: Opportunities for Additive Manufacturing of Nuclear Fuels

One goal of the Transformational Challenge Reactor (TCR) program is to fully assess the opportunity space for additive manufacturing (AM) to overcome the material limitations restricting development and deployment of nuclear power. This report presents an initial analysis of the opportunities for using AM to significantly impact the nuclear fuel fabrication paradigm. The conventional limitations of the three basic solid fuel types are examined herein, and. possibilities for fabrication technologies to overcome these limitations are outlined. Areas that can be most significantly improved by AM are those in which conventional nuclear fuel fabrication routes are known to limit performance. For monolithic systems, one high priority area is fabrication of heterogeneous duplex and compositionally graded structures— fuel designs that bring significant benefits to fuel utilization. A critical first step in meeting this challenge is to successfully demonstrate AM’s ability to achieve representative density and microstructures for fuel materials such as UO 2 . The area with the highest potential for pairing AM methodology with the challenges of nuclear fuel deployment is advanced matrices for particle and dispersion fuels that are free from conventional processing constraints. AM has the potential to achieve packing fractions and matrix quality that are superior to those available through conventional methods. Finally, the use of digital twin data collection, which is possible through in situ build diagnostics, may greatly reduce uncertainty regarding distribution of uranium in these systems. A subsequent FY19 TCR milestone will examine the suitability of existing commercial systems to address these challenges.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Thermomechanical Microstructural Predictions of Fracture Nucleation of Zircaloy-4 Alloys With δ and ε Hydride Distributions

Abstract A crystalline dislocation-density formulation that was incorporated with a nonlinear finite-element (FE) method was utilized to understand and to predict the thermomechanical behavior of an hexagonal closest packed (h.c.p.) zircaloy system with hydrides with either face-centered cubic (f.c.c.) or body-centered cubic (b.c.c.) hydrides. This formulation was then used with a recently developed fracture methodology that is adapted for finite inelastic strains and multiphase crystalline systems to understand how different microstructurally based fracture modes nucleate and propagate. The interrelated microstructural characteristics of the different crystalline hydride and matrix phases with the necessary orientation relationships (ORs) have been represented, such that a detailed physical understanding of fracture nucleation and propagation can be predicted for the simultaneous thermomechanical failure modes of hydride populations and the matrix. The effects of volume fraction, morphology, crystalline structure, and orientation and distribution of the hydrides on simultaneous and multiple fracture modes were investigated for radial, circumferential, and mixed distributions. Another key aspect was accounting for temperatures changes due to the effects of thermal conduction and dissipated plastic work and their collective effects on fracture. For hydrided aggregates subjected to high temperatures, thermal softening resulted in higher ductility due to increased dislocation-density activity, which led to higher shear strain accumulation and inhibited crack nucleation and growth. The predictions provide validated insights into why circumferential hydrides are more fracture-resistant than radial hydrides for different volume fractions and thermomechanical loading conditions.

Engineering↗

Fabrication of Al{sub 2}O{sub 3}-20 vol.% Al nanocomposite powders using high energy milling and their sinterability

Highlights: Black-Right-Pointing-Pointer Al{sub 2}O{sub 3}/Al nanocomposite powders were prepared via high energy ball milling. After 20 h milling, the size of Al{sub 2}O{sub 3}-20 vol.% Al nanocomposite particles was in the range of 23-29 nm. A uniform distribution of nanosized Al reinforcement throughout the Al{sub 2}O{sub 3} matrix, coating the particles was successfully obtained. Black-Right-Pointing-Pointer There was no any sign of phase changes during the milling. A competition between the cold welding mechanism and the fracturing mechanism were found during milling and finally the above two mechanisms reached an equilibrium. Black-Right-Pointing-Pointer The highest value of relative density was obtained for the sintered bodies at 1500 Degree-Sign C. Black-Right-Pointing-Pointer The harness of the sintered composite was decreased while the fracture toughness was improved after addition Al into alumina. -- Abstract: In this study, alumina-based matrix nanocomposite powders reinforced with Al particles were fabricated and investigated. The sinterability of the prepared nanocomposite powder at different firing temperature was also conducted. Their mechanical properties in terms of hardness and toughness were tested. Alumina and aluminum powder mixtures were milled in a planetary ball mill for various times up to 30 h in order to produce Al{sub 2}O{sub 3}-20% Al nanocomposite. The phase composition, morphological and microstructural changes during mechanical milling of the nanocomposite particles were characterized by X-ray diffraction (XRD), transmission electron microscope (TEM), scanning electron microscope (SEM) techniques, respectively. The crystallite size and internal strain were evaluated by XRD patterns using Scherrer methods. A uniform distribution of the Al reinforcement in the Al{sub 2}O{sub 3} matrix was successfully obtained after milling the powders. The results revealed that there was no any sign of phase changes during the milling. The crystal size decreased with the prolongation of milling times, while the internal strain increased. A simple model is presented to illustrate the mechanical alloying of a ductile-brittle component system. A competition between the cold welding mechanism and the fracturing mechanism were found during powder milling and finally the above two mechanisms reached an equilibrium. The maximum relative density was obtained at 1500 Degree-Sign C. The harness of the sintered composite was decreased while the fracture toughness was improved after addition Al into alumina.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A Melting Layer Model for Passive/Active Microwave Remote Sensing Applications: Model Formulation and Comparison with Observations - Part 1

In this study, a 1-D steady-state microphysical model which describes the vertical distribution of melting precipitation particles is developed. The model is driven by the ice-phase precipitation distributions just above the freezing level at applicable gridpoints of "parent" 3-D cloud-resolving model (CRM) simulations. It extends these simulations by providing the number density and meltwater fraction of each particle in finely separated size categories through the melting layer. The depth of the modeled melting layer is primarily determined by the initial material density of the ice-phase precipitation. The radiative properties of melting precipitation at microwave frequencies are calculated based upon different methods for describing the dielectric properties of mixed phase particles. Particle absorption and scattering efficiencies at the Tropical Rainfall Measuring Mission Microwave Imager frequencies (10.65 to 85.5 GHz) are enhanced greatly for relatively small (approx. 0.1) meltwater fractions. The relatively large number of partially-melted particles just below the freezing level in stratiform regions leads to significant microwave absorption, well-exceeding the absorption by rain at the base of the melting layer. Calculated precipitation backscatter efficiencies at the Precipitation Radar frequency (13.8 GHz) increase in proportion to the particle meltwater fraction, leading to a "bright-band" of enhanced radar reflectivities in agreement with previous studies. The radiative properties of the melting layer are determined by the choice of dielectric models and the initial water contents and material densities of the "seeding" ice-phase precipitation particles. Simulated melting layer profiles based upon snow described by the Fabry-Szyrmer core-shell dielectric model and graupel described by the Maxwell-Garnett water matrix dielectric model lead to reasonable agreement with radar-derived melting layer optical depth distributions. Moreover, control profiles that do not contain mixed-phase precipitation particles yield optical depths that are systematically lower than those observed. Therefore, the use of the melting layer model to extend 3-D CRM simulations appears justified, at least until more realistic spectral methods for describing melting precipitation in high-resolution, 3-D CRM's are implemented.

Olson, William S.↗

Scatter in Carbon/Silicon Carbide (C/SiC) Composites Quantified

Carbon-fiber-reinforced silicon carbide matrix (C/SiC) composites processed by chemical vapor infiltration are candidate materials for aerospace thermal structures. Carbon fibers can retain properties at very high temperatures, but they are known to have poor oxidation resistance in adverse, high-temperature environments. Nevertheless, the combination of CVI-SiC matrix with higher stiffness and oxidation resistance, the interfacial coating, and additional surface-seal coating provides the necessary protection to the carbon fibers, and makes the material viable for high-temperature space applications operating under harsh environments. Furthermore, C/SiC composites, like other ceramic matrix composites (CMCs), exhibit graceful non-catastrophic failure because of various inherent energy dissipating mechanisms. The material exhibits nonlinearity in deformation even at very low stress levels. This is the result of the severe matrix microcracking present in the as processed composite because of large differences between the coefficients of thermal expansion of the fiber and the matrix. Utilization of these advanced composites in next generation space vehicles will require innovative structural configurations, updated materials, and refined analyses. Structural safety issues for these vehicles are in direct competition with performance and cost. One would have to quantify the uncertainties associated with the design using formal probabilistic methods. Specifically four fundamental aspects on which analyses are based-- (1) loading conditions, (2) material behavior, (3) geometrical configurations, and (4) structural connections between the composite components and baseline structure--are stochastic in nature. A direct way to formally account for uncertainties is to develop probabilistic structural analysis methods where all participating variables are described by appropriate probability density functions. The present work, however, focuses on analyzing the stochastic material behavior of these advanced composites using formal probabilistic analysis methods. Often, some of the desirable property characteristics that allow composites to offer advantages over conventional structural materials (like tailoring of composite properties) and the complexity are in fact responsible for their greater statistical variability and the requirements for more characterization tests. Composite properties are anisotropic as well, having different properties in different directions. This means that characterization of a property such as stiffness--which will vary greatly depending on the orientation of the fiber relative to the direction of the testing--must be repeated for several different directions and loading conditions. The fabrication process for composites also introduces statistical variations in properties and geometry. A composite part is produced in a number of steps, each of which introduces statistical variability. The matrix is usually produced from a combination of raw materials; and the fiber, which has its own set of properties, is often coated or surface treated, introducing yet another source of variability.

Murthy, Pappu L. N.↗

Efficient real space formalism for hybrid density functionals

We present an efficient real space formalism for hybrid exchange-correlation functionals in generalized Kohn–Sham density functional theory (DFT). In particular, we develop an efficient representation for any function of the real space finite-difference Laplacian matrix by leveraging its Kronecker product structure, thereby enabling the time to solution of associated linear systems to be highly competitive with the fast Fourier transform scheme while not imposing any restrictions on the boundary conditions. We implement this formalism for both the unscreened and range-separated variants of hybrid functionals. We verify its accuracy and efficiency through comparisons with established planewave codes for isolated as well as bulk systems. In particular, we demonstrate up to an order-of-magnitude speedup in time to solution for the real space method. We also apply the framework to study the structure of liquid water using ab initio molecular dynamics, where we find good agreement with the literature. Overall, the current formalism provides an avenue for efficient real-space DFT calculations with hybrid density functionals.

Chemistry↗

Bootstrapping 2d ϕ 4 theory with Hamiltonian truncation data

We combine the methods of Hamiltonian Truncation and the recently proposed generalisation of the S-matrix bootstrap that includes local operators to determine the two-particle scattering amplitude and the two-particle form factor of the stress tensor at s > 0 in the 2d Φ 4 theory. We use the form factor of the stress tensor at s ≤ 0 and its spectral density computed using Lightcone Conformal Truncation (LCT), and inject them into the generalized S-matrix bootstrap set-up. The obtained results for the scattering amplitude and the form factor are fully reliable only in the elastic regime. We independently construct the “pure” S-matrix bootstrap bounds (bootstrap without including matrix elements of local operators), and find that the sinh-Gordon model and its analytic continuation the “staircase model” saturate these bounds. Surprisingly, the Φ 4 two-particle scattering amplitude also very nearly saturates these bounds, and moreover is extremely close to that of the sinh-Gordon/staircase model.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Accurate additive manufacturing of lightweight and elastic carbons using plastic precursors

Despite groundbreaking advances in the additive manufacturing of polymers, metals, and ceramics, scaled and accurate production of structured carbons remains largely underdeveloped. This work reports a simple method to produce complex carbon materials with very low dimensional shrinkage from printed to carbonized state (less than 4%), using commercially available polypropylene precursors and a fused filament fabrication-based process. The control of macrostructural retention is enabled by the inclusion of fiber fillers regardless of the crosslinking degree of the polypropylene matrix, providing a significant advantage to directly control the density, porosity, and mechanical properties of 3D printed carbons. Using the same printed plastic precursors, different mechanical responses of derived carbons can be obtained, notably from stiff to highly compressible. This report harnesses the power of additive manufacturing for producing carbons with accurately controlled structure and properties, while enabling great opportunities for various applications.

36 MATERIALS SCIENCE↗

Reconstruction of apparent orthotropic conductivity tensor image using magnetic resonance electrical impedance tomography

Magnetic resonance electrical impedance tomography visualizes current density and/or conductivity distributions inside an electrically conductive object. Injecting currents into the imaging object along at least two different directions, induced magnetic flux density data can be measured using a magnetic resonance imaging scanner. Without rotating the object inside the scanner, we can measure only one component of the magnetic flux density denoted as B{sub z}. Since the biological tissues such as skeletal muscle and brain white matter show strong anisotropic properties, the reconstruction of anisotropic conductivity tensor is indispensable for the accurate observations in the biological systems. In this paper, we propose a direct method to reconstruct an axial apparent orthotropic conductivity tensor by using multiple B{sub z} data subject to multiple injection currents. To investigate the anisotropic conductivity properties, we first recover the internal current density from the measured B{sub z} data. From the recovered internal current density and the curl-free condition of the electric field, we derive an over-determined matrix system for determining the internal absolute orthotropic conductivity tensor. The over-determined matrix system is designed to use a combination of two loops around each pixel. Numerical simulations and phantom experimental results demonstrate that the proposed algorithm stably determines the orthotropic conductivity tensor.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Barrier height formation in organic blends/metal interfaces: Case of tetrathiafulvalene-tetracyanoquinodimethane/Au(111)

The interface between the tetrathiafulvalene/tetracyanoquinodimethane (TTF-TCNQ) organic blend and the Au(111) metal surface is analyzed by Density Functional Theory calculations, including the effect of the charging energies on the molecule transport gaps. Given the strong donor and acceptor characters of the TTF and TCNQ molecules, respectively, there is a strong intermolecular interaction, with a relatively high charge transfer between the two organic materials, and between the organic layer and the metal surface. We find that the TCNQ LUMO peak is very close to the Fermi level; due to the interaction with the metal surface, the organic molecular levels are broadened, creating an important induced density of interface states (IDIS). We show that the interface energy level alignment is controlled by the charge transfer between TTF, TCNQ, and Au, and by the molecular dipoles created in the molecules because of their deformations when adsorbed on Au(111). A generalization of the Unified-IDIS model, to explain how the interface energy levels alignment is achieved for the case of this blended donor/acceptor organic layer, is presented by introducing matrix equations associated with the Charge Neutrality Levels of both organic materials and with their intermixed screening properties.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Quantitative Uncertainty Analysis for Diffraction

Current software uses a least squares statistical method. For each model parameter, current software yields a point estimate and estimated standard error for each point estimate. In contrast, we have employed a Bayesian approach with sampling algorithm to provide a full posterior probability density function on each model parameter. This yields a different type of solution (probabilistic) with better fidelity. The Quantitative Uncertainty Analysis for Diffraction (QUAD) package returns a matrix of affect model parameters as a function of iteration through the MCMC analysis. QUAD is a Python code that allows for analysis of X-ray and neutron diffraction data to infer the structure of materials with quantifiable certainty. The code uses Bayesian statistics and Markov chain sampling algorithms to create posterior probability distributions on all material structure parameters that are modeled by a researcher. The statistical packages use the open source GSASII project as a library to calculate model diffraction patterns.

Fancher, ChristopherM↗

A Review of Tungsten Heavy Alloy Utilization in Isotope Transport Containers - 13380

A common requirement for radioisotope transport containers is that they provide both durable and efficient shielding of penetrating gamma radiation. This is the case for transport of both spent nuclear fuel as well as intentionally created radioisotopes for medical or other uses. Tungsten heavy alloy (WHA) provides a unique engineering property set for such shielding - easily surpassing more commonly used lead alloys in both strength and attenuation. This family of alloys contains typically 90-98 wt.% W in combination with transition metals such as Ni and Fe. WHA is manufactured in near net shape blanks by liquid phase sintering of compacted powder shapes to full metallurgical density parts. This powder metallurgy approach is described in its ability to provide excellent material utilization and affords efficient manufacturing of various shapes required for gamma shields or collimators. WHAs offer very high density (approaching 19 g/cc) in combination with relatively high thermal conductivity, low thermal expansion, ambient corrosion resistance, and can be provided with mechanical properties comparable to many medium carbon steels. As such, they can be machined to complex, damage resistant geometries using common metal cutting tools and methods. WHA additionally provides a lower toxicity alternative to Pb- or U-based gamma shielding. Given the specialty nature of WHA, specific metallurgical characteristics are reviewed to assist shielding designers who may otherwise encounter difficulties locating important alloy selection and fabrication details. Contained within this materials and applications overview are guidelines for WHA component design, alloy selection, and practical machining, finishing, and assembly considerations. The microstructure of WHA is that of a metal matrix composite. This factor has specific implications in the design of components for stress service as well as their protection in the presence of electrolytes. WHA is also discussed in the broader context of materials compatibility, as it is rarely used in isolated monolithic shapes. Alloy selection for new applications is often made primarily on the basis of density. An alternative strategy to this selection approach is presented which proposes that mechanical requirements for a given shielding use be the primary selection criterion over density. Standard commercial grades of WHA for radiation shielding are defined by specifications such as AMS-T-21014 and ASTM B777. These specifications define 4 density classes of as-sintered WHA. Compositional options, as well as post-sinter processing of WHA, are discussed for shielding components that must exhibit higher levels of ductility or very low magnetic permeability. In addition to the mechanical advantages over Pb-based shielding, the higher linear attenuation of energetic photons for various grades of WHA (as calculated by the NIST XCOM routine) are presented for selected photon energies of interest to illustrate the shield volume reduction generally possible through the use of tungsten-based shielding. While W provides inefficient attenuation of neutrons in a mixed radiation environment, its secondary role in shielding gamma radiation produced as a result of neutron capture is also described. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Sol–gel synthesis of long-lasting phosphors CdSiO{sub 3}: Mn{sup 2+}, RE{sup 3+} (RE = Tb, Eu, Nd) and luminescence mechanism research

Highlights: ► New long-lasting CdSiO{sub 3}: Mn{sup 2+}, RE{sup 3+} (RE = Tb, Eu, Nd) phosphors were synthesized by a sol–gel method. ► The afterglow performance of the CdSiO{sub 3}: Mn{sup 2+}, Eu{sup 3+} phosphor was the best. ► The role of RE{sup 3+} co-doped into the CdSiO{sub 3}: Mn{sup 2+} matrix was discussed in this paper. -- Abstract: Mn{sup 2+} and RE{sup 3+} (RE = Tb, Eu, Nd) co-doped CdSiO{sub 3} orange phosphors were prepared at 1050 °C by a sol–gel method. The phase and crystallinity of the synthesized materials were investigated by X-ray powder diffraction (XRD) and scanning electron microscopy (SEM). The luminescence characteristics were analyzed using photoluminescence (PL) spectra, afterglow decay curves, long-lasting phosphorescence spectra, and thermoluminescence (TL) spectra. Due to the difference in co-doped rare earth ionic radii, it varied greatly in trap density and trap depth caused by the different defects deriving from RE{sup 3+} ions co-doping into the CdSiO{sub 3}: Mn{sup 2+} host. The afterglow intensity and time for these samples increased as follows: CdSiO{sub 3}: Mn{sup 2+}0.2%, Nd{sup 3+}0.8% < CdSiO{sub 3}: Mn{sup 2+}0.4%, Tb{sup 3+}0.8% < CdSiO{sub 3}: Mn{sup 2+}0.4%, Eu{sup 3+}0.3%. CdSiO{sub 3}: Mn{sup 2+}0.4%, Eu{sup 3+}0.3% had the best afterglow properties, which could be due to the proper traps formed by Eu{sup 3+} ions co-doping into the host. The role of RE{sup 3+} co-doped into the CdSiO{sub 3}: Mn{sup 2+} matrix and the possible long-lasting phosphorescence process was also discussed in this paper.

36 MATERIALS SCIENCE↗

Quantitative imaging of trace B in Si and SiO{sub 2} using ToF-SIMS

Changes in the oxidation state of an element can result in significant changes in the ionization efficiency and hence signal intensity during secondary ion mass spectrometry (SIMS) analysis; this is referred to as the SIMS matrix effect [Secondary Ion Mass Spectrometry: A Practical Handbook for Depth Profiling and Bulk Impurity Analysis, edited by R. G. Wilson, F. A. Stevie, and C. W. Magee (Wiley, New York, 1990)]. The SIMS matrix effect complicates quantitative analysis. Quantification of SIMS data requires the determination of relative sensitivity factors (RSFs), which can be used to convert the as measured intensity into concentration units [Secondary Ion Mass Spectrometry: A Practical Handbook for Depth Profiling and Bulk Impurity Analysis, edited by R. G. Wilson, F. A. Stevie, and C. W. Magee (Wiley, New York, 1990)]. In this manuscript, the authors report both: RSFs which were determined for quantification of B in Si and SiO{sub 2} matrices using a dual beam time of flight secondary ion mass spectrometry (ToF-SIMS) instrument and the protocol they are using to provide quantitative ToF-SIMS images and line scan traces. The authors also compare RSF values that were determined using oxygen and Ar ion beams for erosion, discuss the problems that can be encountered when bulk calibration samples are used to determine RSFs, and remind the reader that errors in molecular details of the matrix (density, volume, etc.) that are used to convert from atoms/cm{sup 3} to other concentration units will propagate into errors in the determined concentrations.

42 ENGINEERING↗

Response theory of the ergodic many-body delocalized phase: Keldysh Finkel’stein sigma models and the 10-fold way

We derive the finite temperature Keldysh response theory for interacting fermions in the presence of quenched short-ranged disorder, as applicable to any of the 10 Altland–Zirnbauer classes in an Anderson delocalized phase with at least a U(1) continuous symmetry. In this formulation of the interacting Finkel’stein nonlinear sigma model, the statistics of one-body wave functions are encoded by the constrained matrix field, while physical correlations follow from the hydrodynamic density or spin response field, which decouples the interactions. Integrating out the matrix field first, we obtain weak (anti) localization and Altshuler–Aronov quantum conductance corrections from the hydrodynamic response function. This procedure automatically incorporates the correct infrared cutoff physics, and in particular gives the Altshuler–Aronov–Khmelnitsky (AAK) equations for dephasing of weak (anti)localization due to electron–electron collisions. We explicate the method by deriving known quantumcorrections in two dimensions for the symplectic metal class AII, as well as the spin-SU(2) invariant superconductor classes C and CI. We show that quantum conductance corrections due to the special modes at zero energy in nonstandard classes are automatically cut off by temperature, as previously expected, while the Wigner–Dyson class Cooperon modes that persist to all energies are cut by dephasing. We also show that for short-ranged interactions, the standard self-consistent solution for the dephasing rate is equivalent to a particular summation of diagrams via the self-consistent Born approximation. This should be compared to the corresponding AAK solution for long-ranged Coulomb interactions, which exploits the Markovian noise correlations induced by thermal fluctuations of the electromagnetic field. We discuss prospects for exploring the many-body localization transition as a dephasing catastrophe in short-range interacting models, as encountered by approaching from the ergodic side. - Highlights: • The Keldysh response theory for disordered, interacting fermions is derived. • The theory applies to any Altland–Zirnbauer class with at least a U(1) symmetry. • The correct infrared physics at finite temperature is automatically incorporated. • The Altshuler–Aronov–Khmelnitsky equations for dephasing are systematically derived. • A strategy to explore the MBL transition as a dephasing catastrophe is proposed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Many-body localization as a large family of localized ground states

Many-body localization (MBL) addresses the absence of thermalization in interacting quantum systems, with nonergodic high-energy eigenstates behaving as ground states, only area-law entangled. However, computing highly excited many-body eigenstates using exact methods is very challenging. Instead, we show that one can address high-energy MBL physics using ground-state methods, which are much more amenable to many efficient algorithms. Here, we find that a localized many-body ground state of a given interacting disordered Hamiltonian H 0 is a very good approximation for a high-energy eigenstate of a parent Hamiltonian, close to H 0 but more disordered. Lastly, this construction relies on computing the covariance matrix, easily achieved using density-matrix renormalization group for disordered Heisenberg chains up to L=256 sites.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Study of Thermal and Coherent A{sub 1g} Phonons in Bismuth Telluride

Totally symmetric A{sub 1g} phonons are studied for the equilibrium and coherent states of a Bi{sub 2}Te{sub 3} lattice. Equilibrium phonons were investigated in the frequency domain by the method of spontaneous Raman scattering, whereas coherent phonons were studied by the method of active femtosecond spectroscopy in the time domain. In the latter case, femtosecond laser pulses were used both for generating and detecting coherent A{sub 1g} phonons having a well-defined phase allowing the selective optical control of the lattice dynamics. A comparison of the results obtained in the frequency and time domains suggests that diagonal and nondiagonal elements of the density matrix of lattice excitations relax with the same characteristic time to the equilibrium and zero values, respectively.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗