Direct perturbation theory in terms of energy derivatives: scalar-relativistic treatment up to sixth order

A formulation of sixth-order direct perturbation theory (DPT) to treat relativistic effects in quantum-chemical calculations is presented in the framework of derivative theory. Detailed expressions for DPT6 are given at the Hartree-Fock level in terms of the third derivative of the energy with respect to the relativistic perturbation parameter defined as λ(rel)=c(-2). They were implemented for the computation of scalar-relativistic energy corrections. The convergence of the scalar-relativistic DPT expansion is studied for energies and first-order properties such as dipole moment and electric-field gradient within the series of the hydrogen halides (HX, X = F, Cl, Br, I, and At). Comparison with spin-free Dirac-Coulomb calculations indicates that the DPT series exhibits a smooth and monotonic convergence. The rate of convergence, however, depends on the charge of the involved nuclei and significantly slows down for heavy-element compounds.     

Relativistic effects on nuclear magnetic shielding constants in HX and CH3X (X=Br,I) based on the linear response within the elimination of small component approach

Numerical calculations of relativistic effects on nuclear magnetic shielding constants sigma corresponding to all one-body operators obtained within a formalism developed in previous work were carried out. In this formalism, the elimination of small component scheme is applied to evaluate all quantities entering a four-component RSPT(2) expression of magnetic molecular properties. HX and CH3X (X=Br,I) were taken as model compounds. Calculations were carried out at the Hartree-Fock level for first-order quantities, and at the random-phase approximation (RPA) level for second- and third-order ones. It was found that values of sigma(X) are largely affected by several relativistic corrections not previously considered in the bibliography. sigma Values of the H nucleus are in close agreement with four-component RPA ones. Overall relativistic effects on the shift of sigma(X) from HX to CH3X are smaller than the nonrelativistic shifts.     

Direct perturbation theory in terms of energy derivatives: fourth-order relativistic corrections at the Hartree-Fock level

In this work, the quantum-chemical treatment of relativistic effects by means of direct perturbation theory is extended from its lowest order, DPT2, to the next higher order, DPT4. The required theory is given in terms of energy derivatives with the DPT4 energy correction defined as the corresponding second derivative with respect to the relativistic perturbation parameter λ(rel) = c(2) and c as the speed of light. To facilitate the implementation in standard quantum-chemical program packages, a general formulation of DPT starting from a nonrelativistic Lagrangian is developed, thereby expanding both wave function and operators in terms of λ(rel). The corresponding expressions, which incorporate in an additive manner scalar-relativistic and spin-orbit contributions, are given at the Hartree-Fock level and have been implemented in the CFOUR program package using the available analytic second-derivative techniques. The accuracy of the DPT4 corrections at the HF level is investigated by comparison with rigorous four-component calculations. Scalar-relativistic and spin-orbit contributions are analyzed individually and the importance of the various terms to those corrections is discussed. Furthermore, the basis-set dependence of the computed DPT4 corrections is investigated.     

Systematic DFT study of gas phase and solvated uranyl and neptunyl complexes [AnO2X4]n (An = U, Np; X = F, Cl, OH, n = -2; X = H2O, n = +2)

The six-valent uranyl and neptunyl complexes [An(VI)O2X4]n (An = U, Np; X = F, Cl, OH, n = -2; X = H2O, n = +2) have been studied within the framework of density functional theory. The relative stabilities of the cis and trans isomers, structural properties, charge distribution, and ligand binding energies have been determined using the modified Perdew-Burke-Ernzerhof functional at the all-electron scalar relativistic level. Uranyl and neptunyl complexes with different ligands have been compared in a systematic fashion, demonstrating close similarity of these actinides in oxidation state VI. In addition, the effect of an aqueous solution has been taken into account with the polarizable continuum model COSMO. Computed averaged ligand binding energies permit one to rationalize the observed different stabilities of the title species in aqueous media.     

Dimetalloendofullerene U(2)@C(60) has a U-U multiple bond consisting of sixfold one-electron-two-center bonds

