Search NASA⌕ Search

SEARCH · Search NASA

Results for “SHEAR STRAIN”

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 37 records · Page 2

Predicting the proximity to macroscopic failure using local strain populations from dynamic in situ X-ray tomography triaxial compression experiments on rocks

Predicting the proximity of large-scale dynamic failure is a critical concern in the engineering and geophysical sciences. Here we use evolving contractive, dilatational, and shear strain deformation preceding failure in dynamic X-ray tomography experiments to examine which strain components best predict the proximity to failure. We develop machine learning models to predict the proximity to failure using time series of three-dimensional local incremental strain tensor fields acquired in rock deformation experiments under stress conditions of the upper crust. Three-dimensional scans acquired in situ throughout triaxial compression experiments provide a distribution of density contrasts from which we estimate the three-dimensional incremental strain that accumulates between each scan acquisition. Training machine learning models on multiple experiments of six rock types provides suites of feature importance that indicate the predictive power of each feature. Comparing the average importance of groups of features that include information about each strain component quantifies the ability of the contractive, dilatational and shear strain to predict the proximity of macroscopic failure. A total of 24 models of four machine learning algorithms with six rock types indicate that 1) the dilatational strain provides the best predictive power of the strain components, and 2) the intermediate values (25th-75th percentile) of the strain population provide the best predictive power of the statistics of the strain populations. In addition, the success of the predictions of models trained on one rock type and tested on other rock types quantifies the similarities and differences of the precursory strain accumulation process in the six rock types. These similarities suggest the potential existence of a unified theory of brittle rock deformation for a range of rock types.

58 GEOSCIENCES↗

Survival estimates across five life stages of redfin ( Perca fluviatilis ) exposed to simulated pumped-storage hydropower stressors

The global prevalence of pumped-storage hydropower (PSH) is expected to grow exponentially as countries transition to renewable energy sources. Compared to conventional hydropower, little is currently known regarding PSH impacts on aquatic biota. This study estimated the survival of five life stages (egg, two larval stages, juvenile and adult) of redfin (European) perch (Perca fluviatilis) following passage through a PSH facility during the pumping phase. This was achieved by simulating the individual stressors expected to occur during passage through a 2000-MW PSH facility using laboratory-simulated (shear strain and extreme compression) and modelling (blade strike, BS) approaches. Our results indicate that redfin could survive the shear, pressure and BS stressors expected within the PSH facility, but impacts varied among life stages. Juvenile survival was >70% across all shear strain rates, while the survival of eggs and larvae declined markedly as strain rate increased. All life stages had high survival when exposed to rapid compression and BS. The high survival of redfin to the stressors tested suggests the PSH facility could facilitate the passage of redfin during the pumping phase from the lower to the higher elevation reservoir. This outcome would be welcomed in situations where the species is native, but could have adverse implications for the conservation of native biota where the species is considered a pest.

59 BASIC BIOLOGICAL SCIENCES↗

Intensity measure adequacy assessment for nonlinear site response using Information Theory

Nonlinear site response modeling is a crucial aspect of Probabilistic Seismic Hazard Analysis. Site amplification models routinely rely on a rock intensity measure to characterize the strength of the bedrock motion. However, the adequacy of such intensity measures towards predicting amplifications across the oscillator period range has not been investigated in the literature. This paper analyzes the adequacy of rock intensity measures using state of the art criterion established in Performance-Based Earthquake Engineering and techniques from Information Theory. The efficiency and the sufficiency of several rock intensity measure are assessed. It was found that spectral accelerations at low periods usually are adequate for predicting amplifications across the period range. This supports the current practice of using Peak Ground Acceleration in Ground Motion Models. However, for extremely soft sites, which demonstrate nonlinear effects well into the period range, it is better practice to ensure that amplification and spectral acceleration share the same oscillator period. Lastly, for predicting the peak shear strain (an important parameter that controls nonlinearity of site response), Peak Ground Velocity is generally adequate, and this conclusion is in line with the definition of proxy shear strain.

42 ENGINEERING↗

Mass transport in a highly immiscible alloy on extended shear deformation

