Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “crystal structure prediction”

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 631 records · Page 35

Improving nanoparticle superlattice stability with deformable polymer gels

The self-assembly of colloidal nanoparticles into ordered superlattices typically uses dynamic interactions to govern particle crystallization, as these non-permanent bonds prevent the formation of kinetically trapped, disordered aggregates. However, while the use of reversible bonding is critical in the formation of highly ordered particle arrangements, dynamic interactions also inherently make the structures more prone to disassembly or disruption when subjected to different environmental stimuli. Thus, there is typically a trade-off between the ability to initially form an ordered colloidal material and the ability of that material to retain its order under different conditions. Here, we present a method for embedding colloidal nanoparticle superlattices into a polymer gel matrix. This encapsulation strategy physically prevents the nanoparticles from dissociating upon heating, drying, or the introduction of chemicals that would normally disrupt the lattice. However, the use of a gel as the embedding medium still permits further modification of the colloidal nanoparticle lattice by introducing stimuli that deform the gel network (as this deformation in turn alters the nanoparticle lattice structure in a predictable manner). Moreover, encapsulation of the lattice within a gel permits further stabilization into fully solid materials by removing the solvent from the gel or by replacing the solvent with a liquid monomer that can be photopolymerized. This embedding method therefore makes it possible to incorporate ordered colloidal arrays into a polymer matrix as either dynamic or static structures, expanding their potential for use in responsive materials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A Trichomonas vaginalis C2-XYPPX-repeat protein with a structured C2 domain displaying dampened flexibility upon binding calcium

C2 domains are ubiquitous membrane-binding modules of ∼130 residues in eukaryotes that are often associated with proteins involved in membrane trafficking and lipid modification. The genome of Trichomonas vaginalis, the most common, non-viral, sexually transmitted human pathogen, encodes eight genes that contain a N-terminal C2 module linked to a XYPPX-repeat domain of more than four XYPPX repeats (C2-XYPPX). While the function of the XYPPX-repeat domain remains unknown, its multiple association with C2 domains in T. vaginalis suggests it is important. Here, the C2 domain from one of these C2-XYPPX-repeat proteins, Tv-C2-1, was structurally and physically characterized using X-ray crystallography and NMR spectroscopy. The crystal structure for Tv-C2-1 shows that this domain shares a fold common to all C2 domains, a compact Greek-key motif composed of eight anti-parallel β-strands in the type-2 topology. An NMR chemical shift perturbation study with Ca 2+ showed that Tv-C2-1 bound two Ca 2+ atoms primarily via two loops (loop-1 and loop-3) on the predicted calcium binding face of the protein with K d s of 58.0 ± 0.1 μM and 232 ± 6 μM. Estimations of the overall rotational correlation time, τ c , in the apo (11.1 ns) and Ca 2+ -bound (9.2 ns) state suggests the protein becomes more compact upon Ca 2+ binding, consistent with a decrease in dynamics in loop-3 and marginally in loop-1 suggested by amide 15 N heteronuclear steady-state { 1 H}- 15 N NOEs. Showing Tv-C2-1 binds calcium and adopts a compact Greek-key motif structure, two primary features of C2 domains, suggests understanding the function of the XYPPX-repeat domain may be warranted.

NMR spectroscopy↗

Discovery and Characterization of a Metastable Cubic Interstitial Nickel–Carbon System with an Expanded Lattice