Endohedral metallofullerenes (EMFs) have been extensively studied since their discovery in 1985. Metal-metal bonds, nevertheless, have never been explicitly observed in EMFs synthesized so far. In this contribution, we show by means of all-electron relativistic density functional computations that the dimetalloendofullerene, U(2)@C(60), has an unprecedented U-U multiple bond consisting solely of sixfold ferromagnetically coupled one-electron-two-center bonds with the electronic configuration (5fpi(u))(2)(5fsigma(g))(1)(5fdelta(g))(2)(5fphi(u))(1), which are dominated by the uranium 5f atomic orbitals. This bonding scheme is completely distinct from the metal-metal bonds discovered thus far in the d- and f-block polynuclear metal complexes. This finding initiates a connection of the metal-metal multiple bonding chemistry and the fullerene chemistry.     

Simulated annealing and density functional theory calculations of structural and energetic properties of the ammonium chloride clusters (NH4Cl)n, (NH4+)(NH4Cl)n, and (Cl-)(NH4Cl)n, n = 1-13

Simulated annealing Monte Carlo conformer searches using the "mag-walking" algorithm are employed to locate the global minima of molecular clusters of ammonium chloride of the types (NH(4)Cl)(n), (NH(4)(+))(NH(4)Cl)(n), and (Cl(-))(NH(4)Cl)(n) with n = 1-13. The M06-2X density functional theory method is used to refine and predict the structures, energies, and thermodynamic properties of the neutral, cation, and anion clusters. For selected small clusters, the resulting structures are compared to those obtained from a variety of models and basis sets, including RI-MP2 and B3LYP calculations. M06-2X calculations predict enhanced stability of the (NH(4)(+))(NH(4)Cl)(n) clusters when n = 3, 6, 8, and 13. This prediction corresponds favorably to anomalies previously observed in thermospray mass spectroscopy experiments. The (NH(4)Cl)(n) clusters show alternations in stability between even and odd values of n. Clusters of the type (Cl(-))(NH(4)Cl)(n) display a magic number distribution different from that of the cation clusters, with enhanced stability predicted for n = 2, 6, and 11. None of the observed cluster structures resemble the room-temperature CsCl structure of NH(4)Cl(s), which is consistent with previous work. Numerous clusters have structures reminiscent of the higher-temperature, rock-salt phase of the solid ammonium halides.     

Calculation of EPR g tensors for transition-metal complexes based on multiconfigurational perturbation theory (CASPT2).

The computation of the electronic g tensor by two multireference methods is presented and applied to a selection of molecules including CN, BO, AlO, GaO, InO, ZnH, ZnF, O(2), H(2)O(+), O(3) (-), and H(2)CO(+) (group A) as well as TiF(3), CuCl(4) (2-), Cu(NH(3))(4) (2+), and a series of d(1)-MOX(4) (n-) compounds, with M=V, Cr, Mo, Tc, W, Re and X=F, Cl, Br (group B). Two approaches are considered, namely, one in which spin-orbit coupling and the Zeeman effect are included using second-order perturbation theory and another one in which the Zeeman effect is added through first-order degenerate perturbation theory within the ground-state Kramers doublet. The two methods have been implemented into the MOLCAS quantum chemistry software package. The results obtained for the molecules in group A are in good agreement with experiment and with previously reported calculated g values. The results for the molecules in group B vary. While the g values for the d(1) systems are superior to previous theoretical results, those obtained for the d(9) systems are too large compared to the experimental values.     

Synthesis, characterization, and cytotoxicity of platinum(IV) carbamate complexes

The synthesis, characterization, and cytotoxicity of eight new platinum(IV) complexes having the general formula cis,cis,trans-[Pt(NH(3))(2)Cl(2)(O(2)CNHR)(2)] are reported, where R = tert-butyl (4), cyclopentyl (5), cyclohexyl (6), phenyl (7), p-tolyl (8), p-anisole (9), 4-fluorophenyl (10), or 1-naphthyl (11). These compounds were synthesized by reacting organic isocyanates with the platinum(IV) complex cis,cis,trans-[Pt(NH(3))(2)Cl(2)(OH)(2)]. The electrochemistry of the compounds was investigated by cyclic voltammetry. The aryl carbamate complexes 7-11 exhibit reduction peak potentials near -720 mV vs Ag/AgCl, whereas the alkyl carbamate complexes display reduction peak potentials between -820 and -850 mV vs Ag/AgCl. The cyclic voltammograms of cis,cis,trans-[Pt(NH(3))(2)Cl(2)(O(2)CCH(3))(2)] (1), cis,cis,trans-[Pt(NH(3))(2)Cl(2)(O(2)CCF(3))(2)] (2), and cis-[Pt(NH(3))(2)Cl(4)] (3) were measured for comparison. Density functional theory studies were undertaken to investigate the electronic structures of 1-11 and to determine their adiabatic electron affinities. A linear correlation (R(2) = 0.887) between computed adiabatic electron affinities and measured reduction peak potentials was discovered. The biological activity of 4-11 and, for comparison, cisplatin was evaluated in human lung cancer A549 and normal MRC-5 cells by the MTT assay. The compounds exhibit comparable or slightly better activity than cisplatin against the A549 cells. In MRC-5 cells, all are equally or slightly less cytotoxic than cisplatin, except for 4 and 5, which are more toxic.     

