Evaluation of coupled perturbed and density functional methods of computing the parity-violating energy difference between enantiomers.

We present new coupled-perturbed Hartree-Fock (CPHF) and density functional theory (DFT) computations of the parity-violating energy difference (PVED) between enantiomers for H(2)O(2) and H(2)S(2). Our DFT PVED computations are the first for H(2)S(2) and the first with the new HCTH and OLYP functionals. Like other "second generation" PVED computations, our results are an order of magnitude larger than the original "first generation" uncoupled-perturbed Hartree-Fock computations of Mason and Tranter. We offer an explanation for the dramatically larger size in terms of cancellation of contributions of opposing signs, which also explains the basis set sensitivity of the PVED, and its conformational hypersensitivity (addressed in the following paper). This paper also serves as a review of the different types of "second generation" PVED computations: we set our work in context, comparing our results with those of four other groups, and noting the good agreement between results obtained by very different methods. DFT PVEDs tend to be somewhat inflated compared to the CPHF values, but this is not a problem when only sign and order of magnitude are required. Our results with the new OLYP functional are less inflated than those with other functionals, and OLYP is also more efficient computationally. We therefore conclude that DFT computation offers a promising approach for low-cost extension to larger biosystems, especially polymers. The following two papers extend to terrestrial and extra-terrestrial amino acids respectively, and later work will extend to polymers.

Parity-violating energy shifts of Murchison L-amino acids are consistent with an electroweak origin of meteorite L-enantiomeric excesses.

In 1996, four alpha-methyl amino acids in the Murchison meteorite--L-isovaline, L-alpha-methylnorvaline, L-alpha-methyl-allo-isoleucine and L-alpha-methyl-isoleucine--were found to show significant enantiomeric excesses of the L form, ranging from 2% to 9%. Their deuterium to hydrogen isotope ratios suggest they formed in the pre-solar interstellar gas cloud rather than during a later aqueous processing phase on the asteroid parent body. In this paper we apply the techniques of the preceding two papers to compute the parity-violating energy shifts of these amino acids. We find that, in the gas phase, the PVESs of the neutral L forms of all four Murchison alpha-methyl amino acids are decisively negative, and there is even some correlation between the magnitudes of the L-excesses and the magnitudes of the PVESs--all of which is at least consistent with an electroweak origin of the Murchison enantiomeric excesses.

Electroweak parity-violating energy shifts of amino acids: the "conformation problem".

The preceding paper described our coupled-perturbed Hartree-Fock (CPHF) and density functional theory (DFT) methods of computing the parity-violating energy shift (PVES). This paper addresses the "conformation problem"-the difficulty determining which hand of amino acids in solution is favoured by the weak force due to the difficulty determining the solution conformation. We attempt to resolve this by using the methods of the preceding paper to compute the PVES of solution and gas-phase amino acid structures determined by other groups from high level optimizations that include solvation. We conclude that the conformational hypersensitivity of the PVES still precludes a definite conclusion as to the sign of the PVES of L-alanine in solution, but that there is no problem in the gas phase: the PVES of gas-phase L-alanine is decisively negative. We show that the PVES is very sensitive to certain torsion angles, but is not hypersensitive to bondlengths or bond angles. In determining structures for PVES computations, there is therefore no need for expensive full optimizations: one can just optimize the crucial torsion angles. We present new computations of gas-phase amino acids PVESs, using partial optimizations with small basis sets, and the results agree well with those from higher level techniques. In the following paper we apply these less costly techniques to larger amino acids. The "conformation problem" has led some to dismiss the PVES as the source of life's handedness, but we believe this is premature: we show here that amino acids are a special case because their favoured conformations are almost achiral.