Search NASA⌕ Search

SEARCH · Search NASA

Results for “ionic liquid”

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 307 records · Page 17

Recent Improvements to the NWChem COSMO Module

This work presents recent developments in the COSMO solvation model implementation in NWChem. A new cavity construction approach, based on the solvent-excluding surface (SES) and utilizing the well-established GEPOL algorithm, has been introduced. Additionally, a straightforward procedure to merge surface segments that are too close─often a source of numerical artifacts─has been implemented. The available methods for correcting outlying surface charges have also been reviewed and improved. To validate the new implementation, we computed dielectric solvation energies for a chemically diverse set of approximately 100 molecules, including neutral species, small ions, and common ionic liquid components. Results were compared to those from GAMESS using the double-cavity method as a reference. Although the double-cavity approach can be regarded as more accurate, the simpler correction schemes available in NWChem─based on scaling factors or Lagrange multipliers─can achieve excellent agreement if the potential is also properly corrected, with mean unsigned deviations of around 0.15 kcal/mol. Predictions of typical vapor–liquid and liquid–liquid equilibria using a COSMO-SAC variant based on NWChem also showed very promising results.

Cavities↗

Intermolecular Proton Transfer Enabled Reactive CO 2 Capture by the Malononitrile Anion

Task-specific ionic liquids (ILs) employing carbanions represent a new class of ILs for carbon capture. The deprotonated malononitrile carbanion, [CH(CN) 2 ] - , has shown close to equimolar capacity for reactive CO 2 capture. Although the formation of the [C(CN) 2 COOH] - carboxylic acid was found to be the final product, how the hydrogen atom on the [CH(CN) 2 ] - carbanion transfers to the carboxylate group as a proton has not been fully understood. In this work, we employ density functional theory calculations with an implicit solvation model to investigate the proton transfer mechanisms in forming carboxylic acid from the reaction of the [CH(CN) 2 ] - carbanion with CO 2 . We find that the intramolecular proton-transfer pathway in [CH(CN) 2 COO] - to form [C(CN) 2 COOH] - is unlikely due to the high energy barrier of 152 kJ/mol. Instead, the intermolecular proton transfer pathway between two [CH(CN) 2 COO] - anions is more feasible to form two molecules of [C(CN) 2 COOH] - , with a significantly lower activation energy of 50 kJ/mol. Moreover, the [C(CN) 2 COOH] - dimer is further stabilized by the intermolecular hydrogen bonds of the two –COOH groups in the Z-configuration of the π-conjugated planar geometry. This insight of reactive CO 2 capture enabled by intermolecular proton transfer will be useful in designing novel carbanions and ILs for carbon capture and conversion.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Catalytic Promotion of Transition-Metal-Doped Graphene Cathodes in Li-CO 2 Batteries

The Li-CO 2 battery is a promising energy storage system with impressive theoretical specific energy and discharge capacity. Graphene-based single-atom catalysts (SACs) provide high surface area and long-term electrochemical reactivity and stability, making SACs among the most promising cathode catalysts for these batteries. However, current Li-CO 2 systems have high reaction barriers, slowing the reaction and greatly increasing the overpotential. Improvement of the discharge/charge energetics requires atomic-level innovations in cathode design, such as alterations to the catalyst chemical structure. In this paper, we propose enhancing the SAC by using a Ti metal center, which is found to deliver the highest electrochemical Li + CO 2 activity among 3d transition metal candidates. Furthermore, we propose cathode surface coating with ionic liquids, since these environments promote the formation of reaction intermediates in the electrochemical conversion process. Here, our work provides insights to optimize electrode design for high-performance Li-CO 2 batteries, which can open new avenues to recycle greenhouse gases and achieve enhanced renewable energy storage.

25 ENERGY STORAGE↗

Challenges in Pulsed-Field Gradient Nuclear Magnetic Resonance on Magnetically Heterogeneous Interfaces: Sequence and Field-Dependent Apparent Diffusion Coefficients