Trans or (Unusual) Cis Geometry in d(2) Octahedral Dioxo Complexes. A DFT Study

Geometry optimization of the cis and the trans isomers of several octahedral dioxo complexes of d(2) electronic configuration are performed using the gradient-corrected density functional theory (B3LYP and, for some key structures, BP86). With only monodentate sigma donor ligands (ReO(2)(NH(3))(4)(+), 7), the usual energy order is found (i.e., the trans isomer is the most stable). Complexes with a chelating bidentate ligand, OsO(2)(OCH(2)CH(2)O)(NH(3))(2) (10) and ReO(2)(HN=CHCH=NH)(NH(3))(2)(+) (11), are used as models for the experimental complexes 5 and 2 in which the arrangement of the O=M=O unit is trans and cis, respectively. Our calculations actually show an inversion of the relative energy of the two isomers in going from 10 to 11: while the trans isomer is found to be the most stable in 10, the unusual cis diamagnetic isomer is favored by about 29 kcal mol(-)(1) in 11. This result is traced to the geometric and electronic properties of the bidentate ligand, in particular an acute bite angle and good pi acceptor character. In complex 14 with a bipyridine chelating ligand (weaker pi acceptor than diaza-1,4-butadiene in 11), this energy difference is, however, reduced to 7.5 kcal mol(-)(1) (partial geometry optimization).     

Theoretical study of the potential energy surface for the interaction of cisplatin and their aquated species with water

In this work, a systematic study of the interaction of neutral cisplatin ([Pt(NH(3))(2)Cl(2)]) and their charged aquated species ([Pt(NH(3))(2)Cl(H(2)O)](+) and [Pt(NH(3))(2)(H(2)O)(2)](2+)) with water was carried out. The potential energy surface (PES) was analyzed by considering 35 spatial orientations for the interacting species. The calculations were performed at various levels of theory including Moller-Plesset fourth order perturbation theory and density functional theory (DFT-B3LYP) using extended basis sets. Lennard-Jones (12-6) plus Coulomb classical potential was also used to assess the repulsion-dispersion and electrostatic contributions. The effect of atomic charges on the interaction energies is discussed using Mulliken, charges from electrostatic potential grid method and natural bond orbital schemes. The outcomes show that the electrostatic term plays a primary role on the calculation of interaction energies, with the absolute values of atomic charges from different approaches significantly affecting the overall interaction. Unusual results were revealed by basis set superposition error calculations for the structures located on the platinum-water PES.     

Crystal structures, lattice potential energies, and thermochemical properties of crystalline compounds (1-C(n)H(2n+1)NH3)2ZnCl4(s) (n = 8, 10, 12, and 13)

As part of our ongoing project involving the study of (1-C(n)H(2n+1)NH(3))(2)MCl(4)(s) (where M is a divalent metal ion and n = 8-18), we have synthesized the compounds (1-C(n)H(2n+1)NH(3))(2)ZnCl(4)(s) (n = 8, 10, 12, and 13), and the details of the structures are reported herein. All of the compounds were crystallized in the monoclinic form with the space group P2(1)/n for (1-C(8)H(17)NH(3))(2)ZnCl(4)(s), P21/c for (1-C(10)H(21)NH(3))(2)ZnCl(4)(s), P2(1)/c for (1-C(12)H(25)NH(3))(2)ZnCl(4)(s), and P2(1)/m for (1-C(13)H(27)NH(3))(2)ZnCl(4)(s). The lattice potential energies and ionic volumes of the cations and the common anion of the title compounds were obtained from crystallographic data. Molar enthalpies of dissolution of the four compounds at various molalities were measured at 298.15 K in the double-distilled water. According to Pitzer's theory, molar enthalpies of dissolution of the title compounds at infinite dilution were obtained. Finally, using the values of molar enthalpies of dissolution at infinite dilution (Δ(s)H(m)(∞)) and other auxiliary thermodynamic data, the enthalpy change of the dissociation of [ZnCl(4)](2-)(g) for the reaction [ZnCl(4)](2-)(g)→ Zn(2+)(g) + 4Cl(-)(g) was obtained, and then the hydration enthalpies of cations were calculated by designing a thermochemical cycle.     

