Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “relaxation”

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

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

At least 1,495 records · Page 83

A collision operator for describing dissipation in noncanonical phase space

The phase space of a noncanonical Hamiltonian system is partially inaccessible due to dynamical constraints (Casimir invariants) arising from the kernel of the Poisson tensor. When an ensemble of noncanonical Hamiltonian systems is allowed to interact, dissipative processes eventually break the phase space constraints, resulting in a thermodynamic equilibrium described by a Maxwell–Boltzmann distribution. However, the time scale required to reach Maxwell–Boltzmann statistics is often much longer than the time scale over which a given system achieves a state of thermal equilibrium. Examples include diffusion in rigid mechanical systems, as well as collisionless relaxation in magnetized plasmas and stellar systems, where the interval between binary Coulomb or gravitational collisions can be longer than the time scale over which stable structures are self-organized. Here, we focus on self-organizing phenomena over spacetime scales such that particle interactions respect the noncanonical Hamiltonian structure, but yet act to create a state of thermodynamic equilibrium. We derive a collision operator for general noncanonical Hamiltonian systems, applicable to fast, localized interactions. This collision operator depends on the interaction exchanged by colliding particles and on the Poisson tensor encoding the noncanonical phase space structure, is consistent with entropy growth and conservation of particle number and energy, preserves the interior Casimir invariants, reduces to the Landau collision operator in the limit of grazing binary Coulomb collisions in canonical phase space, and exhibits a metriplectic structure. We further show how thermodynamic equilibria depart from Maxwell–Boltzmann statistics due to the noncanonical phase space structure, and how self-organization and collisionless relaxation in magnetized plasmas and stellar systems can be described through the derived collision operator.

Boltzmann equation↗

A unified non-equilibrium phase change model for injection flow modeling

The homogenous relaxation model (HRM) is one of the most widely used models to describe the liquid- gas phase transition. However, in its original formulation, it is unable to handle multispecies vapor-liquid equilibrium (VLE), which limits its applicability to single-component fluids. In this work, a unified non-equilibrium phase change model that considers the VLE of multicomponent mixtures is proposed building upon the HRM's structure. A time factor is introduced to mimic the effect of different phase change timescales due to different mechanisms, e.g., cavitation, flash-boiling, and evaporation. Here to assess the model's performance, computational fluid dynamics simulations of the internal and near-nozzle injection flow of the Engine Combustion Network's Spray G injector were performed using the nine-component PACE-20 fuel with both the unified model and the original HRM. The predicted fuel density in the near-nozzle region matched well with X-ray tomography measurements. The simulation results indicated that, whereas the HRM failed to capture the vaporization due to convective mixing between the fuel and ambient gas, the unified model performed well in predicting the mixing-driven vaporization and the corresponding evaporative cooling. Further comparisons using the nine-component fuel formula and a single-component fuel surrogate demonstrated the unified model's ability to predict preferential vaporization, which affects the predictions of local mixture composition and rate of vaporization. Finally, it is shown that the unified model is capable of representing multiple phase change mechanisms, and the relaxation time factor plays an important role in determining the degree of phase change due to the different mechanisms.

33 ADVANCED PROPULSION SYSTEMS↗

Residual elastic strain evolution due to thermal cycling of a ceramic-metal composite (WC-Cu) via high energy X-ray diffraction and analytical modeling

Residual stress, when superimposed with in-service loading, can significantly reduce the lifetime and performance of a component. Ceramic-metal composites are susceptible to residual stresses due to the thermal expansion mismatch of the ceramic and metallic phases. The WC-Cu composite explored in the present study provides a promising combination of thermal conductivity and strength properties, while exhibiting counterintuitive improvements in strength and ductility after thermal cycling. Further, this work quantifies the evolution of the residual elastic strains as a result of processing and cyclic thermal loading in a co-continuous WC-Cu composite through experimental high energy X-ray diffraction and kinetics-based modeling. Both analyses indicate that processing-induced residual tensile stress in the copper phase is relieved upon subsequent thermal cycling, with kinetics modeling revealing the cyclic-dependent nature of the active power-law creep mechanisms. The results indicate that, through stress relaxation, this material system maintains structural stability during thermal cycling. The illustrated kinetics of relaxation can inform general material processors and designers of ceramic-metal composites to minimize detrimental residual stress and improve performance of these material systems.

