Force constants and dipole moment derivatives of molecules from perturbed Hartree-Fock calculations II. Applications to limited basis set SCF-MO wavefunctions

The perturbed Hartree–Fock theory developed in the preceding paper is applied to LiH, BH, and HF, using limited basis‐set SCF–MO wavefunctions derived by previous workers. The calculated values for the force constant ke and the dipole‐moment derivative μ(1) are (experimental values in parentheses): LiH, ke  =  1.618(1.026)mdyn/Å,μ(1)  =  −18.77(−2.0±0.3)D/ÅBH,ke  =  5.199(3.032)mdyn/Å,μ(1)  =  −1.03(−)D/Å;HF,ke  =  12.90(9.651)mdyn/Å,μ(1)  =  −2.15(+1.50)D/Å. The values of the force on the proton were calculated exactly and according to the Hellmann–Feynman theorem in each case, and the discrepancies show that none of the wavefunctions used are close to the Hartree–Fock limit, so that the large errors in ke and μ(1) are not surprising. However no difficulties arose in the perturbed Hartree–Fock calculation, so that the application of the theory to more accurate wavefunctions appears quite feasible.

The calculation of force constants and normal co-ordinates II: methyl fluoride

The mathematical difficulties which can arise in the force constant refinement procedure for calculating force constants and normal co-ordinates are described and discussed. The method has been applied to the methyl fluoride molecule, using an electronic computer. The best values of the twelve force constants in the most general harmonic potential field were obtained to fit twenty-two independently observed experimental data, these being the six vibration frequencies, three Coriolis zeta constants and two centrifugal stretching constants DJ and DJK, for both CH3F and CD3F. The calculations have been repeated both with and without anharmonicity corrections to the vibration frequencies. All the experimental data were weighted according to the reliability of the observations, and the corresponding standard errors and correlation coefficients of the force constants have been deduced. The final force constants are discussed briefly, and compared with previous treatments, particularly with a recent Urey-Bradley treatment for this molecule.

The two-dimensional anharmonic oscillator I: vibration-rotation constants of fulminic acid, HCNO

The lowest-wavenumber vibration of HCNO and DCNO, ν5, is known to involve a largeamplitude low-frequency anharmonic bending of the CH bond against the CNO frame. In this paper the anomalous vibrational dependence of the observed rotational constants B(v5, l5), and of the observed l-doubling interactions, is interpreted according to a simple effective vibration-rotation Hamiltonian in which the appropriate vibrational operators are averaged in an anharmonic potential surface over the normal coordinates (Q5x, Q5y). All of the data on both isotopes are interpreted according to a single potential surface having a minimum energy at a slightly bent configuration of the HCN angle ( 170°) with a maximum at the linear configuration about 2 cm−1 higher. The other coefficients in the Hamiltonian are also interpreted in terms of the structure and the harmonic and anharmonic force fields; the substitution structure at the “hypothetical linear configuration” determined in this way gives a CH bond length of 1.060 Å, in contrast to the value 1.027 Å determined from the ground-state rotational constants. We also discuss the difficulties in rationalizing our effective Hamiltonian in terms of more fundamental theory, as well as the success and limitations of its use in practice.

The calculation of force constants and normal coordinates III: constrained force fields for the methyl halides

A method is discussed for imposing any desired constraint on the force field obtained in a force constant refinement calculation. The application of this method to force constant refinement calculations for the methyl halide molecules is reported. All available data on the vibration frequencies, Coriolis interaction constants and centrifugal stretching constants of CH3X and CD3X molecules were used in the refinements, but despite this apparent abundance of data it was found that constraints were necessary in order to obtain a unique solution to the force field. The results of unconstrained calculations, and of three different constrained calculations, are reported in this paper. The constrained models reported are a Urey—Bradley force field, a modified valence force field, and a constraint based on orbital-following bond-hybridization arguments developed in the following paper. The results are discussed, and compared with previous results for these molecules. The third of the above models is found to reproduce the observed data better than either of the first two, and additional reasons are given for preferring this solution to the force field for the methyl halide molecules.