It is well known that the internal gradient (gi) that exists within pores haunts the diffusion coefficient (D) as measured by the pulsed-field gradient (PFG) nuclear magnetic resonance (NMR). Several PFG-NMR methods developed to determine an accurate D were not successful. Then, the steady-state diffusion coefficient (Dapp,8) for the cation [C4mim]+ of [C4mim][Tf2N]; [1-butyl-3-methylimidazolium][bis(trifluoromethylsulfonyl)imde] ionic liquid confined in ordered mesoporous carbon (OMC) were determined by comparing Dapp,8 obtained from 1H PFG-NMR performed with three different stimulated echo sequences: STE, APFG, and MPFG under the two external magnetic field strength, B0 = 9.4 and 14.1 Tesla. The measured Dapp,8 which is an order of magnitude smaller than D of bulk [C4mim][Tf2N], is in good agreement between APFG and MPFG both in B0 = 9.4 and 14.1 Tesla. However, the strong gi artifact, which caused apparent diffusion coefficient (Dapp) depending strongly and weakly on B0 and temperature, respectively, in diffusion-time dependent Dapp, Dapp(?) obtained from a sequence with monopolar gradients (STE) was suppressed by using sequences employing bipolar gradients (APFG and MPFG) in the region of steady-state diffusion. But incompletely suppressed gi artifact resulting in the different behaviors of the early part of Dapp(?) between the sequences leads a ˜ 0.6 and 0.9 in MPFG and APFG, respectively, in the relationship between mean squared displacement and diffusion time: = 2Dta, where a = 0.5 and 1 for 1-dimensional single file diffusion and 3-dimensional bulk diffusion, respectively. The above observations clearly show that the diffusion behavior of ions/molecules within the pores and pore structure, such as the surface-to-volume ratio? (D?_app (?)=D_0 [1-4/(9vp) S/V v(D_0 ?)]) and tortuosity (T = D0/Dapp,8), are possible to be misunderstood, especially in the systems with a non-negligible gi. This work demonstrates that it may be necessary to test several PFG sequences under multiple external magnetic fields for the correct determination of the diffusion behavior of ions/molecules in the pores with a larger internal gradient, gi.

Han, Kee Sung↗

Hybrid Bismuth Halide with Rich Polymorphism and Second Harmonic Generation Response

Hybrid structures have emerged as a promising class of optical materials, due to their ability to couple the robustness of inorganic and the tunability of organic compounds. However, their application in nonlinear optics (NLO) remains limited, largely due to underexplored factors that drive the formation of noncentrosymmetric structures and NLO property characterization. In this work, we explore the formation, structural and temperature polymorphism, and optical properties of the (Et 3 NH) 3 Bi 2 Br 9 composition, which crystallizes as either noncentrosymmetric or centrosymmetric polymorph. Structural analysis showed that the alignment of [Bi 2 Br 9 ] 3– units dictates the symmetry of phases, as well as the nonlinear optical properties, with the triclinic polymorph exhibiting a second harmonic generation (SHG) response both in visible and IR regions (1.29 × KH 2 PO 4 and 0.08 × AgGaS 2 ). Thermal analysis reveals polymorphic phase transitions and low melting points, making them melt-processable and ionic liquid candidates.

36 MATERIALS SCIENCE↗

Ambient Rare Earth Metal Electrodeposition in Nitrogen-Coordinated Silylamide Electrolyte

The transition to a sustainable, low-carbon economy demands energy-efficient and environmentally benign methods for rare earth metal (REM) production. Furthermore, while high-temperature molten salt electrolysis remains energy-intensive, corrosive, and environmentally unfriendly, emerging room-temperature processes based on conventional ionic liquids are also hindered by high viscosity and chemical instability. In this study, a nitrogen-coordinated, water-, oxygen-, and fluorine-free silylamide-based electrolyte is presented as a promising system for room-temperature REM electrodeposition. Derived from commercially available lithium silylamide precursors, the system enables facile synthesis, broad electrochemical stability, and tunable metal–ligand interactions. Electrochemical analysis reveals high voltammetric stripping reversibility, stable cathodic and anodic potentials, and selective neodymium (Nd) deposition at appreciable current densities (>1 mA/cm 2 ), with minimal parasitic reactions. Bulk experiments produced high-purity Nd with reproducible performance across multiple batches. Scaled deposition yielded over 1 g of Nd with near-theoretical mass efficiency (0.40 mg/C) and >90% purity. Using sacrificial dysprosium (Dy) and Nd metal anodes, the system maintained constant Nd loading during extended deposition and enabled cathodic codeposition of stripped anode material, demonstrating the electrolyte’s dual functionality for REM electrorefining. Collectively, this silylamide platform offers a compelling combination of electrochemical robustness, chemical resilience, and process scalability for sustainable ambient REM recovery.