Forced mixing to a single phase or supersaturated solid solution (SSS) and its prerequisite microstructure evolution in immiscible systems has been a focus of research for both fundamental science and a variety of applications. Controlling the formation of SSS by shear deformation assisted processing could enable a material design beyond conventional equilibrium microstructure in immiscible systems. Here, a highly immiscible (mixing enthalpy of ~ 20 kJ·mol -1 ) Cu-50 at. % Cr binary alloy was employed to investigate the microstructure evolution and localized tendencies of SSS during severe shear deformation. The present results demonstrate distinctive microstructure refinement process in each phase and how they lead to localized SSS as a function of shear strain. Preferential dynamic recrystallization occurs in the softer Cu phase due to strain localization, leading to substantial grain refinement. The refinement of Cr phase in the top-most layer, however, is enabled by the progressive evolution of grain lamination, splitting, and fragmentation as a function of shear strain. The eventual SSS is found to be strongly dependent on the local environments that affect the dislocation activity, including the level of microstructure refinement, the interfacial orientation relationship, hardness difference, and supposed stability of oxidation. Ab initio simulations comfirm that it is more favorable to oxidize Cr than Cu at incoherent Cu/Cr interfaces which then limits the mass transport on an incoherent boundary. Our results shed light on the underpinning mechanism for non-equilibrium mass transport in immiscible systems upon severe deformation that can be applicable to a variety of processing techniques aimed at producing immiscible alloys with superior mechanical properties.

36 MATERIALS SCIENCE↗

Microstructure, Texture, and Strength Development during High-Pressure Torsion of CrMnFeCoNi High-Entropy Alloy

The equiatomic face-centered cubic high-entropy alloy CrMnFeCoNi was severely deformed at room and liquid nitrogen temperature by high-pressure torsion up to shear strains of about 170. Its microstructure was analyzed by X-ray line profile analysis and transmission electron microscopy and its texture by X-ray microdiffraction. Microhardness measurements, after severe plastic deformation, were done at room temperature. It is shown that at a shear strain of about 20, a steady state grain size of 24 nm, and a dislocation density of the order of 1016 m−2 is reached. The dislocations are mainly screw-type with low dipole character. Mechanical twinning at room temperature is replaced by a martensitic phase transformation at 77 K. The texture developed at room temperature is typical for sheared face-centered cubic nanocrystalline metals, but it is extremely weak and becomes almost random after high-pressure torsion at 77 K. The strength of the nanocrystalline material produced by high-pressure torsion at 77 K is lower than that produced at room temperature. The results are discussed in terms of different mechanisms of deformation, including dislocation generation and propagation, twinning, grain boundary sliding, and phase transformation.

36 MATERIALS SCIENCE↗

Strong increase in ultrasound attenuation below T c in Sr 2 RuO 4 : Possible evidence for domains

Recent experiments suggest that Sr 2 RuO 4 has a two-component superconducting order parameter (OP). A two-component OP has multiple degrees of freedom in the superconducting state that can result in low-energy collective modes or the formation of domain walls—a possibility that would explain a number of experimental observations including the smallness of the signature of time reversal symmetry breaking at T c and telegraph noise in critical current experiments. We use resonant ultrasound spectroscopy to perform ultrasound attenuation measurements across the superconducting T c of Sr 2 RuO 4 . We find that compressional sound attenuation increases by a factor of 7 immediately below T c , in sharp contrast with what is found in both conventional (s-wave) and high-T c (d-wave) superconductors. Our observations are most consistent with the presence of domain walls that separate different configurations of the superconducting OP. Here, the fact that we only observe an increase in sound attenuation for compressional strains, and not for shear strains, suggests an inhomogeneous superconducting state formed of two distinct, accidentally degenerate superconducting OPs that are not related to each other by symmetry. Whatever the mechanism, a factor of 7 increase in sound attenuation is a singular characteristic that must be reconciled with any potential theory of superconductivity in Sr 2 RuO 4 .

36 MATERIALS SCIENCE↗

Analyzing the Validity of Brazilian Testing Using Digital Image Correlation and Numerical Simulation Techniques