Metastable, i.e., kinetically favored but thermodynamically not stable, interstitial solid solutions of carbon in iron are well-understood. Carbon can occupy the interstitial atoms of the host metal, altering its properties. Alloying of the host metal results in the stabilization of the FeC x phases, widening its application. Pure nickel finds niche applications, mainly focusing on catalysis, while nickel alloys are widely applied, e.g., in gas turbines, reactors, and seawater piping. Nickel carbide (Ni 3 C) is the well-known stable Ni–C system displaying a rhombohedral (R3̅c) crystal structure. Some reports describe an elusive cubic Ni–C system, observed during certain catalytic reactions occurring on nickel and formed by the occupation of the interstitials of the metal with carbon: to date, the stabilization and characterization of this phase have not been accomplished. Hereby, we report on the synthesis of a cubic metastable NiC x phase using chemical vapor deposition of methane on supported nickel nanoparticles. The structure was predicted by DFT/ReaxFF, synthesized and monitored with in situ time-resolved synchrotron XRD, and experimentally confirmed by Rietveld refinement and (S)TEM-EELS under ambient conditions. The results show an Fm3̅m phase with a lattice parameter of a = 3.749 ± 0.037 Å at room temperature, with the highest ever reported atomic percentage of carbon occupying the octahedral interstices of 23.1%, resulting in a NiC 0.3 phase. The degree of occupation of the interstitial voids by carbon can be controlled, enabling the tuning of the host metal’s d-spacing and composition, highlighting the applicability of this synthesis route for catalytic nanoparticle preparation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Sas20 is a highly flexible starch-binding protein in the Ruminococcus bromii cell-surface amylosome

Ruminococcus bromii is a keystone species in the human gut that has the rare ability to degrade dietary resistant starch (RS). This bacterium secretes a suite of starch-active proteins that work together within larger complexes called amylosomes that allow R. bromii to bind and degrade RS. Starch adherence system protein 20 (Sas20) is one of the more abundant proteins assembled within amylosomes, but little could be predicted about its molecular features based on amino acid sequence. Here, we performed a structure–function analysis of Sas20 and determined that it features two discrete starch-binding domains separated by a flexible linker. We show that Sas20 domain 1 contains an N-terminal β-sandwich followed by a cluster of α-helices, and the nonreducing end of maltooligosaccharides can be captured between these structural features. Furthermore, the crystal structure of a close homolog of Sas20 domain 2 revealed a unique bilobed starch-binding groove that targets the helical α1,4-linked glycan chains found in amorphous regions of amylopectin and crystalline regions of amylose. Affinity PAGE and isothermal titration calorimetry demonstrated that both domains bind maltoheptaose and soluble starch with relatively high affinity (K d ≤ 20 μM) but exhibit limited or no binding to cyclodextrins. Finally, small-angle X-ray scattering analysis of the individual and combined domains support that these structures are highly flexible, which may allow the protein to adopt conformations that enhance its starch-targeting efficiency. Taken together, we conclude that Sas20 binds distinct features within the starch granule, facilitating the ability of R. bromii to hydrolyze dietary RS.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Structures of two novel crystal forms of Aspergillus oryzae alpha amylase (taka-amylase)

The structures of Aspergillus oryzae α-amylase were determined in a tetragonal crystal, having one molecule as asymmetric unit, and a monoclinic crystal with two molecules as asymmetric unit. Both crystal forms were obtained from trace contaminants of an old commercial lipase preparation. Structures were determined and refined to 1.65 Å and 1.43 Å resolution respectively. The latter crystal has a non-crystallographic (NCS) twofold axis within the asymmetric unit. Glycosylation at Asn197 is evident, and in the tetragonal crystal can be seen to include three, partially disordered sugar residues following the initial N-acetyl glucosamine (NAG). Superposition of the tetragonal crystal model on the α-amylases from Bacillus subtilis (PDB:1BAG), pig pancreas (PDB:3L2L), and barley (PDB:1AMY), show a high degree of coincidence, particularly for the (β/α)8-barrel domains, and especially within the active site. Using this structural agreement between amylases, we extrapolated the binding model of a six residue, limit dextrin found in pig pancreas α-amylase to the A. oryzae enzyme model, which predicts substrate interacting amino acid residues.

59 BASIC BIOLOGICAL SCIENCES↗

First-principles investigation of phase stability in layered Na x ⁢CrO 2

Layered oxide intercalation compounds continue to attract interest as electrode materials for Na-ion batteries. However, many of these materials undergo complex phase transitions during cycling that influence battery performance but are still not completely understood. We have conducted a first-principles study of layered Na x ⁢CrO 2 (0 ≤ x ≤ 1) to assess phase stability between various Na-vacancy ordered phases in the O3 and P3 host structures. We predict that many of the low-energy phases belong to families of Na orderings that follow specific patterns. At high x, we identify families of vacancy row orderings in O3, which may also couple to magnetic orderings of the Cr spins. We predict similar orderings at intermediate x in P3 that contain antiphase boundaries between regions of the x = 1/2 ordering. In both cases, the average spacing between rows/boundaries is set by the overall composition. At x = 0, we find a strong preference for charge disproportionation and migration of Cr to tetrahedral sites in the intercalation layers. Finally, we obtain generally good agreement with experimental observations and rationalize key discrepancies.

