Metal Oxide Nanomaterials in Nanomedicine: Applications in Photodynamic Therapy and Potential Toxicity.

Metal oxide nanomaterials have exhibited excellent performance as nanomedicines in photodynamic therapy (PDT) for cancer and infection treatment. Their unique and tunable physicochemical properties advance them as promising alternatives in drug delivery, early diagnosis, imaging, and treatment against various tumors and infectious diseases. Moreover, the implementation of nanophototherapy in deep tissue sites is enhanced by advancements in photosensitization technology. Notwithstanding the progress made in emerging metal oxide nanomaterials-derived PDT, the potential toxicity towards adjunct tissues associated with this approach remains challenging. Regulation and legislation have also been recommended and subsequently enacted in response to public concerns related to large-scale production, transportation, use, and disposal of those nanomaterials. Consequently, a quantitative structure-activity relationship (QSAR) paradigm has been adopted and is widely used in evaluating and predicting the side effects of nanomedicines, thus influencing their design and fabrication. This article briefly reviews the application of metal oxide nanomaterials in PDT and their associated adverse impacts as reported in recent publications. The future trends and implications of this platform in nanomedicine are also highlighted. However, more studies and efforts have to be carried out for developing novel nano-therapeutics with high selectivity, sensitivity, biocompatibility, and minimal side effects in PDT.

Using experimental data of Escherichia coli to develop a QSAR model for predicting the photo-induced cytotoxicity of metal oxide nanoparticles.

A quantitative structure-activity relationship (QSAR) study of seventeen metal oxide nanoparticles (MNPs), in regard to their photo-induced toxicity to bacteria Escherichia coli, was developed by using quantum chemical methods. A simple and statistically significant QSAR model (F=33.83, R(2)=0.87) was successfully developed for the dark group based on two descriptors, absolute electronegativity of the metal and the metal oxide. Similarly, a best correlation (F=20.51, R(2)=0.804) was obtained to predict the photo-induced toxicity of MNPs by using two descriptors, molar heat capacity and average of the alpha and beta LUMO (lowest unoccupied molecular orbital) energies of the metal oxide. Revelation of these influential molecular descriptors may be useful in elucidating the mechanisms of nanotoxicity and for predicting the environmental risk associated with release of the MNPs. In addition, the developed model may have a role in the future design and manufacture of safe nanomaterials.

Binding of O2 and NO to heme in heme-nitric oxide/oxygen-binding (H-NOX) proteins. A theoretical study.

The binding of O2 and NO to heme in heme-nitric oxide/oxygen-binding (H-NOX) proteins has been investigated with DFT as well as dispersion-corrected DFT methods. The local protein environment was accounted for by including the six nearest surrounding residues in the studied systems. Attention was also paid to the effects of the protein environment, particularly the distal Tyr140, on the proximal iron-histidine (Fe-His) binding. The Heme-AB (AB = O2, NO) and Fe-His binding energies in iron porphyrin FeP(His)(AB), myoglobin Mb(AB), H-NOX(AB), and Tyr140 → Phe mutated H-NOX[Y140F(AB)] were determined for comparison. The calculated stabilization of bound O2 is even higher in H-NOX than that in a myoglobin (Mb), consistent with the observation that the H-NOX domain of T. tengcongensis has a very high affinity for its oxygen molecule. Among the two different X-ray crystal structures for the Tt H-NOX protein, the calculated results for both AB = O2 and NO appear to support the crystal structure with the PDB code 1XBN , where the Trp9 and Asn74 residues do not form a hydrogen-bonding network with Tyr140. A hydrogen bond interaction from the polar residue does not have obvious effects on the Fe-His binding strength, but a dispersion contribution to Ebind(Fe-His) may be significant, depending on the crystal structure used. We speculate that the Fe-His binding strength in the deoxy form of a native protein could be an important factor in determining whether the bond of His to Fe is broken or maintained upon binding of NO to Fe.