Characterizing the mechanical behavior of rocks plays a crucial role to optimize the fracturing process in unconventional reservoirs. However, due to the intrinsic anisotropy and heterogeneity in unconventional resources, fracture process prediction remains the most significant challenge for sustainable and economic hydrocarbon production. During the deformation tracking under compression, deploying conventional methods (strain gauge, extensometer, etc.) is insufficient to measure the deformation since the physical attachment of the device is restricted to the size of the sample, monitoring limited point-wise deformation, producing difficulties in data retrieval, and a tendency to lose track in failure points, etc. Where conventional methods are limited, the application of digital image correlation (DIC) provides detailed and additional information of strain evolution and fracture patterns under loading. DIC is an image-based optical method that records an object with a camera and monitors the random contrast speckle pattern painted on the facing surface of the specimen. To overcome the existing limitations, this paper presents numerical modeling of Brazilian disc tests under quasi-static conditions to understand the full-field deformation behaviors and finally, it is validated by DIC. As the direct tensile test has limitations in sample preparation and test execution, the Brazilian testing principle is commonly used to evaluate indirectly the tensile strength of rocks. The two-dimensional numerical model was built to predict the stress distribution and full-field deformation on Brazilian disc under compression based on the assumptions of a homogenous, isotropic and linear elastic material. The uniaxial compression test was conducted using the DIC technique to determine the elastic properties of Spider Berea sandstone, which were used as inputs for the simulation model. The model was verified by the analytical solution and compared with the digital image correlation. The numerical simulation results showed that the solutions matched reasonably with the analytical solutions where the maximum deviation of stress distribution was obtained as 14.59%. The strain evolution (normal and shear strains) and displacements along the central horizontal and vertical planes were investigated in three distinguishable percentages of peak loads (20%, 40%, and 90%) to understand the deformation behaviors in rock. The simulation results demonstrated that the strain evolution contours consistently matched with DIC generated contours with a reasonable agreement. The changes in displacement along the central horizontal and vertical planes showed that numerical simulation and DIC generated experimental results were repeatable and matched closely. In terms of validation, Brazilian testing to measure the indirect tensile strength of rocks is still an issue of debate. The numerical model of fracture propagation supported by digital image correlation from this study can be used to explain the fracturing process in the homogeneous material and can be extended to non-homogeneous cases by incorporating heterogeneity, which is essential for rock mechanics field applications.

58 GEOSCIENCES↗

A method for characterization of multiple dynamic constitutive parameters of FRCs

We propose a method to measure multiple dynamic material constitutive parameters of unidirectional fiber reinforced composites (FRCs) in a single experiment. Dynamic short-beam shear (DSBS) experiments were performed on a modified Kolsky compression bar, with integration of high-speed imaging and digital image correlation (DIC). The unidirectional FRCs investigated were S-2 glass fiber reinforced matrix of TGDDM-Jeffamine® D230 with monoamine functionalized partially reacted substructures (mPRS) and commercially available SC-15. Analytical solutions of normal and shear strains of a composite beam were derived based on Timoshenko beam theory, assuming material to be transversely isotropic and have different moduli in tension and compression in each principle material orientation. Tensile and compressive moduli were inversely computed through monitoring normal strain slope when specimen was constantly loaded at a speed of ~7.3 m/s within a small deflection. Non-linear shear stress-strain behavior of the composite was described via Ramberg–Osgood equation. Finite element (FE) analysis was conducted in ABAQUS, simultaneously defining via user subroutine UMAT the transverse isotropy of material, bi-modulus constitutive model, and non-linear shear stress-strain relation. The method proposed in this work was validated by comparing strain distributions computed by FE model and DIC measurements. Comparing with traditional dynamic tensile, compressive, and shear experiments on FRCs, this method significantly simplifies the specimen preparation and design of complicated gripping fixtures for multiple experiments. Here, systematic errors resulting from variations of specimen geometry and dimension, loading direction, and instrumentation are reduced, thereby providing compatible data for numerical studies on impact behavior of composites.

36 MATERIALS SCIENCE↗

A finite element formulation for deformation twinning induced strain localization in polycrystal magnesium alloys