36 - MATERIALS SCIENCE↗

Synergistic Effects of the Electric Field Induced by Imidazolium Rotation and Hydrogen Bonding in Electrocatalysis of CO 2

The roles of the ionic liquid, 1-ethyl-3-methylimidazolium tetrafluoroborate ([EMIM][BF₄]), and water in controlling the mechanism, energetics, and electrocatalytic activity of CO₂ reduction to CO on silver in non-aqueous electrolytes were investigated. The first electron transfer occurs to CO₂ at reduced overpotentials when it is trapped between the planes of [EMIM]⁺ ring and electrode surface due to cation re-orientation as determined from voltammetry, in-situ surface-enhanced Raman spectroscopy, and density functional theory calculations. Within this interface, water up to 0.5 M does not induce significant faradaic activity, opposing the notion of it as a free proton source. Instead, water acts as a hydrogen bond donor and the proton is sourced from [EMIM]⁺. Furthermore, this study demonstrates that alcohols with varying acidity tune the hydrogen bonding network in the interfacial microenvironment to lower the energetics required for CO₂ reduction. Here, the hydrogen bonding suppresses the formation of inactive carboxylate species, thus preserving the catalytic activity of [EMIM]⁺. The ability to tune the hydrogen bonding network opens new avenues for advancing IL-mediated electrocatalytic reactions in non-aqueous electrolytes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Catalytic Borylation of Poly(vinyl chloride) Produces Adhesive Materials

Postpolymerization functionalization of polymers can create new applications for existing materials, while retaining their most favorable, intrinsic properties. Polyvinyl chloride (PVC) is a widely used, commodity polymer that is particularly challenging to modify. We report a copper-catalyzed protocol that replaces a small fraction of the C-Cl bonds in PVC with C-B bonds to boronic esters. The reaction occurs with an inexpensive catalyst comprising copper(II) chloride and an NHC ligand derived from a common ionic liquid that is distinct from the N-heterocyclic carbene (NHC) used for the borylation of small alkyl halides. The resulting materials adhere strongly to common surfaces, such as glass and metals, even more strongly than do commercial glues.

D’Angelo, Kyan A↗

Thermodynamics and transport in molten chloride salts and their mixtures

Relationship between thermophysical properties and phonon mean free path. Heat capacity, viscosity, and thermal conductivity in ionic liquids decrease as mean free path decreases and dynamics become less “solid-like” and collect motion diminishes.