36 MATERIALS SCIENCE↗

Structure-performance relationships in lignin-based transesterification vitrimers: The role of lignin structural features

Lignin has been hailed as an ideal renewable alternative for petrochemical-based prepolymers in material synthesis for a sustainable and circular economy, due to its abundant aromatic network and high carbon content. However, the properties and performance of lignin-derived macromolecules are strongly influenced by the lignin itself. While numerous studies have explored the impact of lignin content on the thermomechanical performance of lignin-based vitrimers, literature on how the inherent structural features of lignin affect these properties is scanty. In this study, hardwood organosolv lignin was fractionated in ethyl acetate, ethanol, and acetone to obtain lignin fractions with varying structural characteristics. These fractions were then modified through carboxylation and crosslinked with epoxidized soybean oil (ESO) at a hydroxyl to epoxy group ratio of 1:1 to produce lignin-based transesterification vitrimers (LVs). The thermal properties (i.e. glass transition temperature and thermal stability), tensile strength, storage modulus, and stress relaxation behavior of the LVs were studied and carefully related to the structural features of lignin. The results revealed a positive relationship between strong hydroxyl content in modified lignin and the tensile strength (5.10–9.71 MPa), storage modulus (1099.4 – 1372.8 MPa), crosslinking density, and stress relaxation of the LVs. Additionally, both the storage modulus and tensile strength exhibited a positive relationship with the ratio of rigid linkages in modified lignin, while lignin molecular weight was found to significantly impact the thermal properties of LVs (i.e Tg and thermal stability). This study not only highlights the valorization of lignin in vitrimer synthesis but also provide insights for designing lignin-based materials with tailored properties for specific applications.

Bio-based polymer↗

Magnetic anomalies associated with domain wall freezing and coupled electron hopping in magnetite nanorods

Magnetite has fascinated researchers for decades, due to its wide range of applications from spintronics to biomedicine. Despite a large body of works aimed at its magnetic properties, no consensus has been reached on the physical origin of the low temperature magnetic anomalies observed in magnetite. Although, a lot of work has been done in studying magnetite nanoparticles, but studies on the low temperature anomalies in those nanoparticles still remains unresearched. We report on the observation of the low temperature magnetic anomalies in highly crystalline, stoichiometric Fe 3 O 4 nanorods and relate them to the coupled electron hopping relaxation process and domain wall motion. Both DC and AC susceptibility show the presence of a hump around 35 K, which is associated with the relaxation of Fe +2 extra electrons. Radiofrequency transverse susceptibility (TS) measurements indicated a noticeable increase in anisotropy field below ~ 25 K, which is attributed to the rearrangement of Fe +2 electrons in the octahedral sites. TS experiments also revealed the domain wall freezing below ~35 K. Our combined DC, AC and TS susceptibility studies shed light on the complex nature of the lowtemperature magnetic behavior in nanostructured magnetite.

36 MATERIALS SCIENCE↗

The effect of annealing on optical transmittance and structure of ZLANI fluorozirconate glass thin films

In this work, we report the effect of thermal annealing on structure, composition, optical transmittance and thickness of a novel fluorozirconate glass (ZLANI) containing Zr, La, Al, Na and In fluorides. In this work, pulsed laser depo- sition was used to grow thin films of ZLANI, and thermal annealing at different temperatures was performed on the films. Annealing did not change the composition, but a clear structural evolution of the ZLANI glass was observed by transmission electron microscopy (TEM), showing that we can control microstructure independent of composition. An increase in transmittance after the film was subject to a 100 °C thermal anneal was ascribed to the removal of defects and structural relaxation in the amorphous state. Following an anneal of 200 °C, the transmittance decreases due to heterogeneous formation of crystalline nuclei and changes in the local bonding. After the final annealing at 300 °C, a wider-scale crystallization occurred, with some major crystal phases formed as Zr 2 F 8 (H 2 O) 6 and ZrO 2 , which alters the shape of the transmittance curve. The crystalline content of the crystal phases that form in the annealed films was quantified using hollow cone dark field TEM imaging. The 100 °C or 200 °C annealing decreases the film thickness by inducing structural relaxation and densification of the amorphous films, while the thickness increase of the 300 °C annealed film resulted from the formed large crystals. These results provide insights for the design of multilayer nanocomposites with a ZLANI glass matrix, which have potential applications as up-/down-conversion luminescent materials and X-ray storage phosphors.

