The SMFA program for quantum chemistry calculations on large molecules

SMFA is a general program package for performing quantum chemistry calculations on large molecules, using an energy-based fragmentation approach. The program can calculate electronic energies, energy gradients and second derivatives; perform geometry optimization; find first order saddle points (transition states); perform energy optimized scans along a user-defined path; and evaluate various molecular properties. The program can use any of the following quantum chemistry packages: GAMESS(US), GAUSSIAN, NWChem and Q-Chem. In addition, SMFA provides a number of utility programs that, inter alia, calculate vibrational frequencies and infrared spectra with isotopic substitutions, the electrostatic potential on the solvent-accessible-surface, and isodesmic and higher order near-iso-energetic reaction schemes. Calculations of the electronic energy and related properties can be carried out using a scheme that provides a computation time that is linearly dependent on the size of the molecule or, if the user has enough processing units available, in a walltime that is independent of the size of the molecule

The practical Pomeron for high energy proton collimation

We present a model which describes proton scattering data from ISR to Tevatron energies, and which can be applied to collimation in high energy accelerators, such as the LHC and FCC. Collimators remove beam halo particles, so that they do not impinge on vulnerable regions of the machine, such as the superconducting magnets and the experimental areas. In simulating the effect of the collimator jaws it is crucial to model the scattering of protons at small momentum transfer t, as these protons can subsequently survive several turns of the ring before being lost. At high energies these soft processes are well described by Pomeron exchange models. We study the behaviour of elastic and single-diffractive dissociation cross sections over a wide range of energy, and show that the model can be used as a global description of the wide variety of high energy elastic and diffractive data presently available. In particular it models low mass diffraction dissociation, where a rich resonance structure is present, and thus predicts the differential and integrated cross sections in the kinematical range appropriate to the LHC. We incorporate the physics of this model into the beam tracking code MERLIN and use it to simulate the resulting loss maps of the beam halo lost in the collimators in the LHC

Genetic Variation and De Novo Mutations in the Parthenogenetic Caucasian Rock Lizard Darevskia unisexualis

Unisexual all-female lizards of the genus Darevskia that are well adapted to various habitats are known to reproduce normally by true parthenogenesis. Although they consist of unisexual lineages and lack effective genetic recombination, they are characterized by some level of genetic polymorphism. To reveal the mutational contribution to overall genetic variability, the most straightforward and conclusive way is the direct detection of mutation events in pedigree genotyping. Earlier we selected from genomic library of D. unisexualis two polymorphic microsatellite containg loci Du281 and Du215. In this study, these two loci were analyzed to detect possible de novo mutations in 168 parthenogenetic offspring of 49 D. unisexualis mothers and in 147 offspring of 50 D. armeniaca mothers . No mutant alleles were detected in D. armeniaca offspring at both loci, and in D. unisexualis offspring at the Du215 locus. There were a total of seven mutational events in the germ lines of four of the 49 D. unisexualis mothers at the Du281 locus, yielding the mutation rate of 0.1428 events per germ line tissue. Sequencing of the mutant alleles has shown that most mutations occur via deletion or insertion of single microsatellite repeat being identical in all offspring of the family. This indicates that such mutations emerge at the early stages of embryogenesis. In this study we characterized single highly unstable (GATA)n containing locus in parthenogenetic lizard species D. unisexualis. Besides, we characterized various types of mutant alleles of this locus found in the D. unisexualis offspring of the first generation. Our data has shown that microsatellite mutations at highly unstable loci can make a significant contribution to population variability of parthenogenetic lizards

Search for New Physics in e mu X Data at D0 Using Sleuth: A Quasi-Model-Independent Search Strategy for New Physics

We present a quasi-model-independent search for the physics responsible for electroweak symmetry breaking. We define final states to be studied, and construct a rule that identifies a set of relevant variables for any particular final state. A new algorithm ("Sleuth") searches for regions of excess in those variables and quantifies the significance of any detected excess. After demonstrating the sensitivity of the method, we apply it to the semi-inclusive channel e mu X collected in 108 pb^-1 of ppbar collisions at sqrt(s) = 1.8 TeV at the D0 experiment during 1992-1996 at the Fermilab Tevatron. We find no evidence of new high p_T physics in this sample.Comment: 23 pages, 12 figures. Submitted to Physical Review

Limits on WWZ and WW\gamma couplings from p\bar{p}\to e\nu jj X events at \sqrt{s} = 1.8 TeV