Effects of local protein environment on the binding of diatomic molecules to heme in myoglobins. DFT and dispersion-corrected DFT studies.

The heme-AB binding energies (AB = CO, O2) in a wild-type myoglobin (Mb) and two mutants (H64L, V68N) of Mb have been investigated in detail with both DFT and dispersion-corrected DFT methods, where H64L and V68N represent two different, opposite situations. Several dispersion correction approaches were tested in the calculations. The effects of the local protein environment were accounted for by including the five nearest surrounding residues in the calculated systems. The specific role of histidine-64 in the distal pocket was examined in more detail in this study than in other studies in the literature. Although the present calculated results do not change the previous conclusion that the hydrogen bonding by the distal histidine-64 residue plays a major role in the O2/CO discrimination by Mb, more details about the interaction between the protein environment and the bound ligand have been revealed in this study by comparing the binding energies of AB to a porphyrin and the various myoglobins. The changes in the experimental binding energies from one system to another are well reproduced by the calculations. Without constraints on the residues in geometry optimization, the dispersion correction is necessary, since it improves the calculated structures and energetic results significantly.

Factors that distort the heme structure in Heme-Nitric Oxide/OXygen-Binding (H-NOX) protein domains. A theoretical study.

DFT and dispersion-corrected DFT calculations were carried out to probe the factors that distort the heme structure in Heme-Nitric oxide/OXygen-binding (H-NOX) protein domains. Various model systems that include heme, heme+surrounding residues, and heme+surrounding residues+additional protein environment were examined; the latter system was calculated with a quantum mechanics/molecular mechanics (QM/MM) method. The computations were extended to a myoglobin (Mb) protein, in which the heme structure is quite planar, in contrast to that in H-NOX. The natural tendency of the heme is to be planar. The strong structural distortion in H-NOX is mainly brought about by the intermolecular interactions between the whole heme molecule (heme ring plus its peripheral substituents) and the surrounding residues, among which the polar residues (Tyr140, Pro115, Mse98) play major roles in distorting the heme structure. The two peripheral propionate substituents that are oriented on the same side of the heme plane can also make the molecule distort, but the distortion caused by this factor is not significant. In Mb, the surrounding residues considered are all nonpolar and do not cause a structural distortion. The different structural features of the heme macrocycle in the different proteins (H-NOX and Mb) are reproduced by the calculations. The dispersion correction is necessary, since it improves the calculated structures. The effects of the distortion on the binding affinity of the axial ligand to the heme were also examined.

Theoretical study of triatomic silver (Ag3) and its ions with coupled-cluster methods and correlation-consistent basis sets.

Coupled-cluster calculations including non-iterative effects of triple excitations (CCSD(T)) have been made with correlation-consistent basis sets to study a range of properties of Ag3, Ag3(-), and Ag3(+). The methodology was tested on atomic and diatomic silver systems. The accuracy achieved for these systems suggests that predictive-quality results can be expected for the triatomic systems. Properties of the triatomic systems studied include structures and energies of ground and excited states (bent and linear geometries), dissociation energies, ionization energies, and vertical electron detachment energies. In the absence of experimental structural data, the present calculated data serve as reliable predictions. The calculated dissociation energies are near the middle of the experimental range (253 ± 13 kJ mol(-1)). The calculated ionization energies clearly favor the lower of the two experimental measurements (i.e. 5.66 eV). Our calculated vertical electron detachment energies of linear Ag3(-) match the observed photoelectron data fairly well. Overall the higher levels of theory used in the present study achieve consistent accuracy for a range of properties. In contrast, some prior density-functional theory studies provide good results for certain properties, but are lacking in accuracy for others. Another feature of the present study is that it has not been necessary to empirically correct results on the triatomic systems because of known deficiencies in the results on atomic and diatomic systems. The majority of prior wave function-based studies cannot make this claim.

Theoretical investigation into the structural, thermochemical, and electronic properties of the decathio[10]circulene.