Deformation twinning induces shear strain localization in hexagonal close-packed crystals and is critical for the material’s ductility and failure. Cracks often occur at twin-twin or twin-grain boundary intersections and propagate along twin bands. However, most crystal plasticity models for deformation twinning are based on a “pseudo-slip” approach and do not capture the localized deformation associated with the formation of each discrete twin band. The few exceptions are discrete twin models that involve very complex numerical algorithms and are often compromised in accuracy due to the numerical convergence. These factors make the discrete twin models hard to adopt. This paper proposes a modification to the conventional finite element weak form, to fully incorporate a twin-induced heterogeneous deformation that does not depend on the “pseudo-slip” assumption. The model starts by splitting the deformation gradient into elastic-slip-twinning components. The twin-induced deformation gradient component is computed separately by solving a microstructural evolution problem and then implemented into finite element weak form by constructing a global “twin-force” vector. The constitutive update (e.g., in the user-defined material subroutine, or UMAT, for ABAQUS) therefore avoids dealing with the twinning and recovers to the form of a regular slip-based crystal plasticity model. The results presented here indicate that the twin-induced strain localization and the associated stress-reversal phenomena near the twin band were naturally captured in the model, which was validated against an in-situ synchrotron X-ray micro-diffraction experiment.

36 MATERIALS SCIENCE↗

Bending and shear improvements in 3D-printed core sandwich composites through modification of resin uptake in the skin/core interphase region

Fused deposition modeling (FDM) printed polymers are rarely used as a structural material due to anisotropic and low mechanical properties compared with conventional composites. In recent years, greater need has been expressed for recycling of materials, such as recyclable FDM, at the end of service life to reduce environmental pollution and manufacture cost. However, how the amount of resin uptake in the skin and skin/core interphase affects the bending and shear performance of the sandwich composites when replacing the low strength and ductile core (conventional core) with a high strength and brittle core (FDM printed PLA (polylactic acid) core) still remains unclear. A new manufacturing routine is needed to improve the incorporation of FDM printed polymers in composite structures. In this work, FDM printed PLA was used as core material and sandwiched between two unidirectional glass fiber reinforced polymer (GFRP) skins to form a sandwich composite by compression-molding (CM) process, which provides a good manufacturing strategy for skin/core interphase modification. The significance of the CM process is proved by investigating the effect of resin uptake on bending and in-plane/out-of-plane shear performances. Current first order shear deformation (FSDT) theory lacks a direct connection between the in-plane shear stress and out-of-shear stress in the core region of sandwich composites. With the help of DIC, a connection between the in-plane shear and the out-of-plane shear strain was built and in-plane shear properties can acquire through out-of-plane shear properties, hence reducing the redundancy of sample preparation or the need for simulation. A significant improvement was found compared with the optimized resin uptake (Optimized resin uptake range: 20.43%–22.86 wt%) 3D-printed PLA core sandwich composite and lowest performance sandwich composite (Improvement: in-plane shear strength (~34%)/modulus (~29%), out-of-plane shear strength (~25%)/modulus (~31%), specific peak bending load (~19%)). Finally, compared with balsa core sandwich composites, the 3D-printed cores are suitable for use in composite sandwich structures in many applications with a satisfactory strength-to-weight ratio.

3D-printed core sandwich composite↗

Estimating epistemic uncertainty in soil parameters for nonlinear site response analyses: Introducing the Latin Hypercube Sampling technique

This study quantifies the effects of epistemic uncertainty in soil parameters on nonlinear (NL) site response analysis (SRA) results, validated against the data recorded at a well-instrumented geotechnical downhole array located in Japan. To this end, a one-dimensional soil column model of the Service Hall Array (SHA) near the Kashiwazaki-Kariwa Nuclear Power Plant (KKNPP) is developed using the finite element (FE) program LS-DYNA. The dynamic stress–strain relationship is characterized by a modified two-stage hyperbolic (MTH) NL backbone curve formulation capable of capturing soil behavior at both small- and large-shear strains. The model is then validated against the ground motion recordings to capture the model bias. The uncertainties associated with the shear-wave velocity profile (a small-strain soil property) and soil shear strength (a large-strain soil property) are incorporated in NL SRA to quantify their separate and joint randomization effects on the results. This study proposes using the Latin Hypercube Sampling (LHS) method as an efficient alternative to commonly used methods, such as Standard Monte Carlo (SMC), to account for uncertainty propagation in such reliability analysis. Both low-intensity and design-level records from the recordings at the SHA are applied to study the contribution of the small- and large-strain NL dynamic soil properties. Results from 46,200 NL FE analyses (23,100 per input motion) are presented. Measured and predicted site response, using recorded ground motions at this downhole array, is compared to assess the significance of soil parameter uncertainty on the observed ground motion dispersions. It is demonstrated that increasing the number of soft realizations and implementing higher level earthquake intensity lead to higher ground motion dispersion. Unlike past studies in randomization of Vs profiles with the SMC method, the LHS method is shown to have no significant effect on the predicted median surface response spectra and amplification factors (AFs) for this case study.