Signs of l-doubling constants in NF3

Computed infrared band contours are presented for the two degenerate fundamentals of the NF3 molecule, using the l-resonance band contour program described by Cartwright and Mills, with values of the vibration-rotation constants determined from the microwave spectrum by Morino and co-workers. Computed contours are presented for both possible signs of the l-doubling constant, and comparison with the observed band contours leads to the conclusion that q3 = −121.4 MHz and q4 = +51.4 MHz.

The calculation of force constants and normal co-ordinates--VII Normal co-ordinates in the E species of allene

Force constant and normal co-ordinate calculations are reported for the E species vibrations of the allene molecule. Data on the fundamental vibration frequencies of allene-h4, allene-d4 and allene-1.1-d2 and on the five experimentally determined Coriolis zeta constants of C3H4 and C3D4, were used in a force constant refinement procedure. Allowing for product and sum rules this gives 21 independent data which were used to refine to the most general harmonic force field (10 parameters) with one constraint (in the absence of any constraints the refinement was not satisfactory). The results have been used to calculate the complete ζz Coriolis interaction matrix for the allene-1.1-d2 molecule, and hence to calculate the expected rotational structure of the perpendicular bending vibrations of this molecule; the good agreement obtained with the observed spectra is a check on our results.

Efficient and cost-effective experimental determination of kinetic constants and data: the success of a Bayesian systematic approach to drug transport, receptor binding, continuous culture and cell transport kinetics

Details about the parameters of kinetic systems are crucial for progress in both medical and industrial research, including drug development, clinical diagnosis and biotechnology applications. Such details must be collected by a series of kinetic experiments and investigations. The correct design of the experiment is essential to collecting data suitable for analysis, modelling and deriving the correct information. We have developed a systematic and iterative Bayesian method and sets of rules for the design of enzyme kinetic experiments. Our method selects the optimum design to collect data suitable for accurate modelling and analysis and minimises the error in the parameters estimated. The rules select features of the design such as the substrate range and the number of measurements. We show here that this method can be directly applied to the study of other important kinetic systems, including drug transport, receptor binding, microbial culture and cell transport kinetics. It is possible to reduce the errors in the estimated parameters and, most importantly, increase the efficiency and cost-effectiveness by reducing the necessary amount of experiments and data points measured. (C) 2003 Federation of European Biochemical Societies. Published by Elsevier B.V. All rights reserved.

Direct time-resolved study of the gas-phase reactions of germylene with ethyl- and diethylgermane: absolute rate constants, temperature dependences, and mechanism

Time-resolved studies of germylene, GeH2, generated by the 193 nm laser flash photolysis of 3,4-dimethyl-1-germacyclopent-3-ene, have been carried out to obtain rate constants for its bimolecular reactions with ethyl- and diethylgermanes in the gas phase. The reactions were studied over the pressure range 1-100 Torr with SF6 as bath gas and at five temperatures in the range 297-564 K. Only slight pressure dependences were found for GeH2 + EtGeH3 (399, 486, and 564 K). The high pressure rate constants gave the following Arrhenius parameters: for GeH2 + EtGeH3, log A = -10.75 +/- 0.08 and E-a = -6.7 +/- 0.6 kJ mol(-1); for GeH2 + Et2GeH2, log A = -10.68 +/- 0.11 and E-a = -6.95 +/- 0.80 kJ mol(-1). These are consistent with fast, near collision-controlled, association processes at 298 K. RRKM modeling calculations are, for the most part, consistent with the observed pressure dependence of GeH2 + EtGeH3. The ethyl substituent effects have been extracted from these results and are much larger than the analogous methyl substituent effects in the SiH2 + methylsilane reaction series. This is consistent with a mechanistic model for Ge-H insertion in which the intermediate complex has a sizable secondary barrier to rearrangement.

