1
|
Welsh BA, Urbina AS, Ho TA, Rempe SL, Slipchenko LV, Zwier TS. Capturing CO 2 in Quadrupolar Binding Pockets: Broadband Microwave Spectroscopy of Pyrimidine-(CO 2) n, n = 1,2. J Phys Chem A 2024; 128:1124-1133. [PMID: 38306293 DOI: 10.1021/acs.jpca.3c07930] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/04/2024]
Abstract
Pyrimidine has two in-plane CH(δ+)/N̈(δ-)/CH(δ+) binding sites that are complementary to the (δ-/2δ+/δ-) quadrupole moment of CO2. We recorded broadband microwave spectra over the 7.5-17.5 GHz range for pyrimidine-(CO2)n with n = 1 and 2 formed in a supersonic expansion. Based on fits of the rotational transitions, including nuclear hyperfine splitting due to the two 14N nuclei, we have assigned 313 hyperfine components across 105 rotational transitions for the n = 1 complex and 208 hyperfine components across 105 rotational transitions for the n = 2 complex. The pyrimidine-CO2 complex is planar, with CO2 occupying one of the quadrupolar binding sites, forming a structure in which the CO2 is stabilized in the plane by interactions with the C-H hydrogens adjacent to the nitrogen atom. This structure is closely analogous to that of the pyridine-CO2 complex studied previously by (Doran, J. L. J. Mol. Struct. 2012, 1019, 191-195). The fit to the n = 2 cluster gives rotational constants consistent with a planar cluster of C2v symmetry in which the second CO2 molecule binds in the second quadrupolar binding pocket on the opposite side of the ring. The calculated total binding energy in pyrimidine-CO2 is -13.7 kJ mol-1, including corrections for basis set superposition error and zero-point energy, at the CCSD(T)/ 6-311++G(3df,2p) level, while that in pyrimidine-(CO2)2 is almost exactly double that size, indicating little interaction between the two CO2 molecules in the two binding sites. The enthalpy, entropy, and free energy of binding are also calculated at 300 K within the harmonic oscillator/rigid-rotor model. This model is shown to lack quantitative accuracy when it is applied to the formation of weakly bound complexes.
Collapse
Affiliation(s)
- Blair A Welsh
- Gas Phase Chemical Physics, Sandia National Laboratories, Livermore, California 94550, United States
| | - Andres S Urbina
- Department of Chemistry, Purdue University, West Lafayette, Indiana 47907-1393, United States
| | - Tuan A Ho
- Geochemistry Department, Sandia National Laboratories, Albuquerque, New Mexico 87185, United States
| | - Susan L Rempe
- Center for Integrated Nanotechnologies, Sandia National Laboratories, Albuquerque, New Mexico 87185, United States
| | - Lyudmila V Slipchenko
- Department of Chemistry, Purdue University, West Lafayette, Indiana 47907-1393, United States
| | - Timothy S Zwier
- Gas Phase Chemical Physics, Sandia National Laboratories, Livermore, California 94550, United States
| |
Collapse
|
2
|
Hernandez-Castillo AO, Calabrese C, Fritz SM, Uriarte I, Cocinero EJ, Zwier TS. Bond Length Alternation and Internal Dynamics in Model Aromatic Substituents of Lignin. Chemphyschem 2022; 23:e202100808. [PMID: 35102679 DOI: 10.1002/cphc.202100808] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/12/2021] [Revised: 12/20/2021] [Indexed: 11/07/2022]
Abstract
Broadband chirped-pulse microwave spectra were recorded over the 2-18 GHz frequency range for a series of four model aromatic components of lignin; namely, guaiacol ( ortho -methoxy phenol, G ), syringol (2,6-dimethoxy phenol, S ), 4-methyl guaiacol ( MG ), and 4-vinyl guaiacol ( VG ), under jet-cooled conditions in the gas phase. Using a combination of 13 C isotopic data and electronic structure calculations, distortions of the phenyl ring by the substituents on the ring are identified. In all four molecules, the r C(1)-C(6) bond between the two substituted C-atoms lengthens, leading to clear bond alternation that reflects an increase in the phenyl ring resonance structure with double bonds at r C(1)-C(2) , r C(3)-C(4) and r C(5)-C(6) . Syringol, with its symmetric methoxy substituents, possesses a microwave spectrum with tunneling doublets in the a -type transitions associated with H-atom tunneling. These splittings were fit to determine a barrier to hindered rotation of the OH group of 1975 cm -1 , a value nearly 50% greater than that in phenol, due to the presence of the intramolecular OH … OCH 3 H-bonds at the two equivalent planar geometries. In 4-methyl guaiacol, methyl rotor splittings are observed and used to confirm and refine an earlier measurement of the three-fold barrier V 3 = 67 cm -1 . Finally, 4-vinyl guaiacol shows transitions due to two conformers differing in the relative orientations of the vinyl and OH groups.
Collapse
Affiliation(s)
- A O Hernandez-Castillo
- Department of Chemistry, Purdue University, West Lafayette, IN 47907-2084, USA
- Department of Molecular Physics, Fritz Haber Institute of the Max Planck Society, Faradayweg 4-6, Berlin, Germany, 14195
| | - Camilla Calabrese
- Department of Physical Chemistry, Faculty of Science and Technology, University of the Basque Country (UPV/EHU), B° Sarriena, S/N, Leioa, 48940, Spain
- Instituto Biofisika (UPV/EHU-CSIC), University of the Basque Country (UPV/EHU), B° Sarriena S/N, Leioa, 48940, Spain
- Fundación Biofísica Bizkaia/Biofisika Bizkaia Fundazioa (FBB), E-48940, Leioa, Spain
| | - Sean M Fritz
- Department of Chemistry, Purdue University, West Lafayette, IN 47907-2084, USA
| | - Iciar Uriarte
- Department of Physical Chemistry, Faculty of Science and Technology, University of the Basque Country (UPV/EHU), B° Sarriena, S/N, Leioa, 48940, Spain
| | - Emilio J Cocinero
- Department of Physical Chemistry, Faculty of Science and Technology, University of the Basque Country (UPV/EHU), B° Sarriena, S/N, Leioa, 48940, Spain
- Instituto Biofisika (UPV/EHU-CSIC), University of the Basque Country (UPV/EHU), B° Sarriena S/N, Leioa, 48940, Spain
| | - Timothy S Zwier
- Department of Chemistry, Purdue University, West Lafayette, IN 47907-2084, USA
- Gas Phase Chemical Physics, Livermore, CA 94550, USA
| |
Collapse
|
3
|
Seifert NA, Prozument K, Davis MJ. Computational optimal transport for molecular spectra: The fully discrete case. J Chem Phys 2021; 155:184101. [PMID: 34773953 DOI: 10.1063/5.0069681] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The use of computational optimal transport is investigated as a tool for comparing two molecular spectra. Unlike other techniques for comparing molecular spectra in a pattern-recognition framework, transport distances simultaneously encode information about line positions and intensities. In addition, it is shown that transport distances are a useful alternative to Euclidean distances as Euclidean distances are based on line-by-line comparisons, while transport distances reflect broader features of molecular spectra and adequately compare spectra with different resolutions. This paper includes a tutorial on the use of optimal transport and investigates several well-chosen examples to illustrate the utility of computational optimal transport for comparing molecular spectra.
Collapse
Affiliation(s)
- Nathan A Seifert
- Chemical Sciences and Engineering Division, Argonne National Laboratory, Lemont, Illinois 60439, USA
| | - Kirill Prozument
- Chemical Sciences and Engineering Division, Argonne National Laboratory, Lemont, Illinois 60439, USA
| | - Michael J Davis
- Chemical Sciences and Engineering Division, Argonne National Laboratory, Lemont, Illinois 60439, USA
| |
Collapse
|
4
|
Mishra P, Fritz SM, Herbers S, Mebel AM, Zwier TS. Gas-phase pyrolysis of trans 3-pentenenitrile: competition between direct and isomerization-mediated dissociation. Phys Chem Chem Phys 2021; 23:6462-6471. [PMID: 33729262 DOI: 10.1039/d1cp00104c] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/25/2023]
Abstract
The flash pyrolysis of trans 3-pentenenitrile (3-PN, CH3-CH[double bond, length as m-dash]CH-CH2-CN) was studied by combining the results of VUV photoionization mass spectra with broadband microwave spectra recorded as a function of the temperature of the pyrolysis tube. The two separated functional groups (vinyl and nitrile) open up isomerization as an initial step in competition with unimolecular dissociation. Primary products were detected by keeping the 3-PN concentration low and limiting reaction times to the traversal time of the gas in the pyrolysis tube (∼100 μs). The reaction is quenched and products are cooled by expansion into vacuum before interrogation over the 8-18 GHz region using chirped-pulse broadband methods. 118 nm VUV photoionization of the same reaction mixture provides a means of detecting all products with ionization potentials below 10.5 eV with minimal fragmentation. These results are combined with a detailed computational investigation of the C5H7N and related potential energy surfaces, leading to a consistent picture of the unimolecular decomposition of 3-PN. Loss of two H-atoms to form a 79 amu product is proven from its microwave transitions to contain trans-Z-2,4-pentadienenitrile, while no pyridine is observed. Methyl loss, HCN loss, and breaking the central C(2)-C(3) bond all occur following isomerization of the position of the double bond, thereby opening up low-energy pathways to these decomposition channels.
Collapse
Affiliation(s)
- Piyush Mishra
- Department of Chemistry, Purdue University, West Lafayette, IN 47907-1393, USA
| | | | | | | | | |
Collapse
|
5
|
Herbers S, Fritz SM, Mishra P, Nguyen HVL, Zwier TS. Local and global approaches to treat the torsional barriers of 4-methylacetophenone using microwave spectroscopy. J Chem Phys 2020; 152:074301. [PMID: 32087663 DOI: 10.1063/1.5142401] [Citation(s) in RCA: 16] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Sven Herbers
- Department of Chemistry, Purdue University, 560 Oval Drive, West Lafayette, Indiana 47907, USA
| | - Sean M. Fritz
- Department of Chemistry, Purdue University, 560 Oval Drive, West Lafayette, Indiana 47907, USA
| | - Piyush Mishra
- Department of Chemistry, Purdue University, 560 Oval Drive, West Lafayette, Indiana 47907, USA
| | - Ha Vinh Lam Nguyen
- Laboratoire Interuniversitaire des Systèmes Atmosphériques (LISA), CNRS UMR 7583, Université Paris-Est Créteil, Université de Paris, Institute Pierre Simon Laplace, 61 avenue du Général de Gaulle, Créteil, France
| | - Timothy S. Zwier
- Department of Chemistry, Purdue University, 560 Oval Drive, West Lafayette, Indiana 47907, USA
| |
Collapse
|
6
|
Herbers S, Fritz SM, Mishra P, Kim Y, Slipchenko L, Zwier TS. The unusual symmetry of hexafluoro-o-xylene-A microwave spectroscopy and computational study. J Chem Phys 2020; 152:064302. [PMID: 32061218 DOI: 10.1063/1.5142169] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The rotational constants and quartic centrifugal distortion coefficients of hexafluoro-o-xylene were precisely derived from the 8 GHz to 18 GHz gas phase microwave spectrum. In addition, the rotational constants of all singly substituted 13C isotopologues were determined. Instead of the intuitively expected symmetry of C2v, as in o-xylene, calculations with a variety of methods (B3LYP, CAM-B3LYP, ωB97XD, MP2, and coupled-cluster singles, doubles, and perturbative triples) predict a C2 symmetry structure in which the two CF3 groups rotate in opposite directions by about 16°. While the experimental results in this study are not capable of proving this unusual symmetry, they can support future microwave, gas phase electron diffraction or nuclear magnetic resonance studies. From the presented data, a preliminary r0 structure was determined, reproducing the experimental rotational constants with deviations of no more than 15 kHz. Analysis of the interactions between the two CF3 groups using an effective fragment potential approach identifies two major contributions to their interaction, due to exchange-repulsion and electrostatic repulsion, with electrostatic repulsion responsible for the barrier at the C2v geometry.
Collapse
Affiliation(s)
- Sven Herbers
- Department of Chemistry, Purdue University, 560 Oval Drive, West Lafayette, Indiana 47907, USA
| | - Sean M Fritz
- Department of Chemistry, Purdue University, 560 Oval Drive, West Lafayette, Indiana 47907, USA
| | - Piyush Mishra
- Department of Chemistry, Purdue University, 560 Oval Drive, West Lafayette, Indiana 47907, USA
| | - Yongbin Kim
- Department of Chemistry, Purdue University, 560 Oval Drive, West Lafayette, Indiana 47907, USA
| | - Lyudmila Slipchenko
- Department of Chemistry, Purdue University, 560 Oval Drive, West Lafayette, Indiana 47907, USA
| | - Timothy S Zwier
- Department of Chemistry, Purdue University, 560 Oval Drive, West Lafayette, Indiana 47907, USA
| |
Collapse
|
7
|
Mishra P, Fritz SM, Hays BM, Mehta-Hurt DN, Jawad KM, Zwier TS. Broadband rotational spectroscopy of trans 3-pentenenitrile and 4-pentenenitrile. Phys Chem Chem Phys 2019; 21:23651-23662. [PMID: 31625538 DOI: 10.1039/c9cp04328d] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022]
Abstract
Titan, a moon of Saturn, has a nitrogen- and methane-rich atmosphere that is similar to prebiotic earth, and is replete with organic nitriles. Pentenenitriles have not yet been detected in Titan's atmosphere or in molecular clouds, but are potential precursors to hetero-aromatic compounds such as pyridine. We performed broadband microwave studies in the 8-18 GHz range on the trans isomer of 3-pentenenitrile (3-PN) and 4-pentenenitrile (4-PN) under jet-cooled conditions. Strong-field coherence breaking (SFCB) was used to selectively modulate the intensities of microwave transitions in a conformer-specific manner for 3-PN, aiding analysis. Two conformers of 3-PN and five conformers of 4-PN were identified and the rotational transitions were assigned. Evidence for methyl internal rotation splitting was observed for both the conformers of 3-PN, and the barrier heights of both conformers was determined experimentally. Comparison is made of the conformational preferences, stability and isomerization barriers through the acquired rotational spectra and potential energy surface (PES) calculations of the structural isomers 3-PN and 4-PN.
Collapse
Affiliation(s)
- Piyush Mishra
- Department of Chemistry, Purdue University, 560 Oval Drive, West Lafayette, IN 47907, USA.
| | | | | | | | | | | |
Collapse
|
8
|
Hernandez-Castillo AO, Abeysekera C, Robicheaux F, Zwier TS. Propagating molecular rotational coherences through single-frequency pulses in the strong field regime. J Chem Phys 2019; 151:084312. [PMID: 31470710 DOI: 10.1063/1.5099049] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
In the weak-field limit in which microwave spectroscopy is typically carried out, an application of a single-frequency pulse that is resonant with a molecular transition will create a coherence between the pair of states involved in the rotational transition, producing a free-induction decay (FID) that, after Fourier transform, produces a molecular signal at that same resonance frequency. With the advent of chirped-pulse Fourier transform microwave methods, the high-powered amplifiers needed to produce broadband microwave spectra also open up other experiments that probe the molecular response in the high-field regime. This paper describes a series of experiments involving resonant frequency pulses interrogating jet-cooled molecules under conditions of sufficient power to Rabi oscillate the two-state system through many Rabi cycles. The Fourier-transformed FID shows coherent signal not only at the applied resonant frequency but also at a series of transitions initially connected to the original one by sharing an upper or lower level with it. As the duration of the single-frequency excitation is increased from 250 to 1500 ns, the number of observed off-resonant, but dipole-allowed, molecular coherences grow. The phenomenon is quite general, having been demonstrated in Z-phenylvinylnitrile, E-phenylvinylnitrile (E-PVN), benzonitrile, guaiacol, and 4-pentynenitrile. In E-PVN, the highest power/longest pulse duration, coherent signal is also present at energetically nearby but not directly connected transitions. Even in molecular samples containing more than one independent species, only transitions due to the single species responsible for the original resonant transition are present. We develop a time-dependent model of the molecular/photon system and use it in conjunction with the experiment to test possible sources of the phenomenon.
Collapse
Affiliation(s)
| | - Chamara Abeysekera
- Department of Chemistry, Purdue University, West Lafayette, Indiana 47907-2084, USA
| | - F Robicheaux
- Department of Physics and Astronomy, Purdue University, West Lafayette, Indiana 47907-2084, USA
| | - Timothy S Zwier
- Department of Chemistry, Purdue University, West Lafayette, Indiana 47907-2084, USA
| |
Collapse
|
9
|
Fritz SM, Mishra P, Zwier TS. Strong-field coherence breaking as a tool for identifying methyl rotor states in microwave spectra: 2-hexanone. J Chem Phys 2019; 151:041104. [PMID: 31370529 DOI: 10.1063/1.5109656] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Sean M. Fritz
- Department of Chemistry, Purdue University, 560 Oval Drive, West Lafayette, Indiana 47907-2084, USA
| | - Piyush Mishra
- Department of Chemistry, Purdue University, 560 Oval Drive, West Lafayette, Indiana 47907-2084, USA
| | - Timothy S. Zwier
- Department of Chemistry, Purdue University, 560 Oval Drive, West Lafayette, Indiana 47907-2084, USA
| |
Collapse
|
10
|
Juanes M, Saragi RT, Caminati W, Lesarri A. The Hydrogen Bond and Beyond: Perspectives for Rotational Investigations of Non‐Covalent Interactions. Chemistry 2019; 25:11402-11411. [DOI: 10.1002/chem.201901113] [Citation(s) in RCA: 63] [Impact Index Per Article: 12.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/10/2019] [Revised: 04/15/2019] [Indexed: 12/22/2022]
Affiliation(s)
- Marcos Juanes
- Departamento de Química Física y Química Inorgánica—IU CINQUIMAFacultad de CienciasUniversidad de Valladolid 47011 Valladolid Spain
| | - Rizalina T. Saragi
- Departamento de Química Física y Química Inorgánica—IU CINQUIMAFacultad de CienciasUniversidad de Valladolid 47011 Valladolid Spain
| | - Walther Caminati
- Dipartimento di Chimica “G. Ciamician”Università di Bologna Via Selmi 2 40126 Bologna Italy
| | - Alberto Lesarri
- Departamento de Química Física y Química Inorgánica—IU CINQUIMAFacultad de CienciasUniversidad de Valladolid 47011 Valladolid Spain
| |
Collapse
|
11
|
Hernandez-Castillo AO, Abeysekera C, Stanton JF, Zwier TS. Structural Characterization of Phenoxy Radical with Mass-Correlated Broadband Microwave Spectroscopy. J Phys Chem Lett 2019; 10:2919-2923. [PMID: 31084013 DOI: 10.1021/acs.jpclett.9b00837] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/09/2023]
Abstract
A combination of broadband microwave spectroscopy and VUV photoionization time-of-flight mass spectra has been used to record rotational spectra of the prototypical phenoxy radical, its per-deuterated isotopomers, and the full set of singly 13C-substituted analogues. Rotational parameters associated with the fits to the full set of isotopomers produce a highly accurate r0 structure for the phenoxy radical. High-level ab initio calculations accurately reproduce the rotational constants and spin-rotation parameters. The structure of the phenoxy radical is distinctly quinoidal, with delocalization of the unpaired electron spin density on the oxygen and phenyl ring. The fitted Fermi contact terms for the 13C atoms reflect a weighting of resonance structures that is 27% on the O atom, 21.5% on each of the two ortho C's, and 30% on the para C, providing a quantitative measure of its sites for subsequent reactions that will control its abundances in combustion and atmospheric environments.
Collapse
Affiliation(s)
- A O Hernandez-Castillo
- Department of Chemistry , Purdue University , 560 Oval Drive , West Lafayette , Indiana 47907-2084 , United States
| | - Chamara Abeysekera
- Department of Chemistry , Purdue University , 560 Oval Drive , West Lafayette , Indiana 47907-2084 , United States
| | - John F Stanton
- Department of Chemistry , University of Florida , P.O. Box 117200, Gainesville , Florida 32611 , United States
| | - Timothy S Zwier
- Department of Chemistry , Purdue University , 560 Oval Drive , West Lafayette , Indiana 47907-2084 , United States
| |
Collapse
|
12
|
Blackshaw KJ, Ortega BI, Quartey NK, Fritzeen WE, Korb RT, Ajmani AK, Montgomery L, Marracci M, Vanegas GG, Galvan J, Sarvas Z, Petit AS, Kidwell NM. Nonstatistical Dissociation Dynamics of Nitroaromatic Chromophores. J Phys Chem A 2019; 123:4262-4273. [DOI: 10.1021/acs.jpca.9b02312] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- K. Jacob Blackshaw
- Department of Chemistry, The College of William and Mary, Williamsburg, Virginia 23187-8795, United States
| | - Belinda I. Ortega
- Department of Chemistry and Biochemistry, California State University—Fullerton, Fullerton, California 92834-6866, United States
| | - Naa-Kwarley Quartey
- Department of Chemistry, The College of William and Mary, Williamsburg, Virginia 23187-8795, United States
| | - Wade E. Fritzeen
- Department of Chemistry, The College of William and Mary, Williamsburg, Virginia 23187-8795, United States
| | - Robert T. Korb
- Department of Chemistry, The College of William and Mary, Williamsburg, Virginia 23187-8795, United States
| | - Annalise K. Ajmani
- Department of Chemistry, The College of William and Mary, Williamsburg, Virginia 23187-8795, United States
| | - Lehman Montgomery
- Department of Chemistry, The College of William and Mary, Williamsburg, Virginia 23187-8795, United States
| | - Marcus Marracci
- Department of Chemistry and Biochemistry, California State University—Fullerton, Fullerton, California 92834-6866, United States
| | - Geronimo Gudino Vanegas
- Department of Chemistry and Biochemistry, California State University—Fullerton, Fullerton, California 92834-6866, United States
| | - John Galvan
- Department of Chemistry and Biochemistry, California State University—Fullerton, Fullerton, California 92834-6866, United States
| | - Zach Sarvas
- Department of Chemistry and Biochemistry, California State University—Fullerton, Fullerton, California 92834-6866, United States
| | - Andrew S. Petit
- Department of Chemistry and Biochemistry, California State University—Fullerton, Fullerton, California 92834-6866, United States
| | - Nathanael M. Kidwell
- Department of Chemistry, The College of William and Mary, Williamsburg, Virginia 23187-8795, United States
| |
Collapse
|