Engineering↗

Microstructural evolution and intermetallic formation in Zn-3Mg (wt%) powder mixture processed by high-pressure torsion

Severe plastic deformation (SPD) techniques have been used extensively over the past 40 years for producing strong metals and alloys. High-pressure torsion (HPT) is one of the most promising SPD techniques for achieving high strength through nanoscale grain refinement and phase transformation. In this research, a mixture of pure zinc (Zn) and magnesium (Mg) powders, Zn-3Mg (wt%), was HPT-processed under a pressure of 6 GPa for 1, 5, 10, 20, and 30 turns at room temperature to achieve a high strength biodegradable material. In order to understand the effects of pre-consolidation on the resulting microstructure and hardness, HPT processing was performed on loose powders placed in the die and also on a pre-compacted powder mixture and the resulting HPT disks were characterized by X-ray diffraction, scanning electron microscopy, atom probe tomography, and Vickers microhardness. In both cases, the HPT disk microstructures contained nanoscale grains, and stable and metastable strain-induced intermetallics, but an unusual softening appeared at large shear strains. Grain size, grain morphology, and the formation of different intermetallics were analyzed to explain the unusual hardness distribution, and it was found that an inverse Hall-Petch relationship between hardness and grain size exists. It is suggested that thermally-activated phenomena such as grain boundary sliding contributed to the strain-induced softening of this nano-structured biomaterial due to its low melting point. Finally, the current results are compared with those for HPT-processed cast alloys and hybrids of the same composition.

36 MATERIALS SCIENCE↗

Evolution of Elastic and Mechanical Properties During Fault Shear: The Roles of Clay Content, Fabric Development, and Porosity

Phyllosilicates weaken faults due to the formation of shear fabrics. Although the impacts of clay abundance and fabric on frictional strength, sliding stability, and porosity of faults are well studied, their influence on elastic properties is less known, though they are key factors for fault stiffness. We document the role that fabric and consolidation play in elastic properties and show that smectite content is the most important factor determining whether fabric or porosity controls the elastic response of faults. We conducted a suite of shear experiments on synthetic smectite-quartz fault gouges (10–100 wt% smectite) and sediment incoming to the Sumatra subduction zone. We monitored Vp, Vs, friction, porosity, shear and bulk moduli. We find that mechanical and elastic properties for gouges with abundant smectite are almost entirely controlled by fabric formation (decreasing mechanical and elastic properties with shear). Though fabrics control the elastic response of smectite-poor gouges over intermediate shear strains, porosity is the primary control throughout the majority of shearing. Elastic properties vary systematically with smectite content: High smectite gouges have values of Vp ~ 1,300–1,800 m/s, Vs ~ 900–1,100 m/s, K ~ 1–4 GPa, and G ~ 1–2 GPa, and low smectite gouges have values of Vp ~ 2,300–2,500 m/s, Vs ~ 1,200–1,300 m/s, K ~ 5–8 GPa, and G ~ 2.5–3 GPa. We find that, even in smectite-poor gouges, shear fabric also affects stiffness and elastic moduli, implying that while smectite abundance plays a clear role in controlling gouge properties, other fine-grained and platy clay minerals may produce similar behavior through their control on the development of fabrics and thin shear surfaces.

58 GEOSCIENCES↗

Shear Softening in a Metallic Glass: First-Principles Local-Stress Analysis

Metallic glasses deform elastically under stress. However, the atomic-level origin of elastic properties of metallic glasses remain unclear. In this Letter using ab initio molecular dynamics simulations of the Cu 5 0Zr 50 metallic glass under shear strain, we show that the heterogeneous stress relaxation results in the increased charge transfer from Zr to Cu atoms, enhancing the softening of the shear modulus. Changes in compositional short-range order and atomic position shifts due to the nonaffine deformation are discussed. It is shown that the Zr subsystem exhibits a stiff behavior, whereas the displacements of Cu atoms from their initial positions, induced by the strain, provide the stress drop and softening.