For the first time, a theoretical study has been performed on the prototypical decathio[10]circulene (C(20)S(10)) species, which is an analogue of the novel octathio[8]circulene "Sulflower" molecule (C(16)S(8)). Examinations of the singlet and triplet states of C(20)S(10) were made at the B3LYP/6-311G(d) level. Local minima of C(2) and C(s) symmetry were found for the lowest singlet and triplet states, respectively. The stability of C(20)S(10) was assessed by calculating the ΔH°(f) of C(16)S(8) and C(20)S(10) and the ΔH(o) for their decomposition into C(2)S units. Frontier molecular orbital plots show that structural adjacent steric factors along with the twist and strain orientations of C(20)S(10) do not disturb the aromatic π-delocalizing effects. In fact, C(20)S(10) maintains the same p(z) HOMO character as C(16)S(8). These similarities are further verified by density-of-states characterization. Calculated infrared spectra of C(16)S(8) and C(20)S(10) show broad similarities. Molecular electrostatic potential results reveal that eight of the peripheral sulfur atoms are the most electronegative atoms in the molecule, while the interior ten-membered ring exhibits virtually no electronegativity.

FeP(Im)-AB Bonding Energies Evaluated with A Large Number of Density Functionals (P = porphine, Im = imidazole, AB = CO, NO, and O(2)).

Sixty-four (64) density functionals, ranging from GGA, meta-GGA, hybrid GGA to hybrid meta-GGA, were tested to evaluate the FeP(Im)-AB bonding energies (E(bond)) in the heme model complexes FeP(Im)(AB) (P = porphine, Im = imidazole, AB = CO, NO, and O(2)). The results indicate that an accurate prediction of E(bond) for the various ligands to heme is difficult with the DFT methods; usually a functional successful for one system does not perform equally well for the other system(s). Relatively satisfactory results for the various FeP(Im)-AB bonding energies are obtained with the meta-GGA funtionals BLAP3 and Bmτ1; they yield E(bond) values of ca.1.1, 1.2, and 0.4 eV for AB = CO, NO, and O(2), respectively, which are in reasonable agreement with experimental data (0.78 - 0.85 eV for CO, 0.99 eV for NO, and 0.44 - 0.53 eV for O(2)). The other functionals show more or less deficiency for one or two of the systems. The performances of the various functionals in describing the spin-state energetics of the five-coordinate FeP(Im) complex were also examined.

Structure, bonding, and linear optical properties of a series of silver and gold nanorod clusters: DFT/TDDFT studies.

DFT/TDDFT calculations have been carried out for a series of silver and gold nanorod clusters (Ag(n), Au(n), n = 12-120) whose structures are of cigar-type. Pentagonal Ag(n) clusters with n = 49-121 and hexagonal Au(n) clusters with n = 14-74 were also calculated for comparison. Metal-metal distances, binding energies per atom, ionization potentials, and electron affinities were determined, and their trends with cluster size were examined. The TDDFT calculated excitation energies and oscillator strengths were fit by a Lorentz line shape modification, which gives rise to the simulated absorption spectra. The significant features of the experimental spectra for actual silver and gold nanorod particles are well reproduced by the calculations on the clusters. The calculated spectral patterns are also in agreement with previous theoretical results on different-type Ag(n) clusters. Many differences in the calculated properties are found between the Ag(n) and Au(n) clusters, which can be explained by relativistic effects.

Iron porphyrins with different imidazole ligands. A theoretical comparative study.