36 MATERIALS SCIENCE↗

Chromium enhances the mechanical performance of 3d transition metal high entropy alloys

A comprehensive examination of the compositional effects on the deformation behavior of high entropy alloys (HEAs) was conducted through nanoindentation, indentation creep, and stress relaxation experiments at ambient temperature. The alloys investigated included NiCoFe, NiCoCr, NiCoCrFe, and NiCoCrFeMn, all of which underwent identical thermomechanical processing. From the experimental results, NiCoCr exhibited the highest maximum shear stress for dislocation nucleation (13 GPa) and nanoindentation hardness (3.8 GPa) among the alloys tested. It also showed the lowest stress sensitivity (0.02) and activation volume (3–4 b 3 ) under creep and stress relaxation conditions, indicative of its superior plastic flow properties. Density functional theory (DFT) calculations further revealed an uneven charge density and larger bonding-length variation in NiCoCr due to the presence of Cr, which increased lattice distortion, impeding dislocation movement. Furthermore, in this study, post-deformation microscopy using transmission electron microscopy (TEM) revealed a high density of stacking faults and twins in NiCoCr that effectively enhanced its creep strength. The results highlight the significance of specific elemental effects, particularly from Cr in this study, over configurational entropy effect (i.e., compositional complexity) in governing the deformation microstructure and mechanical properties of HEAs. Finally, these insights will be instrumental for the design and optimization of advanced alloys for load-bearing applications.

36 - MATERIALS SCIENCE↗

Origins of heat transport anisotropy in MoTe2 and other bulk van der Waals materials

Knowledge of how heat flows anisotropically in van der Waals (vdW) materials is crucial for thermal management of emerging 2D materials devices and design of novel anisotropic thermoelectric materials. Despite the importance, anisotropic heat transport in vdW materials is yet to be systematically studied and is often presumably attributed to anisotropic speeds of sound in vdW materials due to soft interlayer bonding relative to covalent in-plane networks of atoms. In this work, we investigate the origins of the anisotropic heat transport in vdW materials, through time-domain thermoreflectance (TDTR) measurements and first-principles calculations of anisotropic thermal conductivity of three different phases of MoTe 2 . MoTe 2 is ideal for the study due to its weak anisotropy in the speeds of sound. We find that even when the speeds of sound are roughly isotropic, the measured thermal conductivity of MoTe 2 along the c-axis is 5–8 times lower than that along the in-plane axes. We derive meaningful characteristic heat capacity, phonon group velocity, and relaxation times from our first principles calculations for selected vdW materials (MoTe 2 , BP, h-BN, and MoS 2 ), to assess the contributions of these factors to the anisotropic heat transport. Interestingly, we find that the main contributor to the heat transport anisotropy in vdW materials is anisotropy in heat capacity of the dominant heat-carrying phonon modes in different directions, which originates from anisotropic optical phonon dispersion and disparity in the frequency of heat-carrying phonons in different directions. Further, the discrepancy in frequency of the heat-carrying phonons also leads to ~2 times larger average relaxation times in the cross-plane direction, and partially explains the apparent dependence of the anisotropic heat transport on the anisotropic speeds of sound. This work provides insight into understanding of the anisotropic heat transport in vdW materials.

36 MATERIALS SCIENCE↗

Dynamics of Glyceline and Interactions of Constituents: A Multitechnique NMR Study