36 MATERIALS SCIENCE↗

Shear-Induced Conformations of Salt-Free Polyelectrolytes in Semidilute Solutions

Here, we present coarse-grained molecular dynamics simulations of salt-free polyelectrolyte chains in semidilute solutions under simple shear flow, with full hydrodynamic interactions and explicit dipolar solvent. At equilibrium, chain orientation statistics follow a pseudo-Voigt distribution, and the structural correlation length and chain end-to-end vector autocorrelation function exhibit scaling behavior consistent with theoretical predictions for polyelectrolytes. Under shear, chains transition from coiled to stretched states and the end-to-end vector autocorrelation function reveals oscillatory dynamics at high Weissenberg numbers. Analysis of the gyration tensor and shear strain distributions identified three distinct chain populations with directional alignment along and against the shear gradient. Compared with their neutral polymer counterparts, polyelectrolytes exhibit stronger shear thinning and enhanced chain alignment under the same Weissenberg number, which is attributed to electrostatic interactions and shear-induced counterion release. These findings provide molecular insight into the distinct flow response of charged polymers and have implications for tailoring the rheological properties of polyelectrolyte-based materials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Programmable nanowrinkle-induced room-temperature exciton localization in monolayer WSe2

Abstract Localized states in two-dimensional (2D) transition metal dichalcogenides (TMDCs) have been the subject of intense study, driven by potential applications in quantum information science. Despite the rapidly growing knowledge surrounding these emitters, their microscopic nature is still not fully understood, limiting their production and application. Motivated by this challenge, and by recent theoretical and experimental evidence showing that nanowrinkles generate strain-localized room-temperature emitters, we demonstrate a method to intentionally induce wrinkles with collections of stressors, showing that long-range wrinkle direction and position are controllable with patterned array design. Nano-photoluminescence (nano-PL) imaging combined with detailed strain modeling based on measured wrinkle topography establishes a correlation between wrinkle properties, particularly shear strain, and localized exciton emission. Beyond the array-induced wrinkles, nano-PL spatial maps further reveal that the strain environment around individual stressors is heterogeneous due to the presence of fine wrinkles that are less deterministic. At cryogenic temperatures, antibunched emission is observed, confirming that the nanocone-induced strain is sufficiently large for the formation of quantum emitters. At 300 K, detailed nanoscale hyperspectral images uncover a wide range of low-energy emission peaks originating from the fine wrinkles, and show that the states can be tightly confined to regions <10 nm, even in ambient conditions. These results establish a promising potential route towards realizing room temperature quantum emission in 2D TMDC systems.

36 MATERIALS SCIENCE↗

Phase transitions in high-purity zirconium under dynamic compression

In this work, we present results from ramp compression experiments on high-purity Zr that show the α → ω, ω → β, as well as reverse β → ω phase transitions. Simulations with a multi-phase equation of state and phenomenological kinetic model match the experimental wave profiles well. While the dynamic α ! ! transition occurs ~ 9 GPa above the equilibrium phase boundary, the ω → β transition occurs within 0.9 GPa of equilibrium. We estimate that the dynamic compression path intersects the equilibrium ω – β line at P = 29:2 GPa, and T = 490 K. The thermodynamic path in the interior of the sample lies ~ 100 K above the isentrope at the point of the ω → β transition. Approximately half of this dissipative temperature rise is due to plastic work, and half is due to the non-equilibrium α → ω transition. The inferred rate of the α → ω transition is several orders of magnitude higher than that measured in dynamic diamond anvil cell (DDAC) experiments in an overlapping pressure range. We discuss a model for the influence of shear stress on the nucleation rate. We find that the shear stress s ji has the same effect on the nucleation rate as a pressure increase δP = cϵ ij s ji /(ΔV/V); where c is a geometric constant ~ 1 and, ϵ ij are the transformation shear strains. The small fractional volume change ΔV/V ≈ 0:1 at the α → ω transition amplifies the effect of shear stress, and we estimate that for this case δP is in the range of several GPa. Correcting our transition rate to a hydrostatic rate brings it approximately into line with the DDAC results, suggesting that shear stress plays a significant role in the transformation rate.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