25 ENERGY STORAGE↗

Controlling magnetic order, magnetic anisotropy, and band topology in the semimetals Sr(Mn 0.9 Cu 0.1 )Sb 2 and Sr(Mn 0.9 Zn 0.1 )Sb 2

Neutron diffraction and magnetic susceptibility studies show that orthorhombic single-crystals of topological semimetals Sr(Mn 0.9 Cu 0.1 ) Sb 2 and Sr(Mn 0.9 Zn 0.1 )Sb 2 undergo three-dimensional C-type antiferromagnetic (AFM) ordering of the Mn 2+ moments at T N = 200 ± 10 and 210 ± 12 K, respectively, significantly lower than that of the parent SrMnSb 2 with T N = 297 ± 3 K. Magnetization versus applied magnetic field (perpendicular to MnSb planes) below T N exhibits slightly modified de Haas van Alphen oscillations for the Zn-doped crystal as compared to that of the parent compound. By contrast, the Cu-doped system does not show de Haas van Alphen magnetic oscillations, suggesting that either Cu substitution for Mn changes the electronic structure of the parent compound substantially, or that the Cu sites are strong scatterers of carriers that significantly shorten their mean free path thus diminishing the oscillations. Density functional theory (DFT) calculations including spin-orbit coupling predict the C-type AFM state for the parent, Cu-, and Zn-doped systems and identify the a -axis (i.e., perpendicular to the Mn layer) as the easy magnetization direction in the parent and 12.5% of Cu or Zn substitutions. In contrast, 25% of Cu content changes the easy magnetization to the b-axis (i.e., within the Mn layer). Here, we find that the incorporation of Cu and Zn in SrMnSb 2 tunes electronic bands near the Fermi level resulting in different band topology and semimetallicity. The parent and Zn-doped systems have coexistence of electron and hole pockets with opened Dirac cone around the Y-point whereas the Cu-doped system has dominant hole pockets around the Fermi level with a distorted Dirac cone. The tunable electronic structure may point out possibilities of rationalizing the experimentally observed de Haas van Alphen magnetic oscillations.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Occupational disorder as the origin of flattening of the acoustic phonon branches in the clathrate Ba 8 ⁢Ga 16 ⁢Ge 30

In the search for high-performance thermoelectrics, materials such as clathrates have drawn attention due to having both glasslike low phonon thermal conductivity and crystal-like high electrical conductivity. Ba 8 ⁢Ga 16 ⁢Ge 30 (BGG) has a loosely bound guest Ba atom trapped inside rigid Ga-Ge cage structures. Avoided crossings between acoustic phonons and the flat guest atom branches have been proposed to be the source of the low lattice thermal conductivity of BGG. Ga-Ge site disorder with Ga and Ge exchanging places in different unit cells has also been reported. We used time-of-flight neutron scattering to measure the complete phonon spectrum in a large single crystal of BGG and compared these results with predictions of density functional theory to elucidate the effect of the disorder on heat-carrying phonons. Experimental results agreed much better with the calculation assuming the disorder than with the calculation assuming the ordered configuration. Although the atomic masses of Ga and Ge are nearly identical, we found that disorder strongly reduces phonon group velocities, which significantly reduces thermal conductivity. Finally, our work points to a path towards optimizing thermoelectrics.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

An operator-based approach to topological photonics