A theoretical comparative study of a series of five- and six-coordinate iron porphyrins, FeP(L) and FeP(L)(O(2)), has been carried out using DFT methods, where P = porphine and L = imidazole (Im), 1-methylimidazole (1-MeIm), 2-methylimidazole (2-MeIm), 1,2-dimethylimidazole (1,2-Me(2)Im), 4-ethylimidazole (4-EtIm), or histidine (His). Two ligated "picket-fence" iron porphyrins, FeTpivPP(2-MeIm) and FeTpivPP(2-MeIm)(O(2)), were also included in the study for comparison. A number of density functionals were employed in the computations to obtain reliable results. The performance of functionals and basis set effects were investigated in detail on FeP, FeP(Im), and FeP(Im)(O(2)), for which certain experimental information is available and there are some previous calculations in the literature for comparison. Many subtle distinctions in the effects of the different imidazole ligands on the structures and energetics of the deoxy- and oxy iron porphyrins are revealed. While FeP(2-MeIm) is identified to be high spin (S = 2), the ground state of FeP(1-MeIm) may be an admixture of a high-spin (S = 2) and an intermediate-spin (S = 1) state. The ground state of FeP(L)(O(2)) may be different with different L. A weaker Fe-L bond more likely leads to an open-shell singlet ground state for the oxy complex. The 2-methyl group in 2-MeIm, which increases steric contact between the ligand and the porphyrinato skeleton, weakens the Fe-O(2) bond, and thus iron porphyrins with 2-MeIm mimic T-state (low affinity) hemoglobin. The calculated FeP(2-MeIm)-O(2) bonding energy is comparable to the FeTpivPP(2-MeIm)-O(2) one, in agreement with the fact that the picket-fence iron porphyrin binds O(2) with affinity similar to that of myoglobin but different from the result obtained by the CPMD scheme. Im and 4-EtIm closely resemble His, the biologically axial base, and so future computations on hemoprotein models can be simplified safely by using Im.

Theoretical characterization of the F(2)O(3) molecule by coupled-cluster methods.

Coupled-cluster calculations with extended basis sets that include noniterative connected triple excitations (CCSD(T)) have been used to study the FOOOF isomer of F(2)O(3). Second-order Moller-Plessett perturbation theory (MP2) and density-functional theory (B3LYP functional) calculations have also been performed for comparison. Two local minima of similar energy, namely, conformers of C(2) and C(s) symmetry have been located. Structures, harmonic vibrational frequencies, and standard enthalpies and free energies of formation have been calculated. The calculated bond lengths of F(2)O(3) are more characteristic of those in F(2)O and a "normal" peroxide than the unusual bond lengths in F(2)O(2). Both conformers have equal F-O and O-O bond lengths, contrary to a recent suggestion of an unsymmetrical structure. The harmonic vibrational frequencies can aid possible identification of gaseous F(2)O(3). The calculated Δ(f)H° and Δ(f)G° are 110 and 173 kJ mol(-1), respectively. These values are based on extrapolation of CCSD(T) results with augmented triple- and quadruple-ζ basis sets and are expected to be within chemical accuracy (i.e., 1 kcal mol(-1) or 4 kJ mol(-1)). F(2)O(3) is calculated to be stable to decomposition to either FO + FOO or F(2) + O(3), but unstable to decomposition to its elements, to F(2)O(2) + (1)/(2)O(2), and to F(2)O + O(2).

Supramolecular interactions of fullerenes with (Cl)Fe- and Mn porphyrins. A theoretical study.

The electronic structure and bonding in the noncovalent, supramolecular complexes of fullerenes (C(60), C(70)) with (Cl)Fe- and Mn porphyrins [(Cl)FeP, MnP] were investigated in detail with DFT methods. A dispersion correction was made for the fullerene-porphyrin binding energy through an empirical approach. Several density functionals were employed in the calculations in order to obtain reliable results. Our calculated results differ from those obtained in a previous paper (J. Phys. Chem. A, 2005, 109, 3704). The ground state of (Cl)FeP*C(60) is predicted to be high spin (S = 5/2), in agreement with the experimental results. MnP*C(70) is calculated to have a high-spin (S = 5/2) ground state as well; this is similar to (Cl)FeP*C(60), but at variance with the assignment of a low-spin (S = 1/2) state for this complex. According to the calculations, C(70) in MnP*C(70) does not have sufficient ligand-field strength to cause a high- to low-spin state change in MnP. An additional calculation on a comparable, high-spin (Py)MnP complex gives support for the calculated results on MnP*C(70). More detailed experimental investigations are desirable, which might help to resolve the question of the MnP*C(70) electronic structure. The estimated dispersion energies (E(disp)) in the fullerene-porphyrin systems are rather large, ranging from 0.6 to 1.0 eV. Including E(disp) improves the calculated binding energy considerably.