In this work, the dynamics of the organic components of the deep eutectic solvent (DES) glyceline are analyzed using an array of complementary nuclear magnetic resonance (NMR) methods. Fast-Field Cycling 1 H relaxometry, Pulsed Field Gradient diffusion, Nuclear Overhauser Effect Spectroscopy (NOESY), 13 C NMR relaxation and pressure dependent NMR experiments are deployed to sample a range of frequencies and modes of motion of the glycerol and choline components of the DES. Generally, translational and rotational diffusion of glycerol are more rapid than those of choline while short range rotational motions observed from 13 C relaxation indicate slow local motion of glycerol at low choline chloride (ChCl) content. The rates of glycerol and choline local motions become more similar at higher ChCl. This result taken together with pressure dependent NMR studies show that the addition of ChCl makes it easier to disrupt glycerol packing. Finally, a relatively slow hydroxyl H-exchange process between glycerol and choline protons is deduced from the data. Consistent with this, NOESY results indicate relatively little direct H-bonding between glycerol and choline. These results suggest that the glycerol H-bonding network is disrupted as choline is added, but primarily in regions where there is intimate mixing of the two components. Thus, the local dynamics of most of the glycerol, resembles that of pure glycerol until substantial choline chloride is present.

13C NMR relaxation↗

Plasmon-Induced Hot-Carrier Excited-State Dynamics in Plasmonic Semiconductor Nanocrystals

The variance of carrier relaxation pathways for WO 3–x plasmonic semiconductor nanocrystals (PSNCs) is monitored by transient absorption spectroscopy following excitation of the localized surface plasmon resonance (LSPR) versus the optical band gap (E g,opt ). Excitation of the LSPR leads to efficient hot carrier population above the Fermi level in WO 3–x via Landau damping, in analogy to noble metal LSPR relaxation mechanisms. Hot carrier depopulation occurs on the femtosecond timescale, observed as the concomitant recovery of an LSPR bleach with the appearance of discrete interband and intraband photoinduced absorption features. Further, by comparison, the direct excitation of E g,opt results in trion recombination at donor–acceptor sites within the WO 3–x NC, consistent with exciton decay dynamics observed for typical wide-band-gap semiconductor NCs. From the analysis of pump power dependency data, a hot-carrier electron–phonon coupling constant of 1.47 × 10 11 J K –1 s –1 cm –3 is extracted. The direct comparison of the decay dynamics following E g,opt versus LSPR excitation confirms that the observed plasmon in trioctylphosphine oxide passivated, spherical WO 3–x is a resonance state in which hot carriers are generated only from excitation on resonance with the LSPR frequency. This study on WO 3–x PSNCs provides a toolset that can be used to evaluate the role of hot carriers following LSPR excitation of n-type, plasmonic transparent conducting oxide NCs, where enhancement of photocatalysis, photovoltaic performance, and optical enhancement has been reported.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Deconvolution of Water-Splitting on the Triple-Conducting Ruddlesden–Popper-Phase Anode for Protonic Ceramic Electrolysis Cells

Triple-conducting materials have been proved to improve the performance of popular protonic ceramic electrolysis cells. However, partially because of the complexity of the water splitting reaction involving three charge carriers, that is, oxygen (O 2– ), proton (H + ), and electron (e – ), the triple-conducting reaction mechanism was not clear, and the reaction conducting pathways have seldom been addressed. In this study, the triple conducting Ruddlesden–Popper phase Pr 1.75 Ba 0.25 NiO 4+δ as an anode on the BaCe 0.7 Zr 0.1 Y 0.1 Yb 0.1 O 3–δ electrolyte was fabricated and its electroresponses were characterized by electrochemical impedance spectroscopy with various atmospheres and temperatures. The impedance spectra are deconvoluted by means of the distribution of the relaxation time method. The surface exchange rate and chemical diffusivity of H + and O 2– are characterized by electrical conductivity relaxation. The physical locations of electrochemical processes are also identified by atomic layer deposition with a surface inhibitor. A microkinetics model is proposed toward conductivities, triple-conducting pathways, reactant dependency, surface exchange and bulk diffusion capabilities, and other relevant properties. Lastly, the rate-limiting steps and suggestions for further improvement of electrode performance are presented.

25 ENERGY STORAGE↗

Proton surface exchange kinetics of perovskite triple conducting thin films for protonic ceramic electrolysis cells: BaPr 0.9 Y 0.1 O 3–δ (BPY) vs. Ba 1–x Co 0.4 Fe 0.4 Zr 0.1 Y 0.1 O 3–δ (BCFZY)