Abstract Recently, the study of topological structures in photonics has garnered significant interest, as these systems can realize robust, nonreciprocal chiral edge states and cavity-like confined states that have applications in both linear and nonlinear devices. However, current band theoretic approaches to understanding topology in photonic systems yield fundamental limitations on the classes of structures that can be studied. Here, we develop a theoretical framework for assessing a photonic structure’s topology directly from its effective Hamiltonian and position operators, as expressed in real space, and without the need to calculate the system’s Bloch eigenstates or band structure. Using this framework, we show that nontrivial topology, and associated boundary-localized chiral resonances, can manifest in photonic crystals with broken time-reversal symmetry that lack a complete band gap, a result that may have implications for new topological laser designs. Finally, we use our operator-based framework to develop a novel class of invariants for topology stemming from a system’s crystalline symmetries, which allows for the prediction of robust localized states for creating waveguides and cavities.

97 MATHEMATICS AND COMPUTING↗

A Platform-Independent X-ray Diffraction Diagnostic for Phase Transition Kinetics in Traditional and Synthetic Microstructure Materials (LDRD Project 213088 Final Report)

Pulsed-power generators using the magnetic loading technique are able to produce well-controlled continuous ramp compression of condensed matter for high-pressure equation-of-state studies. X-ray diffraction (XRD) data from dynamically compressed samples provide direct measurements of the elastic compression of the crystal lattice, onset of plastic flow, strength-strain rate dependence, structural phase transitions, and density of crystal defects such as dislocations. Here, we present a cost effective, compact X-ray source for XRD measurements on pulsed-power-driven ramp-loaded samples. This combination of magnetically-driven ramp compression of materials with single, short-pulse XRD diagnostic will be a powerful capability for the dynamic materials community. The success in fielding this new XRD diagnostic dramatically improves our predictive capability and understanding of rate-dependent behavior at or near phase transition. As Sandia plans the next-generation pulse-power driver platform, a key element needed to deliver new state-of-the-art experiments will be having the necessary diagnostic tools to probe new regimes and phenomena. These diagnostics need to be as versatile, compact, and portable as they are powerful. The development of a platform-independent XRD diagnostic gives Sandia researchers a new window to study the microstructure and phase dynamics of materials under load. This project has paved the way for phase transition research in a variety of materials with mission interest.

36 MATERIALS SCIENCE↗

Data Driven Approach to Dislocation-Based Plasticity Models of Face-Centered Cubic Metals

Dislocation dynamics controls plastic deformation, mechanical strength, and failure of crystalline materials. It also governs fatigue resistance under cyclic loading, creep resistance at elevated-temperature, and radiation resistance for reactor applications. There is a compelling need for understanding fundamental dislocation mechanisms for deformation because virtually all structural metals used in energy systems are fabricated to desired forms and shapes by deformation processes. To date, the most outstanding problem in a physics-based multiscale model of crystal plasticity is the lack of quantitative connections between continuum plasticity (CP) models with the lower scale dislocation models. As a result, existing CP models used in engineering applications are still phenomenological, while evidence continues to mount that they can make inaccurate predictions under realistically complex scenarios. This project takes advantage of the recent advances in high-performance discrete dislocation dynamics (DDD) simulations and data science approaches to establish the first fully connected multiscale plasticity model for pure face-centered cubic (FCC) single crystals.

36 MATERIALS SCIENCE↗

Decoding Zeolite Crystallization and Stage III in Nuclear Waste Glasses by Coupled Modeling and Experiments