Dispersion-corrected DFT calculations on C(60)-porphyrin complexes.

The quality of the newly added, empirical dispersion correction in density functional theory (DFT) calculations is examined for several supramolecular complexes of fullerene (C(60)) with free-base and metal porphyrins (Por). The benzene dimer (C(6)H(6))(2), naphthalene dimer (C(10)H(8))(2), and anthracene dimer (C(14)H(10))(2) were also included in the study for comparison. Three density functionals, two damping functions, and two types of basis sets were employed in the computations. The estimated dispersion energies in the fullerene-porphyrin systems are rather large, ranging from 0.5 eV in C(60).ZnP to 1 eV in C(60).H(2)TPP. Any dispersion-corrected DFT (DFT + E(disp)) method is shown to perform well for C(60).H(2)TPP, C(60).ZnTPP, and C(60).ZnP, where the intermolecular distances are relatively large. But large basis sets, e.g. TZP (triple-zeta + one polarization function), are required in order to obtain reliable results with DFT + E(disp). In the case of C(60).FeP, where the intermolecular distance R is short, the DFT + E(disp) calculated R depends on the damping function as well as on the DFT method, and all the DFT + E(disp) calculations lead to significant changes in the relative energies of the various spin states. The quality of the DFT + E(disp) calculated results on C(60).FeP is hard to judge here without detailed experimental data on a C(60).FePor complex. Owing to error cancellation, the pure DFT calculations with a smaller DZP (double-zeta + one polarization function) basis set without any correction are shown to give quite accurate results.

Theoretical study of the tautomerism, structures, and vibrational frequencies of the phosphaalkenes XP=C(CH3)2 (X = H, F, Cl, Br, OH, Ar(F) (Ar(F) = 2,6-(CF3)2C6H3)).

Ab initio theoretical calculations have been used to study the influence of phosphorus substituents, Y, on the tautomerism between the vinylphosphine XP(H)C(CH(3))=CH(2) and the phosphaalkene XP=C(CH(3))(2) (X = H, F, Cl, Br, OH, and Ar(F); Ar(F) = 2,6-(CF(3))(2)C(6)H(3)) and on the acidity of the aforementioned vinylphosphine. The stabilization of the phosphaalkene and the increased acidity of the vinylphosphine by Ar(F) are possible factors in the successful synthesis of certain isolable phosphaalkenes. In this work, the properties of Ar(F) are assessed theoretically. Density functional theory using the B3LYP functional has been used for all substituents. In addition, coupled-cluster singles and doubles with noniterative triples (CCSD(T)) has been used for X = H, F, Cl, Br, and OH. The phosphaalkene is favored over the vinylphosphine for all substituents, with F having the strongest stabilizing effect. Cl, Br, and OH have stronger stabilizing effects than Ar(F). In contrast, the most acidic vinylphosphine is that with Ar(F). To aid in the interpretation and analysis of future experimental work, CCSD(T) calculations were used to provide structures and vibrational frequencies for the series XP=C(CH(3))(2) (X = H, F, Cl, Br). The influence of the substituent on geometries and C=P and X-P stretching frequencies was examined, and comparisons were made with the CH(2)=PX series.

Coupled-cluster study of isomers of H2SO2.