Protonic ceramic electrolysis cells (PCECs) are an attractive green H 2 production technology, given their intermediate-temperature operating range and ability to produce dry H 2 . However, PCECs will benefit from development of more efficient and durable “triple conducting” anodes where steam is split, H incorporated, and oxygen evolved. In this work, we evaluated the kinetics of the steam-splitting/H incorporation reaction on BaPr 0.9 Y 0.1 O 3–δ (BPY) in comparison to the benchmark Ba 1–x Co 0.4 Fe 0.4 Z r0.1 Y 0.1 O 3–δ (BCFZY) composition, replacing most of the transition metal elements (Co, Fe, Zr) with the lanthanide Pr. We prepared geometrically well-defined perovskite BPY and BCFZY thin films by pulsed laser deposition and performed simultaneous optical transmission relaxation and electrical conductivity relaxation measurements at 400–500 °C in 0.21 atm O 2 during switching of the steam partial pressure to isolate and compare their proton surface exchange coefficients (k). The k values of BPY were comparable to those of BCFZY and more stable over time. According to angle-resolved XPS and STEM-EDS mapping of FIB cross-sections, the surface of BPY exhibited Ba enrichment, Pr deficiency, and Si contamination. In contrast, BCFZY exhibited Ba deficiency throughout, no obvious surface segregation, and less Si contamination. The Ba segregation on the BPY film appears to have promoted steam splitting/H incorporation kinetics even though the more basic surface reacted with the acidic environmental SiO x H y . Faster kinetics observed on stoichiometric BCFZY vs. Ba-deficient BCFZY confirmed the benefit of a high A-site Ba concentration. This result contrasts with most work on perovskites applied in solid oxide electrolysis cell anodes, in which A-site segregation is considered deleterious for surface reaction kinetics.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Triboplasma assisted chemical conversion in granular systems: a semi-analytic model

We present a semi-analytic model for a novel plasma-assisted chemical conversion pathway using triboplasmas generated in granular flows. Triboelectric charge relaxation is a well known phenomena where the potential generated from contact charging of particles exceeds the breakdown voltage of the background gas. In this work, we extend the triboelectric charge relaxation theory to include non equilibrium plasma energy and particle balance equations to predict the formation of dissociated and excited species that act as precursors to chemical conversion, for example in plasma-assisted ammonia synthesis. Our example case study with nitrogen background gas and teflon/aluminum tribomaterial system yielded high excited nitrogen species densities per collision that are comparable to current plasma-assisted conversion pathways. We also present a regime diagram for various gases where Paschen breakdown parameters are used to determine whether triboplasmas can be formed for a given effective work-function difference between two materials. Our sensitivity studies indicate particle velocity, particle radius, solids fraction and space charge effects play a critical role in overall plasma densities and excited species production.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Planar XY magnetic glass state in the Gd 2 ScNbO 7 pyrochlore

Here a spin glass system with emergent planar ordered spin clusters is investigated. The mixed B-site pyrochlore Gd 2 ScNbO 7 has been synthesized and characterized through a variety of techniques, including x-ray diffraction, magnetic susceptibility, muon spin relaxation, heat capacity and neutron scattering. Despite a Curie–Weiss temperature of -3.93(3) K, indicating net antiferromagnetic interactions, no signs of long ranged magnetic ordering are found down to T = 0.3 K. Instead, a disordered magnetic state emerges with a small correlation length of 2.1(1) Å of single tetrahedra. A Reverse Monte Carlo analysis of the polarized neutron scattering data reveals short-range antiferromagnetic order with emergent XY spin ordering similar to the parent pyrochlore compounds. Muon spin relaxation, and AC susceptibility measurements confirm that the magnetization condenses into a glass, with 10 % of the potential entropy missing in the specific heat. This magnetic ground state is similar to what is observed in Gd 2 Sn 2 O 7 just above the ordering temperature, without the eventual long-range ordering at low temperature.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

The properties of the nitrogen-vacancy center in milled chemical vapor deposition nanodiamonds