Under specific conditions of pH and temperature, nuclear waste immobilization borosilicate glasses may exhibit a sudden acceleration in their corrosion kinetics (stage III)—a behavior that has been associated with the formation of zeolite crystals. Such accelerated dissolution may compromise the integrity of nuclear wasteforms placed in geological depositories. However, thus far, none of the available models is able to predict the thermodynamic propensity and kinetics of zeolite precipitation as a function of the solution conditions due to (i) a lack of fundamental knowledge regarding the nucleation & growth mechanisms of zeolitic phases, (ii) uncertainty regarding the compositions (types) of zeolites that may form and the rate-limiting step in their precipitation as a function of the solution conditions, and (iii) the complexities that arise due to the vast parametric space (i.e., solution chemistry, temperature, number of secondary phases, etc.) that encompass these systems under conditions of environmental exposure. To resolve these challenges, this project aimed to unambiguously identify the thermodynamic propensity for zeolite precipitation and the kinetics thereof as a function of the solution conditions (composition, pH, and temperature). To achieve this goal: 1) We identified the solution conditions and zeolite phases relevant to nuclear glass dissolution. 2) We performed a series of ab initio molecular dynamics (AIMD) simulations to compute the thermodynamic properties of a group of characteristic zeolites that features a large range of compositions, various hydration levels, a wide range of framework structures, and partial atomic site occupancies. 3) We released a first-of-a-kind self-consistent thermodynamic database that can be used to assess the kinetics and the stability fields of zeolitic phases within a Gibbs energy minimization (GEM) framework. 4) We developed a robust geochemical modeling method allowing us to predict the stability of secondary phases (including zeolites, calcium–silicate–hydrate gels, and clays) upon the dissolution of nuclear waste immobilization glasses. 5) We introduced a model that predicts the dissolution kinetics of a series of borosilicate nuclear waste immobilization glasses in terms of the topology of their atomic network. 6) We investigated the roles of the solution composition on the crystallization kinetics of phillipsite zeolites and tobermorite silicate hydrates. Via PNNL’s collaboration and engagement, this project directly supports DOE’s nuclear waste immobilization activities by offering a technical, science-based foundation that will (i) facilitate predictions of the long-term corrosion rates and extents of existing nuclear waste immobilization glasses to help ensure safe and successful vitrification operations, (ii) inform the development of advanced glass formulations with enhanced durability, and, (iii) enable cost-savings that result from making more decisive and hence less conservative predictions while offering higher levels of nuclear waste embedment in smaller, more compact glass volumes.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Liquid-liquid phase separation in supercooled water from ultrafast heating of low-density amorphous ice

Recent experiments continue to find evidence for a liquid-liquid phase transition (LLPT) in supercooled water, which would unify our understanding of the anomalous properties of liquid water and amorphous ice. These experiments are challenging because the proposed LLPT occurs under extreme metastable conditions where the liquid freezes to a crystal on a very short time scale. Here, we analyze models for the LLPT to show that coexistence of distinct high-density and low-density liquid phases may be observed by subjecting low-density amorphous (LDA) ice to ultrafast heating. We then describe experiments in which we heat LDA ice to near the predicted critical point of the LLPT by an ultrafast infrared laser pulse, following which we measure the structure factor using femtosecond x-ray laser pulses. Consistent with our predictions, we observe a LLPT occurring on a time scale <100 ns and widely separated from ice formation, which begins at times >1 μs.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

In Situ Insights into Cathode Calcination for Predictive Synthesis: Kinetic Crystallization of LiNiO 2 from Hydroxides

Abstract Calcination is a solid‐state synthesis process widely deployed in battery cathode manufacturing. However, its inherent complexity associated with elusive intermediates hinders the predictive synthesis of high‐performance cathode materials. Here, correlative in situ X‐ray absorption/scattering spectroscopy is used to investigate the calcination of nickel‐based cathodes, focusing specifically on the archetypal LiNiO 2 from Ni(OH) 2 . Combining in situ observation with data‐driven analysis reveals concurrent lithiation and dehydration of Ni(OH) 2 and consequently, the low‐temperature crystallization of layered LiNiO 2 alongside lithiated rocksalts. Following early nucleation, LiNiO 2 undergoes sluggish crystallization and structural ordering while depleting rocksalts; ultimately, it turns into a structurally‐ordered layered phase upon full lithiation but remains small in size. Subsequent high‐temperature sintering induces rapid crystal growth, accompanied by undesired delithiation and structural degradation. These observations are further corroborated by mesoscale modeling, emphasizing that, even though calcination is thermally driven and favors transformation towards thermodynamically equilibrium phases, the actual phase propagation and crystallization can be kinetically tuned via lithiation, providing freedom for structural and morphological control during cathode calcination.

36 MATERIALS SCIENCE↗

Synthesizing intermetallic thermoelectric materials using metal flux growth method