Cockrell, C. [Nuclear Futures Institute, Bangor Un↗

Consistent and reproducible computation of the glass transition temperature from molecular dynamics simulations

In many fields, from semiconductors for opto-electronic applications to ionic liquids (ILs) for separations, the glass transition temperature (Tg) of a material is a useful gauge for its potential use in practical settings. As a result, there is a great deal of interest in predicting Tg using molecular simulations. However, the uncertainty and variation in the trend shift method, a common approach in simulations to predict Tg, can be high. This is due to the need for human intervention in defining a fitting range for linear fits of density with temperature assumed for the liquid and glass phases across the simulated cooling. The definition of such fitting ranges then defines the estimate for the Tg as the intersection of linear fits. We eliminate this need for human intervention by leveraging the Shapiro–Wilk normality test and proposing an algorithm to define the fitting ranges and, consequently, Tg. Through this integration, we incorporate into our automated methodology that residuals must be normally distributed around zero for any fit, a requirement that must be met for any regression problem. Consequently, fitting ranges for realizing linear fits for each phase are statistically defined rather than visually inferred, obtaining an estimate for Tg without any human intervention. The method is also capable of finding multiple linear regimes across density vs temperature curves. We compare the predictions of our proposed method across multiple IL and semiconductor molecular dynamics simulation results from the literature and compare other proposed methods for automatically detecting Tg from density–temperature data. We believe that our proposed method would allow for more consistent predictions of Tg. We make this methodology available and open source through GitHub.

Chemistry↗

Unconventional gate-induced superconductivity in transition-metal dichalcogenides

Superconductivity in few-layer semiconducting transition-metal dichalcogenides (TMDs) can be induced by field-effect doping through ionic-liquid gating. While several experimental observations have been collected over the years, a fully consistent theoretical picture is still missing. Here we develop a realistic framework that combines the predictive power of first-principles simulations with the versatility and insight of Bardeen-Cooper-Schrieffer gap equations to rationalize such experiments. The multivalley nature of semiconducting TMDs is taken into account, together with the doping- and momentum-dependent electron-phonon and Coulomb interactions. Consistently with experiments, we find that superconductivity occurs when the electron density is large enough that the 𝑄 valleys get occupied, as a result of a large enhancement of electron-phonon interactions. Despite being phonon driven, the superconducting state is predicted to be sensitive to Coulomb interactions, which can lead to the appearance of a relative sign difference between valleys and thus to a 𝑠 +− character. We discuss qualitatively how such scenario may account for many of the observed physical phenomena for which no microscopic explanation has been found so far, including in particular the presence of a large subgap density of states, and the sample-dependent dome-shaped dependence of 𝑇 𝑐 on accumulated electron density. Our results provide a comprehensive analysis of gate-induced superconductivity in semiconducting TMDs, and introduce an approach that will likely be valuable for other multivalley electronic systems, in which superconductivity occurs at relatively low electron density.

2-dimensional systems↗

Integrated low-temperature PVC and polyolefin upgrading

Polyolefins and their chlorinated derivatives such as polyvinyl chloride (PVC) are among the most prevalent plastics in global production and waste streams. Traditional waste-to-energy methods such as incineration and pyrolysis, as well as most chemical upcycling methods for PVC utilization, require thorough, high-temperature dechlorination to prevent the release of toxic chlorinated compounds. Here, we present here a strategy for upgrading discarded PVC into chlorine-free fuel range hydrocarbons and hydrogen chloride in a single-stage process catalyzed by chloroaluminate ionic liquids. This approach offsets endothermic dechlorination and carbon-carbon bond cleavage with exothermic alkylation and hydrogen transfer by isobutane or isopentane in a low-temperature tandem process. The light isoalkanes are available from refinery processes and partly from recycling of the product stream. This process is suitable for handling real-world mixed and contaminated PVC and polyolefin waste streams.

Zhang, Wei [East China Normal Univ. (ECNU), Shangh↗

Progress and Perspectives: Zirconium Electrodeposition from Different Electrolytes

Zirconium (Zr) possesses outstanding properties, including exceptional chemical resistance and a high melting point, and is therefore desirable for use in a wide variety of challenging environments, such as in nuclear reactors and the chemical processing industry. To minimize the amount of Zr required for any given application, developing methods for generating metallic Zr coatings is highly advantageous. Owing to Zr’s highly negative reduction potential, electrodeposition in traditional solvents at near ambient temperatures remains challenging and thus not well understood. Due to the extreme conditions required to deposit this metal, extensive work has been conducted in molten salt electrolytes, however the broad applicability of this methodology is limited due to its corrosivity. The primary focus of this article is to present an overview of Zr’s electrochemical behavior and to consolidate the efforts of researchers in exploring electrodeposition techniques for Zr involving aqueous, organic, ionic liquid, deep eutectic, and molten salt solvents. With this information, we highlight trends across solvent systems and opportunities for future research.

42 ENGINEERING↗

Science Driven Discovery of Nanoparticle Photocatalysts

Photocatalytic water splitting using suspensions of nanoparticle photocatalysts is a promising route to economically sustainable production of green hydrogen. The principal challenge is to develop photocatalysts with overall solar-to-hydrogen conversion efficiency that exceeds 10 percent. In this project we have developed a new platform for investigating candidate materials for photocatalytic water splitting. Our platform consists of patterned Au electrodes and a Ag/AgCl reference electrode on an insulating substrate onto which we disperse nanoparticle photocatalysts. We then cover the substrate with a thin layer of ionogel containing a protic ionic liquid that dissolves water from the ambient. Using this platform we have demonstrated photoelectrochemical activity mapping for single and small clusters of BiVO 4 nanoparticle photocatalysts and correlated these results to their Raman and photoluminescence spectra. The preliminary results suggest a strong correlation for low efficiency nanoparticles, followed by saturation for those with higher activities, indicating that interface reaction or electrolyte transport become the limiting factor. We anticipate that further application of this platform to investigation of candidate photocatalyst materials will provide useful insights into the mechanisms that limit their performance.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

Microchannel-based Membrane-less Extraction of Li from Unconventional Lithium Sources & the Separation of REE

This final report provides an overview of the Project's entire duration, covering July 1, 2021 to December 31, 2023. It primarily focuses on the achievements, technological developments, and unique challenges the team faced while working on separating and extracting Lithium from produced waters. The project's primary aim was to create an integrated, high-throughput, membrane-less, and modular microfluidic platform that could extract Lithium from unconventional sources. We have successfully met all goals and milestones envisioned in the SOPO document. The most critical primary milestones, including the Go-No-Go milestone (refer to the Gantt chart in the Appendices), were successfully accomplished. We demonstrated phase separation (>90%) and extraction (>85%) performance in the MPSE using synthetic, and representative produced water composition feed at 50 ml/min total flow through MPSE 36. We have also performed a parametric study of the MPSE operations, beyond the scope of SOPO, exploring operating conditions of current and broader interest. The extended investigation of operational parameters is concurrent with our efforts to seek further development of the MPSE technology beyond the scope of the Project. Along these lines of development, we have made efforts to be responsive to DOE calls for technological developments of other types of resources (beyond PW) for the recovery of Critical Materials and higher TRL development (beyond TRL 4). During the work on this Project, we developed and implemented three innovative technical approaches that emerged from our efforts to successfully meet the Project milestones. The innovative & original technical approaches developed and implemented in this Project are now the contributions to process engineering that could be clearly credited to the Project. First, Convergent Design Approach is a comprehensive feedforward & feedback loop of four design phases: i) design for functionality, ii) design for manufacturing, iii) design for sustainability, and iv) design for market. Next was Process Intensification. A major aim of this Project was to create an innovative phase separation & extraction microscale-based technology for Li separation – thus the words microchannel-based in the Project title. A microscale-based technology is intrinsically in the center of the Process Intensification domain as defined by its unique principles. Therefore, Process Intensification was implicitly envisioned in the Project’s SOPO. Lastly, Time Scale Analysis is a novel tool for discovering the needs and directions of Process Intensification implementations in any process technology. This Project is fully credited for developing and implementing the three novel technical approaches mentioned above. These are general contributions to process engineering that emerged from this Project. Beyond the original SOPO scope, the OSU-U.Pitt research group utilized a Convergent Design methodology, integrating first-principles mathematical modeling with experimental validation on the Minimum Development Vehicle. By creating these Digital Twins, the team rapidly assessed manufacturing iterations to support TEA analysis. This framework further enabled the development of advanced Surface Modification Techniques, where hydrophobic and oleophobic coating strategies were optimized via Digital Twin tools and validated through rigorous 100-hour longevity testing. TEA Analysis: The closing efforts of this Project were focused on the TEA analysis. TEA analysis had two primary functions: i) enabling critical assessments of design variations withing 10 the Concurrent Design Approach, thus enabling evolution of the MPSE design to reach faster- better-cheaper alternatives; and ii) to create a bridge between the accomplishments of this Project and future projects of higher TRL, beyond TRL 6 level. It is important to note that the TEA model created in the Project stirred the technological solutions for the recovery of critical materials toward a vision of a very profitable modular plant that has unique zero-waste water discharge signature. More importantly, thanks to our experimental performance data and conservative assumptions, the TEA model predicts minimal technological and investment risks. Low cost of a modular unit of a nominal capacity of [1000 tons of Li 2 CO 3 /year] positions the MPSE based technology within the reach of community investors, thus offering a paradigm shift in the development of critical technologies. The project successfully navigated two primary challenges: solvent selection and manufacturing adaptation. Restricted by the SOPO to existing literature for lithium recovery, the team identified a critical need for a "material excellence program" to develop next-generation solvents, eventually concluding with a preliminary investigation into promising Ionic Liquids (ILs). Simultaneously, COVID-19 supply chain disruptions forced a pivot from traditional manufacturing to advanced additive methods at ATAMI-OSU. By transitioning from stainless steel to 3D-printed polymer substrates, the team achieved a transformative three-order-of- magnitude reduction in manufacturing costs and compressed prototyping timelines from several months to just two days. The MPSE technology offers significant energy, environmental, and economic advantages by overcoming the traditional bottlenecks of phase-separation hardware and contactor size. Unlike conventional mixer-settlers or membrane-based systems, MPSE operates without moving parts or fouling-prone membranes, achieving robust performance even with challenging, viscous, or particulate-heavy feeds. Key performance metrics include an energy intensity reduction of 5–50x (3–40 kJ/m 3 ) compared to incumbent technologies and a dramatic reduction of processing time to under 60 seconds, which drastically reduces the physical plant footprint. These technical efficiencies translate into superior economic outcomes; for a 100 t/year Li 2 CO 3 facility, implementing MPSE is projected to nearly halve contactor CAPEX (from $\$$6.08M to $\$$3.01M) and significantly increase the project's Net Present Value (NPV), derisking new investment and enabling distributed critical-mineral processing configurations. The commercialization of MPSE technology is being spearheaded by Vigsur Dynamics Inc., which has adopted a structured, parallel approach to technical and business development since its formation in January 2026. Following extensive customer discovery and engagement with the Oregon State University accelerator, Vigsur Dynamics is working to establish a business model that transitions from pilot demonstrations to modular hardware sales, ultimately aiming for a "build-own-operate" service strategy. Current technical milestones—including 100 hours of continuous operation, superior energy efficiency, and successful 6-unit modular scale-up— provide a foundation for this transition. Backed by ongoing IP licensing and a growing network of industrial and venture advisors, the company is actively de-risking the platform to replace conventional mixer-settler systems in the critical minerals market.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Development of a New Aggregation Method to Remove Nanoplastics from the Ocean: Proof of Concept Using Mussel Exposure Tests