Gas phase kinetic and quantum chemical studies of the reactions of silylene with the methylsilanes. Absolute rate constants, temperature dependences, RRKM modelling and potential energy surfaces

Time resolved studies of silylene, SiH2, generated by the 193 nm laser. ash photolysis of phenylsilane, have been carried out to obtain rate coefficients for its bimolecular reactions with methyl-, dimethyl- and trimethyl-silanes in the gas phase. The reactions were studied over the pressure range 3 - 100 Torr with SF6 as bath gas and at five temperatures in the range 300 - 625 K. Only slight pressure dependences were found for SiH2 + MeSiH3 ( 485 and 602 K) and for SiH2 + Me2SiH2 ( 600 K). The high pressure rate constants gave the following Arrhenius parameters: [GRAPHICS] These are consistent with fast, near to collision-controlled, association processes. RRKM modelling calculations are consistent with the observed pressure dependences ( and also the lack of them for SiH2 + Me3SiH). Ab initio calculations at both second order perturbation theory (MP2) and coupled cluster (CCSD(T)) levels, showed the presence of weakly-bound complexes along the reaction pathways. In the case of SiH2 + MeSiH3 two complexes, with different geometries, were obtained consistent with earlier studies of SiH2 + SiH4. These complexes were stabilised by methyl substitution in the substrate silane, but all had exceedingly low barriers to rearrangement to product disilanes. Although methyl groups in the substrate silane enhance the intrinsic SiH2 insertion rates, it is doubtful whether the intermediate complexes have a significant effect on the kinetics. A further calculation on the reaction MeSiH + SiH4 shows that the methyl substitution in the silylene should have a much more significant kinetic effect ( as observed in other studies).

Metal complexes of a dipyridine octaazamacrocycle: stability constants, structural and modelling studies

Two 28-membered octaazamacrocycles, [28]py(2)N(6) and Me-2[28]py(2)N(6), have been synthesized. The protonation constants of the N-methyl. derivative and the stability constants of its complexes with Ni2+, Cu2+, Zn2+, Cd2+, and Pb2+ were determined at 25degreesC in 0.10 mol dm(-3) KNO3. The high overall basicity of Me-2[28]py(2)N(6) is ascribed to the weaker repulsion between protonated contiguous charged ammonium sites separated by propyl chains. These studies together with NMR, UV-vis and EPR spectroscopies indicated the presence of mono- and di-nuclear species, The single crystal structure of the complex [Ni-2([28]py(2)N(6))(H2O)(4)]Cl-4.3H(2)O was determined, and showed each nickel centre in a distorted octahedral co-ordination environment. The nickel centres are held within the macrocycle at a large distance of 6.991(g) Angstrom from each other. The formation of mononuclear complexes was evaluated theoretically via molecular mechanics (MM) and molecular dynamics (MD) calculations and showed that these large macrocycles have sufficient flexibility to encapsulate metal ions with different stereo-electronic sizes. Structures for small and large metal ions are proposed.

Metal complexes of cyclen and cyclam derivatives useful for medical applications: a discussion based on thermodynamic stability constants and structural data

A series of the most common chelators used in magnetic resonance imaging ( MRI) and in radiopharmaceuticals for medical diagnosis and tumour therapy, H(4)dota, H(4)teta, H(8)dotp and H(8)tetp, is examined from a chemical point of view. Differences between 12- and 14-membered tetraazamacrocyclic derivatives with methylcarboxylate and methylphosphonate pendant arms and their chelates with divalent first-series transition metal and trivalent lanthanide ions are discussed on the basis of their thermodynamic stability constants, X- ray structures and theoretical studies.

On defining base units in terms of fundamental constants

The definitions of the base units of the international system of units have been revised many times since the idea of such an international system was first conceived at the time of the French revolution. The objective today is to define all our units in terms of 'invariants of nature', i.e. by referencing our units to the fundamental constants of physics, or the properties of atoms, rather than the characteristics of our planet or of artefacts. This situation is reviewed, particularly in regard to finding a new definition of the kilogram to replace its present definition in terms of a prototype material artefact.