Fluorescent nanodiamonds (FNDs) containing negatively charged nitrogen-vacancy (NV − ) centers are vital for many emerging quantum sensing applications from magnetometry to intracellular sensing in biology. However, developing a scalable fabrication method for FNDs hosting color centers with consistent bulk-like photoluminescence (PL) and spin coherence properties remains a highly desired but unrealized goal. Here, we investigate optimized ball milling of single-crystal diamonds produced via chemical vapor deposition (CVD) and containing 2 ppm of substitutional nitrogen and 0.3 ppm of NV − to achieve this goal. The NV charge state, PL lifetime, and spin properties of bulk CVD diamond samples are directly compared to milled CVD FNDs and commercial high-pressure high-temperature (HPHT) FNDs. We find that on average, the relative contribution of the NV − charge state to the total NV PL is lower and the NV PL lifetime is longer in CVD FNDs compared to HPHT FNDs, both likely due to the lower N s 0 concentration in CVD FNDs. The CVD bulk and CVD FNDs on average show similar average T 1 spin relaxation times of 3.2 ± 0.7 ms and 4.7 ± 1.6 ms, respectively, compared to 0.17 ± 0.01 ms for commercial HPHT FNDs. Our results demonstrate that ball milling of CVD diamonds enables the large-scale fabrication of NV ensembles in FNDs with bulk-like T 1 spin relaxation properties.

chemical vapor deposition↗

Hyperbolic Reformulation Approach to Enable Efficient Simulation of Groundwater Flow and Reactive Transport

We apply Cattaneo's relaxation approach to the one-dimensional coupled Boussinesq groundwater flow and advection-diffusion-reaction equations, commonly used in engineering applications to simulate contaminant transport in the subsurface. The diffusion-type governing equations are reformulated as a hyperbolic system, augmented by an equation that can be interpreted as a momentum balance. The hyperbolization enables an efficient unified computation of the primary variable and its gradients, for example piezometric head and unit discharge in the Boussinesq equation. An augmented Roe scheme is used to solve the hyperbolic system. The hyperbolized system of equations is studied in a set of steady state and transient test cases with idealized geometry. These test cases confirm the equivalence of the hyperbolic system to its original formulation. The larger time step size of the hyperbolic equation is verified theoretically by means of a stability analysis and numerically in the test cases. Finally, a reach-scale application of flow and transport across a river meander is considered. This application case shows that the performance of the hyperbolic relaxation approach holds for more realistic groundwater flow and transport problems, relevant to water resources management.

54 ENVIRONMENTAL SCIENCES↗

Origin of the anomalously low Raman exponents in single molecule magnets

The Raman exponent of single molecular magnetic relaxation may take various unexpected values because of rich phonon spectra and spin-phonon coupling of molecular crystals. Here, we systematically examine the origins of different abnormalities and clarify misunderstandings of the past, particularly the appropriateness of the fitting procedures for the exponents. We find that exponential laws raised by optical phonons can yield spurious power laws with low exponents. This observation indicates long-standing misunderstandings of the origins of low Raman exponents in a large bulk of single molecule magnets. Resulting from spin-lattice coupling with optical modes, the presence of these exponents suggests the importance of the local dynamical environment for the magnetic relaxation in this regime.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Magnetization reversal driven by electron localization-delocalization crossover in the inverse spinel Co 2 VO 4

Neutron diffraction, magnetization, and muon spin relaxation measurements, supplemented by density functional theory (DFT) calculations are employed to unravel temperature-driven magnetization reversal in inverse spinel Co 2 VO 4 . All measurements show a second-order magnetic phase transition at T C = 168 K to a collinear ferrimagnetic phase. Neutron diffraction measurements reveal two antiparallel ferromagnetic (FM) sublattices, belonging to magnetic ions on two distinct crystal lattice sites, where the relative balance between the two sublattices determine the net FM moment in the unit cell. As the evolution of the ordered moment with temperature differs between the two sublattices, the net magnetic moment reaches a maximum at T NC = 138 K and reverses its sign at T MR = 65 K. The DFT results suggest that the underlying microscopic mechanism for the reversal is a delocalization of the unfilled 3d -shell electrons on one sublattice just below T C , followed by a gradual localization as the temperature is lowered. Furthermore, this delocalized-localized crossover is supported by muon spectroscopy results, as strong T 1 relaxation observed below T C indicates fluctuating internal fields.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