The overproduction and mismanagement of plastics has led to the accumulation of these materials in the environment, particularly in the marine ecosystem. Once in the environment, plastics break down and can acquire microscopic or even nanoscopic sizes. Given their sizes, microplastics (MPs) and nanoplastics (NPs) are hard to detect and remove from the aquatic environment, eventually interacting with marine organisms. This research mainly aimed to achieve the aggregation of micro- and nanoplastics (MNPs) to ease their removal from the marine environment. To this end, the size and stability of polystyrene (PS) MNPs were measured in synthetic seawater with the different components of the technology (ionic liquid and chitosan). The MPs were purchased in their plain form, while the NPs displayed amines on their surface (PS NP-NH2). The results showed that this technology promoted a significant aggregation of the PS NP-NH2, whereas, for the PS MPs, no conclusive results were found, indicating that the surface charge plays an essential role in the MNP aggregation process. Moreover, to investigate the toxicological potential of MNPs, a mussel species (M. galloprovincialis) was exposed to different concentrations of MPs and NPs, separately, with and without the technology. In this context, mussels were sampled after 7, 14, and 21 days of exposure, and the gills and digestive glands were collected for analysis of oxidative stress biomarkers and histological observations. In general, the results indicate that MNPs trigger the production of reactive oxygen species (ROS) in mussels and induce oxidative stress, making gills the most affected organ. Yet, when the technology was applied in moderate concentrations, NPs showed adverse effects in mussels. The histological analysis showed no evidence of MNPs in the gill’s tissues.

Cid-Samamed, Antonio (ORCID:0000000205073394)↗

Comment on 'Formation of Venusian canali': Considerations of lava types and their thermal behaviors' by T. K. P. Gregg and R. Greeley

Gregg and Greeley (1993) evaluated the mechanisms by which various magma types might have formed the long sinuous channels apparent on the surface of Venus. Among the magma types they evaluated were carbonatites, which are carbonate-rich ionic liquids. Gregg and Greeley (1993) applied thermophysical data for terrestrial carbonatite magmas to evaluate channel formation by carbonatites on Venus. Their analysis is flawed because they used incorrect values for the solidification temperatures of carbonatite magmas.

Treiman, Allan H.↗