A theoretical study has been made on six isomers of H2SO2 using coupled-cluster singles and doubles with noniterative triple excitations (CCSD(T)). The isomers studied are sulfoxylic acid (S(OH)2; C2 and Cs conformers), sulfinic acid (HS(=O)OH; 2 C1 conformers), dihydrogen sulfone (H2SO2; C2v), sulfhydryl hydroperoxide (HSOOH; C1), thiadioxirane (Cs), and dihydrogen persulfoxide (H2SOO; Cs). Molecular geometries, harmonic vibrational frequencies, and infrared intensities of all species were obtained using the CCSD(T) method and the 6-311++G(2d,2p) basis set. All aforementioned species were found to be local minima, with the exception of thiadioxirane, which has one imaginary frequency. A prior possible infrared observation of sulfinic acid was reassessed on the basis of the present data. In agreement with previous MP2 results, the present CCSD(T) data provide support for at most 4 of the 8 observed frequencies. The CCSD(T) frequencies and intensities should be of assistance in future identification of H2SO2 isomers by vibrational spectroscopy. Relative energies were calculated using the CCSD(T) method and several larger basis sets. As found previously, the lowest energy species is C2 S(OH)2, followed by Cs S(OH)2, HS(=O)OH, H2SO2, HSOOH, thiadioxirane, and H2SOO. Expanding the basis set significantly reduces the relative energies of HS(=O)OH and H2SO2. The CCSD(T) method was used with extended basis sets (up to aug-cc-pV(Q+d)Z) and basis set extrapolation in two reaction schemes to calculate the DeltaH degrees t (25 degrees C) of C2 S(OH)2. The two reaction schemes gave -285.8 and -282.7 kJ mol-1, which are quite close to a prior theoretical estimate (-290 kJ mol-1).

Electronic structure of some substituted iron(II) porphyrins. Are they intermediate or high spin?

The electronic structure of some substituted, four-coordinate iron(II) porphyrins has been investigated with DFT methods. These systems include iron tetraphenylporphine (FeTPP), iron octamethyltetrabenzporphine (FeOTBP), iron tetra(alpha,alpha,alpha,alpha-orthopivalamide)phenylporphine (FeTpivPP, also called "picket fence" porphyrin), halogenated iron porphyrins (FeTPPXn, X=F, Cl; n=20, 28), and iron octaethylporphine (FeOEP). A number of density functionals were used in the calculations. Different from the popular, intermediate-spin FeTPP, the ground states of FeOTBP, FeTPPCl28, and FeTPPF20betaCl8 are predicted to be high spin. The calculated result for FeOTBP is in agreement with the early experimental measurement, thereby changing the previous conclusion drawn from the calculations with only the BP functional (J. Chem. Phys. 2002, 116, 3635). But FeTpivPP might have an intermediate-spin ground state, a conclusion that is different from the "experimental" one. With a notably expanded Fe-N bond length, FeOEP might exist as an admixed-spin (S=1, 2) state. We also calculated the electron affinities (EAs) for the various iron porphyrins and compared them to experiment. On the basis of the calculated trends in the EAs and in the orbital energies, the experimental EAs for FeTpivPP, FeTPPF20, and FeTPPCl28 may be too small by 0.4-0.5 eV.

Interaction of metal porphyrins with fullerene C60: a new insight.

The electronic structure and bonding in the noncovalent, supramolecular complexes of fullerene C60 with a series of first-row transition metal porphines MP (M=Fe, Co, Ni, Cu, Zn) have been re-examined with DFT methods. A dispersion correction was made for the C60-MP binding energy through an empirical method (J. Comput. Chem. 2004, 25, 1463). Several density functionals and two types of basis sets were employed in the calculations. Our calculated results are rather different from those obtained in a recent paper (J. Phys. Chem. A 2005, 109, 3704). The ground state of C60.FeP is predicted to be high spin (S=2); the low-spin (S=0), closed-shell state is even higher in energy than the intermediate-spin (S=1) state. With only one electron in the Co-dz2 orbital, the calculated Co-C60 distance is in fact rather short, about 0.1 A longer than the Fe-C60 distance in high-spin C60.FeP. Double occupation of an M-dz2 orbital in MP prevents close association of any axial ligand, and so the Ni-C60, Cu-C60, and Zn-C60 distances are much longer than the Co-C60 one. The evaluated MP-C60 binding energies (Ebind) are 0.8 eV (18.5 kcal/mol) for M=Fe/Co and 0.5 eV (11.5 kcal/mol) for M=Ni/Cu/Zn (Ebind is about 0.2 eV larger in the case of C60-MTPP). They are believed to be reliable and accurate based on our dispersion-corrected DFT calculations that included the counterpoise (CP) correction. The effects of the C60 contact on the redox properties of MP were also examined.