Structure and stability of (TiO2)n, (SiO2)n, and mixed Ti(m)Si(n-m)O(2n) [n = 2-5, m = 1 to (n - 1)] clusters

Structures, energetics, and vibrational spectra are investigated for small pure (TiO(2))(n), (SiO(2))(n), and mixed Ti(m)Si(n-m)O(2n) [n = 2-5, m = 1 to (n - 1)] oxide clusters by density functional theory (DFT). The BP86/ATZP level of theory is employed to obtain constitutional isomers of the oxide clusters. In accordance with previous studies, our calculations show three-dimensional compact structures are preferred for pure (TiO(2))(n) with oxo-stabilized higher hexavalent states, and linear chain structures are favored for pure (SiO(2))(n) with tetravalent states. However, the herein theoretically first reported mixed Ti(m)Si(n-m)O(2n) oxide clusters prefer either three-dimensional compact or linear chain structures depending upon the stoichiometry of the compound. Vibrational analysis of the important modes of some highly stable structures is provided. Coupled-cluster single and double excitation (with triples) [CCSD(T)] computed energy gaps for the TiO(2) dimers compare well with results from previous study. Excitation energies are computed by use of time-dependent (TD) DFT and equation-of-motion coupled-cluster calculations with singles and doubles (EOM-CCSD) for the most stable isomers.     

DFT investigation of the tri(amino)amine N(NH(2))(3)(2+) and the tri(azido)amine N(N(3))(3)(2+) dications and related mixed amino(azido)ammonium ions (N(3))(x)N(NH(2))(4-x)(+) (x = 0-4)(1)

Structures of the tri(amino)amine N(NH(2))(3)(2+) and the tri(azido)amine N(N(3))(3)(2+) dications were calculated at the density functional theory (DFT) B3LYP/6-311+G level. The tri(amino)amine dication (NH(2))(3)N(2+) (1) was found to be highly resonance stabilized with a high kinetic barrier for deprotonation. The structures of diamino(azido)amine dication (NH(2))(2)N(N(3))(2+) (2), amino(diazido)amine dication (NH(2))N(N(3))(2)(2+) (3), and tri(azido)amine dication (N(3))(3)N(2+) (4) were also found to be highly resonance stabilized. The structures and energetics of the related mixed amino(azido)ammonium ions (N(3))(x)N(NH(2))(4-x)(+) (x = 0-4) were also calculated.     

A novel series of materials composed of arrays of vanadium oxide container molecules, [V18O42(X)] (X = H2O,Cl-, Br-): synthesis and characterization of [M2(H2N(CH2)2NH2)5][(M(H2N(CH2)2NH2)2]2V18O42(X)]. 9H2O (M = Zn,Cd)

The synthesis and characterization of the novel systems [Zn(2)(H(2)N(CH(2))(2)NH(2))(5)][(Zn(H(2)N(CH(2))(2)NH(2))(2))(2)V(18)O(42)(H(2)O)].9H(2)O (1), [Cd(2)(H(2)N(CH(2))(2)NH(2))(5)][(Cd(H(2)N(CH(2))(2)NH(2))(2))(2)V(18)O(42)(Br)].9H(2)O (2), and [Zn(2)(H(2)N(CH(2))(2)NH(2))(5)][(Zn(H(2)N(CH(2))(2)NH(2))(2))(2)V(18)O(42)(Cl)].9H(2)O (3) have been described. These materials represent a new class of solids that have been prepared by combining conventional coordination compounds with spherical polyoxovanadate clusters. The isomorphous structures of these hybrid solids consist of two-dimensional arrays of container cluster molecules [V(18)O(42)(X)] (X = H(2)O, Br-, Cl-) interlinked by the transition metal complex moieties [M(H(2)N(CH(2))(2)NH(2))(2)] (M = Zn, Cd). These compounds contain an unprecedented complex cation, [M(2)(H(2)N(CH(2))(2)NH(2))(5)](4+). Crystal data for 1: C(9)H(46)N(9)O(26)V(9)Zn(2), monoclinic space group P2(1)/m (No. 11), a = 12.3723(7) A, b = 20.9837(11) A, c = 15.8379(8) A, beta = 97.3320(10) degrees, Z = 4.     