We present limits on anomalous WWZ and WW-gamma couplings from a search for WW and WZ production in p-bar p collisions at sqrt(s)=1.8 TeV. We use p-bar p -> e-nu jjX events recorded with the D0 detector at the Fermilab Tevatron Collider during the 1992-1995 run. The data sample corresponds to an integrated luminosity of 96.0+-5.1 pb^(-1). Assuming identical WWZ and WW-gamma coupling parameters, the 95% CL limits on the CP-conserving couplings are -0.33<lambda<0.36 (Delta-kappa=0) and -0.43<Delta-kappa<0.59 (lambda=0), for a form factor scale Lambda = 2.0 TeV. Limits based on other assumptions are also presented.Comment: 11 pages, 2 figures, 2 table

Search For Heavy Pointlike Dirac Monopoles

We have searched for central production of a pair of photons with high transverse energies in $p\bar p$ collisions at $\sqrt{s} = 1.8$ TeV using $70 pb^{-1}$ of data collected with the D\O detector at the Fermilab Tevatron in 1994--1996. If they exist, virtual heavy pointlike Dirac monopoles could rescatter pairs of nearly real photons into this final state via a box diagram. We observe no excess of events above background, and set lower 95% C.L. limits of $610, 870, or 1580 GeV/c^2$ on the mass of a spin 0, 1/2, or 1 Dirac monopole.Comment: 12 pages, 4 figure

The Dijet Mass Spectrum and a Search for Quark Compositeness in bar{p}p Collisions at sqrt{s} = 1.8 TeV

Using the DZero detector at the 1.8 TeV pbarp Fermilab Tevatron collider, we have measured the inclusive dijet mass spectrum in the central pseudorapidity region |eta_jet| < 1.0 for dijet masses greater than 200 Gev/c^2. We have also measured the ratio of spectra sigma(|eta_jet| < 0.5)/sigma(0.5 < |eta_jet| < 1.0). The order alpha_s^3 QCD predictions are in good agreement with the data and we rule out models of quark compositeness with a contact interaction scale < 2.4 TeV at the 95% confidence level.Comment: 11 pages, 4 figures, 2 tables, submitted to Phys. Rev. Let

Search for High Mass Photon Pairs in p-pbar --> gamma-gamma-jet-jet Events at sqrt(s)=1.8 TeV

A search has been carried out for events in the channel p-barp --> gamma gamma jet jet. Such a signature can characterize the production of a non-standard Higgs boson together with a W or Z boson. We refer to this non-standard Higgs, having standard model couplings to vector bosons but no coupling to fermions, as a "bosonic Higgs." With the requirement of two high transverse energy photons and two jets, the diphoton mass (m(gamma gamma)) distribution is consistent with expected background. A 90(95)% C.L. upper limit on the cross section as a function of mass is calculated, ranging from 0.60(0.80) pb for m(gamma gamma) = 65 GeV/c^2 to 0.26(0.34) pb for m(gamma gamma) = 150 GeV/c^2, corresponding to a 95% C.L. lower limit on the mass of a bosonic Higgs of 78.5 GeV/c^2.Comment: 9 pages, 3 figures. Replacement has new H->gamma gamma branching ratios and corresponding new mass limit

Ratio of the Isolated Photon Cross Sections at \sqrt{s} = 630 and 1800 GeV

The inclusive cross section for production of isolated photons has been measured in \pbarp collisions at $\sqrt{s} = 630$ GeV with the \D0 detector at the Fermilab Tevatron Collider. The photons span a transverse energy ($E_T$) range from 7-49 GeV and have pseudorapidity $|\eta| < 2.5$. This measurement is combined with to previous \D0 result at $\sqrt{s} = 1800$ GeV to form a ratio of the cross sections. Comparison of next-to-leading order QCD with the measured cross section at 630 GeV and ratio of cross sections show satisfactory agreement in most of the $E_T$ range.Comment: 7 pages. Published in Phys. Rev. Lett. 87, 251805, (2001

Zgamma Production in pbarp Collisions at sqrt(s)=1.8 TeV and Limits on Anomalous ZZgamma and Zgammagamma Couplings

We present a study of Z +gamma + X production in p-bar p collisions at sqrt{S}=1.8 TeV from 97 (87) pb^{-1} of data collected in the eegamma (mumugamma) decay channel with the D0 detector at Fermilab. The event yield and kinematic characteristics are consistent with the Standard Model predictions. We obtain limits on anomalous ZZgamma and Zgammagamma couplings for form factor scales Lambda = 500 GeV and Lambda = 750 GeV. Combining this analysis with our previous results yields 95% CL limits |h{Z}_{30}| < 0.36, |h{Z}_{40}| < 0.05, |h{gamma}_{30}| < 0.37, and |h{gamma}_{40}| < 0.05 for a form factor scale Lambda=750 GeV.Comment: 17 Pages including 2 Figures. Submitted to PR