Thiol and Sulfenic Acid Oxidation of AhpE, the One-Cysteine Peroxiredoxin from Mycobacterium tuberculosis: Kinetics, Acidity Constants, and Conformational Dynamics

Drug resistance and virulence of Mycobacterium tuberculosis are partially related to the pathogen`s antioxidant systems. Peroxide detoxification in this bacterium is achieved by the heme-containing catalase peroxidase and different two-cysteine peroxiredoxins. M. tuberculosis genome also codifies for a putative one-cysteine peroxiredoxin, alkyl hydroperoxide reductase E (MtAhpE). Its expression was previously demonstrated at a transcriptional level, and the crystallographic structure of the recombinant protein was resolved under reduced and oxidized states. Herein, we report that the conformation of MtAhpE changed depending on its single cysteine redox state, as reflected by different tryptophan fluorescence properties and changes in quaternary structure. Dynamics of fluorescence changes, complemented by competition kinetic assays, were used to perform protein functional studies. MtAhE reduced peroxynitrite 2 orders of magnitude faster than hydrogen peroxide (1.9 x 10(7) M(-1) s(-1) vs 8.2 x 10(4) M(-1) s(-1) at pH 7.4 and 25 degrees C, respectively). The latter also caused cysteine overoxidation to sulfinic acid, but at much slower rate constant (40 M(-1) s(-1)). The pK(a) of the thiol in the reduced enzyme was 5.2, more than one unit lower than that of the sulfenic acid in the oxidized enzyme. The pH profile of hydrogen peroxide-mediated thiol and sulfenic acid oxidations indicated thiolate and sulfenate as the reacting species. The formation of sulfenic acid as well as the catalytic peroxidase activity of MtAhpE was demonstrated using the artificial reducing substrate thionitrobenzoate. Taken together, our results indicate that MtAhpE is a relevant component in the antioxidant repertoire of M. tuberculosis probably involved in peroxide and specially peroxynitrite detoxification.

NMR Chemical Shielding and Spin-Spin Coupling Constants of Liquid NH(3): A Systematic Investigation using the Sequential QM/MM Method

The NMR spin coupling parameters, (1)J(N,H) and (2)J(H,H), and the chemical shielding, sigma((15)N), of liquid ammonia are studied from a combined and sequential QM/MM methodology. Monte Carlo simulations are performed to generate statistically uncorrelated configurations that are submitted to density functional theory calculations. Two different Lennard-Jones potentials are used in the liquid simulations. Electronic polarization is included in these two potentials via an iterative procedure with and without geometry relaxation, and the influence on the calculated properties are analyzed. B3LYP/aug-cc-pVTZ-J calculations were used to compute the V(N,H) constants in the interval of -67.8 to -63.9 Hz, depending on the theoretical model used. These can be compared with the experimental results of -61.6 Hz. For the (2)J(H,H) coupling the theoretical results vary between -10.6 to -13.01 Hz. The indirect experimental result derived from partially deuterated liquid is -11.1 Hz. Inclusion of explicit hydrogen bonded molecules gives a small but important contribution. The vapor-to-liquid shifts are also considered. This shift is calculated to be negligible for (1)J(N,H) in agreement with experiment. This is rationalized as a cancellation of the geometry relaxation and pure solvent effects. For the chemical shielding, U(15 N) Calculations at the B3LYP/aug-pcS-3 show that the vapor-to-liquid chemical shift requires the explicit use of solvent molecules. Considering only one ammonia molecule in an electrostatic embedding gives a wrong sign for the chemical shift that is corrected only with the use of explicit additional molecules. The best result calculated for the vapor to liquid chemical shift Delta sigma((15)N) is -25.2 ppm, in good agreement with the experimental value of -22.6 ppm.