The first acetimino complexes of Rh(III). 4-imino-2-methylpentan-2-amino-Rh(III) complexes through metal-assisted intramolecular aldol-type condensation of two acetimino ligands

[Rh(Cp)Cl(mu-Cl)](2) (Cp = pentamethylcyclopentadienyl) reacts (i) with [Au(NH=CMe(2))(PPh(3))]ClO(4) (1:2) to give [Rh(Cp)(mu-Cl)(NH=CMe(2))](2)(ClO(4))(2) (1), which in turn reacts with PPh(3) (1:2) to give [Rh(Cp)Cl(NH=CMe(2))(PPh(3))]ClO(4) (2), and (ii) with [Ag(NH=CMe(2))(2)]ClO(4) (1:2 or 1:4) to give [Rh(Cp)Cl(NH=CMe(2))(2)]ClO(4) (3) or [Rh(Cp)(NH=CMe(2))(3)](ClO(4))(2).H(2)O (4.H(2)O), respectively. Complex 3 reacts (i) with XyNC (1:1, Xy = 2,6-dimethylphenyl) to give [Rh(Cp)Cl(NH=CMe(2))(CNXy)]ClO(4) (5), (ii) with Tl(acac) (1:1, acacH = acetylacetone) or with [Au(acac)(PPh(3))] (1:1) to give [Rh(Cp)(acac)(NH=CMe(2))]ClO(4) (6), (iii) with [Ag(NH=CMe(2))(2)]ClO(4) (1:1) to give 4, and (iv) with (PPN)Cl (1:1, PPN = Ph(3)P=N=PPh(3)) to give [Rh(Cp)Cl(imam)]Cl (7.Cl), which contains the imam ligand (N,N-NH=C(Me)CH(2)C(Me)(2)NH(2) = 4-imino-2-methylpentan-2-amino) that results from the intramolecular aldol-type condensation of the two acetimino ligands. The homologous perchlorate salt (7.ClO(4)) can be prepared from 7.Cl and AgClO(4) (1:1), by treating 3 with a catalytic amount of Ph(2)C=NH, in an atmosphere of CO, or by reacting 4with (PPN)Cl (1:1). The reactions of 7.ClO(4) with AgClO(4) and PTo(3) (1:1:1, To = C(6)H(4)Me-4) or XyNC (1:1:1) give [Rh(Cp)(imam)(PTo(3))](ClO(4))(2).H(2)O (8) or [Rh(Cp)(imam)(CNXy)](ClO(4))(2) (9), respectively. The crystal structures of 3 and 7.Cl have been determined.     

Electronic structure of square planar bis(benzene-1,2-dithiolato)metal complexes [M(L)(2)](z) (z = 2-, 1-, 0; M = Ni, Pd, Pt, Cu, Au): an experimental, density functional, and correlated ab initio study

The three diamagnetic square planar complexes of nickel(II), palladium(II), and platinum(II) containing two S,S-coordinated 3,5-di-tert-butylbenzene-1,2-dithiolate ligands, (L(Bu))(2-), namely [M(II)(L(Bu))(2)](2-), have been synthesized. The corresponding paramagnetic monoanions [M(II)(L(Bu))(L(Bu)(*))](-) (S = (1)/(2)) and the neutral diamagnetic species [M(II)(L(Bu)(*))(2)] (M = Ni, Pd, Pt) have also been generated in solution or in the solid state as [N(n-Bu)(4)][M(II)(L(Bu))(L(Bu)(*))] salts. The corresponding complex [Cu(III)(L(Bu))(2)](-) has also been investigated. The complexes have been studied by UV-vis, IR, and EPR spectroscopy and by X-ray crystallography; their electro- and magnetochemistry is reported. The electron-transfer series [M(L(Bu))(2)](2-,-,0) is shown to be ligand based involving formally one (L(Bu)(*))(-) pi radical in the monoanion or two in the neutral species [M(II)(L(Bu)(*))(2)] (M = Ni, Pd, Pt). Geometry optimizations using all-electron density functional theory with scalar relativistic corrections at the second-order Douglas-Kroll-Hess (DKH2) and zeroth-order regular approximation (ZORA) levels result in excellent agreement with the experimentally determined structures and electronic spectra. For the three neutral species a detailed analysis of the orbital structures reveals that the species may best be described as containing two strongly antiferromagnetically interacting ligand radicals. Furthermore, multiconfigurational ab initio calculations using the spectroscopy oriented configuration interaction (SORCI) approach including the ZORA correction were carried out. The calculations predict the position of the intervalence charge-transfer band well. Chemical trends in the diradical characters deduced from the multiconfigurational singlet ground-state wave function along a series of metals and ligands were discussed.     