With zero energy poverty by 2050 goals, exploring innovative energy technologies is essential. Thermoelectric materials (TEMs) have the potential to harness waste heat and convert it to electricity (e.g., for spacecraft power, sensors, and cryogenic cooling). Two limitations of thermoelectric compounds are their low efficiency and limited pool of materials. Optimizing the efficiency of TEMs is challenging due to the interdependent relationship between electrical and thermal conductivity. A deeper understanding of electronic properties and lattice thermal conductivity is necessary to “understand, predict and ultimately control” the thermoelectric properties of TEMs. Understanding the influence of light-atom-doping and intrinsic defects on the electronic properties of the f-element containing TEMs are important questions that need to be answered to design and synthesize TEMs with optimal parameters. The purpose of this work is to understand how doping f-element based TEMs with light elements, such as carbon and boron, would impact the crystal structure, bond lengths, band gaps, and thermoelectric characteristics. High-quality f¬-element containing intermetallic single crystals are grown in metal flux media. In molten metal flux growth method, one or more low-melting metals such as aluminum, gallium, tin, and bismuth are used in excess as the reaction medium. As these are solution-state reactions, we can use lower temperatures than solid-state reactions and isolate the kinetic products instead of the thermodynamically stable compound. Then these are characterized using X-ray diffraction, scanning electron microscopy, energy dispersive spectroscopy, magnetization, electrical and thermal conductivities, and thermoelectric power.

37 - INORGANIC, ORGANIC, PHYSICAL AND ANALYTICAL C↗

tinyIFD: A High-Throughput Binding Pose Refinement Workflow Through Induced-Fit Ligand Docking

A critical step in structure-based drug discovery is predicting whether and how a candidate molecule binds to a model of a therapeutic target. However, substantial protein side chain movements prevent current screening methods, such as docking, from accurately predicting the ligand conformations and require expensive refinements to produce viable candidates. Here, we present the development of a high-throughput and flexible ligand pose refinement workflow, called “tinyIFD”. The main features of the workflow include the use of specialized high-throughput, small-system MD simulation code mdgx.cuda and an actively learning model zoo approach. We show the application of this workflow on a large test set of diverse protein targets, achieving 66% and 76% success rates for finding a crystal-like pose within the top-2 and top-5 poses, respectively. We also applied this workflow to the SARS-CoV-2 main protease (M pro ) inhibitors, where we demonstrate the benefit of the active learning aspect in this workflow.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Workhorse minimally empirical dispersion-corrected density functional with tests for weakly bound systems: r 2 SCAN + rVV 10

SCAN+rVV10 has been demonstrated to be a versatile van der Waals (vdW) density functional that delivers good predictions of both energetic and structural properties for many types of bonding. Recently, the r 2 SCAN functional was devised as a revised form of SCAN with improved numerical stability. In this work, we refit the rVV10 functional to optimize the r 2 SCAN+rVV10 vdW density functional and test its performance for molecular interactions and layered materials. Our molecular tests demonstrate that r 2 SCAN+rVV10 outperforms its predecessor SCAN+rVV10 in both efficiency (numerical stability) and accuracy. This good performance is also found in lattice-constant predictions. In comparison with benchmark results from higher-level theories or experiments, r 2 SCAN+rVV10 yields excellent interlayer binding energies and phonon dispersions for layered materials.

36 MATERIALS SCIENCE↗

Direct Prediction of Phonon Density of States With Euclidean Neural Networks

Abstract Machine learning has demonstrated great power in materials design, discovery, and property prediction. However, despite the success of machine learning in predicting discrete properties, challenges remain for continuous property prediction. The challenge is aggravated in crystalline solids due to crystallographic symmetry considerations and data scarcity. Here, the direct prediction of phonon density‐of‐states (DOS) is demonstrated using only atomic species and positions as input. Euclidean neural networks are applied, which by construction are equivariant to 3D rotations, translations, and inversion and thereby capture full crystal symmetry, and achieve high‐quality prediction using a small training set of examples with over 64 atom types. The predictive model reproduces key features of experimental data and even generalizes to materials with unseen elements, and is naturally suited to efficiently predict alloy systems without additional computational cost. The potential of the network is demonstrated by predicting a broad number of high phononic specific heat capacity materials. The work indicates an efficient approach to explore materials' phonon structure, and can further enable rapid screening for high‐performance thermal storage materials and phonon‐mediated superconductors.

97 MATHEMATICS AND COMPUTING↗