DFT/TDDFT study of lanthanide(III) mono- and bisporphyrin complexes.

The electronic structure, molecular structure, and electronic spectra of lanthanide(III) mono- and bisporphyrin complexes are investigated using a DFT/TDDFT method. These complexes include YbP(acac), YbP(2), [YbP(2)](+), YbHP(2), and [YbP(2)](-) (where P = porphine and acac = acetylacetonate). To shed some light on the origin of the out-of-plane displacement of Yb in YbP(acac), unligated model systems, namely, planar D(4h) and distorted C(4nu) YbP, were calculated. For comparison, the calculations were also extended to include the C and [Ce(IV)P(2) ](+) systems. Even without an axial ligand, the lanthanide atom lies considerably above the porphyrin plane; the distortion of the YbP molecular structure from a planar D(4h) to the nonplanar C(4nu) symmetry leads to a considerable energy lowering. The axial ligand makes the metal out-of-plane displacement even larger, and it also changes the redox properties of the lanthanide monoporphyrin. The ground-state configurations of YbP(2) and YbHP(2) were determined by considering several possible low-lying states. YbP(2) is confirmed to be a single-hole radical. The special redox properties of the bisporphyrin complexes can well be accounted for by the calculated ionization potentials and electron affinities. The TDDFT results provide a clear description of the UV-vis and near-IR absorption spectra of the various lanthanide porphyrins.

Ground and electronically excited states of methyl hydroperoxide: comparison with hydrogen peroxide.

Equilibrium geometries of the ground states of hydrogen peroxide (H(2)O(2)) and methyl hydroperoxide (CH(3)OOH) have been obtained using quadratic configuration interaction methods with correlation-consistent basis sets. These results are compared with experiments and prior calculations. The dipole moments of the ground states of these two molecules have been calculated. The results illustrate the sensitivity of this quantity to molecular geometry. Several excited states of H(2)O(2) and CH(3)OOH were calculated using the equation-of-motion coupled-cluster singles-and-doubles method. Aside from vertical excitation energies, excited state energies along the O-O, O-H, and C-O dissociation pathways were calculated. The results are expected to be of assistance in resolving discrepancies in the experimental interpretation of the UV absorption spectrum and photodissociation of CH(3)OOH.

Assessment of the performance of density-functional methods for calculations on iron porphyrins and related compounds.

The behaviors of a large number of GGA, meta-GGA, and hybrid-GGA density functionals in describing the spin-state energetics of iron porphyrins and related compounds have been investigated. There is a large variation in performance between the various functionals for the calculations of the high-spin state relative energies. Most GGA and meta-GGA functionals are biased toward lower-spin states and so fail to give the correct ground state for the high-spin systems, for which the meta-GGA functionals show more or less improvement over the GGA ones. The GGA functionals that use the OPTX correction for exchange show remarkably high performance for calculating the high-spin state energetics, but their results for the intermediate-spin states are somewhat questionable. A heavily parameterized GGA functional, HCTH/407, provides results which are in qualitative agreement with the experimental findings for the iron porphyrins [FeP, FeP(Cl), FeP(THF)2], but its relative energies for the high-spin states are probably somewhat too low. The high-spin state relative energies are then even more underestimated by the corresponding meta-GGA functional tau-HCTH. For the hybrid-GGA functionals, the Hartree-Fock (HF)-type (or exact) exchange contribution strongly stabilizes the high-spin states, and so the performance of such functionals is largely dependent upon the amount of the HF exchange admixture in them. The B3LYP, B97, B97-1, and tau-HCTH-hyb functionals are able to provide a satisfactory description of the energetics of all the systems considered.