Cl(3)V(mu-S(CH(3))(2))(3)VCl(3)(2)(-): a first-row,face-shared bioctahedral complex with multiple metal-metal bonding

Density functional theory calculations have been used to investigate the structure and bonding of the d(3)d(3) bioctahedral complexes X(3)V(mu-S(CH(3))(2))(3)VX(3)(2)(-) (X = F(-), Cl(-), OH(-), SH(-), NH(2)(-)). According to geometry optimizations using the broken-symmetry approach and the VWN+B-LYP combination of density functionals, the halide-terminated complexes have a V-V bond order of approximately 2, while complexes featuring OH(-), SH(-), or NH(2)(-) as terminal ligands exhibit full triple bonding between the vanadium atoms. The tendency toward triple bonding in the latter complexes is consistent with an increased covalency of the vanadium-ligand bonds, and the influence of bond covalency is apparent also in the tendency for V-V bond elongation in the complexes with OH(-) and NH(2)(-) terminal ligands. Detailed examination of the composition of molecular orbitals in all of the thioether-bridged V(II) complexes substantiates the conclusion that the strong antiferromagnetic coupling which we have determined for these complexes (-J > 250 cm(-)(1)) is due to direct bonding between metal atoms rather than superexchange through the bridging ligands. As such, these V(II) complexes comprise the first apparent examples of multiple metal-metal bonding in first-transition-row, face-shared dinuclear complexes and are therefore of considerable structural and synthetic interest.     

Fourth-order relativistic corrections to electrical first-order properties using direct perturbation theory

In this work, we present relativistic corrections to first-order electrical properties obtained using fourth-order direct perturbation theory (DPT4) at the Hartree-Fock level. The considered properties, i.e., dipole moments and electrical-field gradients, have been calculated using numerical differentiation techniques based on a recently reported DPT4 code for energies [S. Stopkowicz and J. Gauss, J. Chem. Phys. 134, 064114 (2011)]. For the hydrogen halides HX, X=F, Cl, Br, I, and At, we study the convergence of the scalar-relativistic contributions by comparing the computed DPT corrections to results from spin-free Dirac-Hartree-Fock calculations. Furthermore, since in the DPT series spin-orbit contributions first appear at fourth order, we investigate their magnitude and judge the performance of the DPT4 treatment by means of Dirac-Hartree-Fock benchmark calculations. Finally, motivated by experimental investigations of the molecules CH(2)FBr, CHF(2)Br, and CH(2)FI, we present theoretical results for their halogen quadrupole-coupling tensors and give recommendations concerning the importance of higher-order scalar-relativistic and spin-orbit corrections.     

Structures, energetics, vibrational spectra of NH4+ (H2O)(n=4,6) clusters: Ab initio calculations and first principles molecular dynamics simulations

Important structural isomers of NH(4) (+)(H(2)O)(n=4,6) have been studied by using density functional theory, Moller-Plesset second order perturbation theory, and coupled-cluster theory with single, double, and perturbative triple excitations [CCSD(T)]. The zero-point energy (ZPE) correction to the complete basis set limit of the CCSD(T) binding energies and free energies is necessary to identify the low energy structures for NH(4) (+)(H(2)O)(n=4,6) because otherwise wrong structures could be assigned for the most probable structures. For NH(4) (+)(H(2)O)(6), the cage-type structure, which is more stable than the previously reported open structure before the ZPE correction, turns out to be less stable after the ZPE correction. In first principles Car-Parrinello molecular dynamics simulations around 100 K, the combined power spectrum of three lowest energy isomers of NH(4) (+)(H(2)O)(4) and two lowest energy isomers of NH(4) (+)(H(2)O)(6) explains each experimental IR spectrum.