Geant4 Hadronic Cascade Models and CMS Data Analysis : Computational Challenges in the LHC era

This work belongs to the field of computational high-energy physics (HEP). The key methods used in this thesis work to meet the challenges raised by the Large Hadron Collider (LHC) era experiments are object-orientation with software engineering, Monte Carlo simulation, the computer technology of clusters, and artificial neural networks. The first aspect discussed is the development of hadronic cascade models, used for the accurate simulation of medium-energy hadron-nucleus reactions, up to 10 GeV. These models are typically needed in hadronic calorimeter studies and in the estimation of radiation backgrounds. Various applications outside HEP include the medical field (such as hadron treatment simulations), space science (satellite shielding), and nuclear physics (spallation studies). Validation results are presented for several significant improvements released in Geant4 simulation tool, and the significance of the new models for computing in the Large Hadron Collider era is estimated. In particular, we estimate the ability of the Bertini cascade to simulate Compact Muon Solenoid (CMS) hadron calorimeter HCAL. LHC test beam activity has a tightly coupled cycle of simulation-to-data analysis. Typically, a Geant4 computer experiment is used to understand test beam measurements. Thus an another aspect of this thesis is a description of studies related to developing new CMS H2 test beam data analysis tools and performing data analysis on the basis of CMS Monte Carlo events. These events have been simulated in detail using Geant4 physics models, full CMS detector description, and event reconstruction. Using the ROOT data analysis framework we have developed an offline ANN-based approach to tag b-jets associated with heavy neutral Higgs particles, and we show that this kind of NN methodology can be successfully used to separate the Higgs signal from the background in the CMS experiment.

Visual Analysis of Biomolecular Surfaces

Surface models of biomolecules have become crucially important for the study and understanding of interaction between biomolecules and their environment. We argue for the need for a detailed understanding of biomolecular surfaces by describing several applications in computational and structural biology. We review methods used to model, represent, characterize, and visualize biomolecular surfaces focusing on the role that geometry and topology play in identifying features on the surface. These methods enable the development of efficient computational and visualization tools for studying the function of biomolecules.

Theoretical and computational approaches on heterogeneous nucleation

Nucleation is the first step of a first order phase transition. A new phase is always sprung up in nucleation phenomena. The two main categories of nucleation are homogeneous nucleation, where the new phase is formed in a uniform substance, and heterogeneous nucleation, when nucleation occurs on a pre-existing surface. In this thesis the main attention is paid on heterogeneous nucleation. This thesis wields the nucleation phenomena from two theoretical perspectives: the classical nucleation theory and the statistical mechanical approach. The formulation of the classical nucleation theory relies on equilibrium thermodynamics and use of macroscopically determined quantities to describe the properties of small nuclei, sometimes consisting of just a few molecules. The statistical mechanical approach is based on interactions between single molecules, and does not bear the same assumptions as the classical theory. This work gathers up the present theoretical knowledge of heterogeneous nucleation and utilizes it in computational model studies. A new exact molecular approach on heterogeneous nucleation was introduced and tested by Monte Carlo simulations. The results obtained from the molecular simulations were interpreted by means of the concepts of the classical nucleation theory. Numerical calculations were carried out for a variety of substances nucleating on different substances. The classical theory of heterogeneous nucleation was employed in calculations of one-component nucleation of water on newsprint paper, Teflon and cellulose film, and binary nucleation of water-n-propanol and water-sulphuric acid mixtures on silver nanoparticles. The results were compared with experimental results. The molecular simulation studies involved homogeneous nucleation of argon and heterogeneous nucleation of argon on a planar platinum surface. It was found out that the use of a microscopical contact angle as a fitting parameter in calculations based on the classical theory of heterogeneous nucleation leads to a fair agreement between the theoretical predictions and experimental results. In the presented cases the microscopical angle was found to be always smaller than the contact angle obtained from macroscopical measurements. Furthermore, molecular Monte Carlo simulations revealed that the concept of the geometrical contact parameter in heterogeneous nucleation calculations can work surprisingly well even for very small clusters.

Enhanced Ion Anisotropy by Nonconventional Coordination Geometry: Single-Chain Magnet Behavior for a [{(FeL)-L-II}(2){Nb-IV(CN)(8)}] Helical Chain Compound Designed with Heptacoordinate Fe-II

Nonconventional heptacoordination in combination with efficient magnetic exchange coupling is shown to yield a 1-D heteronuclear {(FeNbIV)-Nb-II} compound with remarkable magnetic features when compared to other Fe(II)-based single chain magnets (SCM). Cyano-bridged heterometallic {3d-4d} and {3d-5d} chains are formed upon assembling Fe(II) bearing a pentadentate macrocycle as the blocking ligand with octacyano metallates, [M(CN)(8)](4-) (M = Nb-IV, Mo-IV, W-IV.) X-ray diffraction (single-crystal and powder) measurements reveal that the [{(H2O)Fe(L-1)}{M(CN)(8)}{Fe(L-1)}](infinity) architectures consist of isomorphous 1-D polymeric structures based on the alternation of {Fe(L-1)}(2+) and {M(CN)(8)}(4-) units (L-1 stands for the pentadentate macrocycle). Analysis of the magnetic susceptibility behavior revealed cyano-bridged {Fe-Nb} exchange interaction to be antiferromagnetic with J = -20 cm(-1) deduced from fitting an Ising model taking into account the noncollinear spin arrangement. For this ferrimagnetic chain a slow relaxation of its magnetization is observed at low temperature revealing a SCM behavior with Delta/k(B) = 74 K and tau(0) = 4.6 x 10(-11) s. The M versus H behavior exhibits a hysteresis loop with a coercive field of 4 kOe at 1 K and reveals at 380 mK magnetic avalanche processes, i.e., abrupt reversals in magnetization as H is varied. The origin of these characteristics is attributed to the combination of efficient {Fe-Nb} exchange interaction and significant anisotropy of the {Fe(L-1)) unit. High field EPR and magnetization experiments have revealed for the parent compound [Fe(L-1)(H2O)(2)]Cl-2 a negative zero field splitting parameter of D approximate to -17 cm(-1). The crystal structure, magnetic behavior, and Mossbauer data for [Fe(L-1)(H2O)(2)]Cl-2 are also reported.

Enhanced Ion Anisotropy by Nonconventional Coordination Geometry: Single-Chain Magnet Behavior for a [{(FeL)-L-II}(2){Nb-IV(CN)(8)}] Helical Chain Compound Designed with Heptacoordinate Fe-II

Nonconventional heptacoordination in combination with efficient magnetic exchange coupling is shown to yield a 1-D heteronuclear {(FeNbIV)-Nb-II} compound with remarkable magnetic features when compared to other Fe(II)-based single chain magnets (SCM). Cyano-bridged heterometallic {3d-4d} and {3d-5d} chains are formed upon assembling Fe(II) bearing a pentadentate macrocycle as the blocking ligand with octacyano metallates, [M(CN)(8)](4-) (M = Nb-IV, Mo-IV, W-IV.) X-ray diffraction (single-crystal and powder) measurements reveal that the [{(H2O)Fe(L-1)}{M(CN)(8)}{Fe(L-1)}](infinity) architectures consist of isomorphous 1-D polymeric structures based on the alternation of {Fe(L-1)}(2+) and {M(CN)(8)}(4-) units (L-1 stands for the pentadentate macrocycle). Analysis of the magnetic susceptibility behavior revealed cyano-bridged {Fe-Nb} exchange interaction to be antiferromagnetic with J = -20 cm(-1) deduced from fitting an Ising model taking into account the noncollinear spin arrangement. For this ferrimagnetic chain a slow relaxation of its magnetization is observed at low temperature revealing a SCM behavior with Delta/k(B) = 74 K and tau(0) = 4.6 x 10(-11) s. The M versus H behavior exhibits a hysteresis loop with a coercive field of 4 kOe at 1 K and reveals at 380 mK magnetic avalanche processes, i.e., abrupt reversals in magnetization as H is varied. The origin of these characteristics is attributed to the combination of efficient {Fe-Nb} exchange interaction and significant anisotropy of the {Fe(L-1)) unit. High field EPR and magnetization experiments have revealed for the parent compound [Fe(L-1)(H2O)(2)]Cl-2 a negative zero field splitting parameter of D approximate to -17 cm(-1). The crystal structure, magnetic behavior, and Mossbauer data for [Fe(L-1)(H2O)(2)]Cl-2 are also reported.

New computational approaches for wrinkled and slack membranes

The static response of thin, wrinkled membranes is studied using both a tension field approximation based on plane stress conditions and a 3D nonlinear elasticityformulation, discretized through 8-noded Cosserat point elements. While the tension field approach only obtains the wrinkled/slack regions and at best a measure of the extent of wrinkliness, the 3D elasticity solution provides, in principle, the deformed shape of a wrinkled/slack membrane. However, since membranes barely resist compression, the discretized and linearized system equations via both the approaches are ill-conditioned and solutions could thus be sensitive to discretizations errors as well as other sources of noises/imperfections. We propose a regularized, pseudo-dynamical recursion scheme that provides a sequence of updates, which are almost insensitive to theregularizing term as well as the time step size used for integrating the pseudo-dynamical form. This is borne out through several numerical examples wherein the relative performance of the proposed recursion scheme vis-a-vis a regularized Newton strategy is compared. The pseudo-time marching strategy, when implemented using 3D Cosserat point elements, also provides a computationally cheaper, numerically accurate and simpler alternative to that using geometrically exact shell theories for computing large deformations of membranes in the presence of wrinkles. (C) 2010 Elsevier Ltd. All rights reserved.

Experimental and Computational Characterization of Alloy--Spinel--Corundum Equilibrium in the System Fe--Co--Al--O at 1873K

The three-phase equilibrium between alloy, spinel solid solution and alpha -Al sub 2 O sub 3 in the Fe--Co--Al--O system at 1873k was fully characterized as a function of alloy composition using both experimental and computational methods. The equilibrium oxygen content of the liquid alloy was measured by suction sampling and inert gas fusion analysis. The O potential corresponding to the three-phase equilibrium was determined by emf measurements on a solid state galvanic cell incorporating (Y sub 2 O sub 3 )ThO sub 2 as the solid electrolyte and Cr + Cr sub 2 O sub 3 as the reference electrode. The equilibrium composition of the spinel phase formed at the interface between the alloy and alumina crucible was measured by electron probe microanalysis (EPMA). The experimental results were compared with the values computed using a thermodynamic model. The model used values for standard Gibbs energies of formation of pure end-member spinels and Gibbs energies of solution of gaseous O in liquid Fe and cobalt available in the literature. The activity--composition relationship in the spinel solid solution was computed using a cation distribution model. The variation of the activity coefficient of O with alloy composition in the Fe--Co--O system was estimated using both the quasichemical model of Jacob and Alcock and Wagner's model along with the correlations of Chiang and Chang and Kuo and Chang. The computed results of spinel composition and O potential are in excellent agreement with the experimental data. Graphs. 29 ref.--AA

Geometry changes induced by negative hyperconjugative interactions involving carbonyl and thiocarbonyl groups

MNDO geometry optimizations have been carried out on a series of acyclic and cyclic unsymmetrically disubstituted carbonyl and thiocarbonyl compounds. The C=X unit shows a consistent and often sizeable tilt towards one of the substituents, following the order O > Snot, vert, similarN > C > B. Reference ab initio calculations and available experimental results support the MNDO results. The effect, which is particularly dramatic in small rings, is attributed primarily to favorable negative hyperconjugative interaction between the lone pair on X and a low lying adjacent σ* orbital. Such an interaction can lead to highly distorted structures, including perhaps to a planar molecule with an inverted sp2 carbon center.

sRNAscanner: A Computational Tool for Intergenic Small RNA Detection in Bacterial Genomes

Background:Bacterial non-coding small RNAs (sRNAs) have attracted considerable attention due to their ubiquitous nature and contribution to numerous cellular processes including survival, adaptation and pathogenesis. Existing computational approaches for identifying bacterial sRNAs demonstrate varying levels of success and there remains considerable room for improvement. Methodology/Principal Findings: Here we have proposed a transcriptional signal-based computational method to identify intergenic sRNA transcriptional units (TUs) in completely sequenced bacterial genomes. Our sRNAscanner tool uses position weight matrices derived from experimentally defined E. coli K-12 MG1655 sRNA promoter and rho-independent terminator signals to identify intergenic sRNA TUs through sliding window based genome scans. Analysis of genomes representative of twelve species suggested that sRNAscanner demonstrated equivalent sensitivity to sRNAPredict2, the best performing bioinformatics tool available presently. However, each algorithm yielded substantial numbers of known and uncharacterized hits that were unique to one or the other tool only. sRNAscanner identified 118 novel putative intergenic sRNA genes in Salmonella enterica Typhimurium LT2, none of which were flagged by sRNAPredict2. Candidate sRNA locations were compared with available deep sequencing libraries derived from Hfq-co-immunoprecipitated RNA purified from a second Typhimurium strain (Sittka et al. (2008) PLoS Genetics 4: e1000163). Sixteen potential novel sRNAs computationally predicted and detected in deep sequencing libraries were selected for experimental validation by Northern analysis using total RNA isolated from bacteria grown under eleven different growth conditions. RNA bands of expected sizes were detected in Northern blots for six of the examined candidates. Furthermore, the 5'-ends of these six Northern-supported sRNA candidates were successfully mapped using 5'-RACE analysis. Conclusions/Significance: We have developed, computationally examined and experimentally validated the sRNAscanner algorithm. Data derived from this study has successfully identified six novel S. Typhimurium sRNA genes. In addition, the computational specificity analysis we have undertaken suggests that similar to 40% of sRNAscanner hits with high cumulative sum of scores represent genuine, undiscovered sRNA genes. Collectively, these data strongly support the utility of sRNAscanner and offer a glimpse of its potential to reveal large numbers of sRNA genes that have to date defied identification. sRNAscanner is available from: http://bicmku.in:8081/sRNAscanner or http://cluster.physics.iisc.ernet.in/sRNAscanner/.

Investigation of the effect of nonlocal scale on ultrasonic wave dispersion characteristics of a monolayer graphene

In this article, an ultrasonic wave propagation in graphene sheet is studied using nonlocal elasticity theory incorporating small scale effects. The graphene sheet is modeled as an isotropic plate of one-atom thick. For this model, the nonlocal governing differential equations of motion are derived from the minimization of the total potential energy of the entire system. An ultrasonic type of wave propagation model is also derived for the graphene sheet. The nonlocal scale parameter introduces certain band gap region in in-plane and flexural wave modes where no wave propagation occurs. This is manifested in the wavenumber plots as the region where the wavenumber tends to infinite or wave speed tends to zero. The frequency at which this phenomenon occurs is called the escape frequency. The explicit expressions for cutoff frequencies and escape frequencies are derived. The escape frequencies are mainly introduced because of the nonlocal elasticity. Obviously these frequencies are function of nonlocal scaling parameter. It has also been obtained that these frequencies are independent of y-directional wavenumber. It means that for any type of nanostructure, the escape frequencies are purely a function of nonlocal scaling parameter only. It is also independent of the geometry of the structure. It has been found that the cutoff frequencies are function of nonlocal scaling parameter (e(0)a) and the y-directional wavenumber (k(y)). For a given nanostructure, nonlocal small scale coefficient can be obtained by matching the results from molecular dynamics (MD) simulations and the nonlocal elasticity calculations. At that value of the nonlocal scale coefficient, the waves will propagate in the nanostructure at that cut-off frequency. In the present paper, different values of e(o)a are used. One can get the exact e(0)a for a given graphene sheet by matching the MD simulation results of graphene with the results presented in this paper. (C) 2010 Elsevier B.V. All rights reserved.

New computational approaches for inelastic x-ray scattering

Inelastic x-ray scattering spectroscopy is a versatile experimental technique for probing the electronic structure of materials. It provides a wealth of information on the sample's atomic-scale structure, but extracting this information from the experimental data can be challenging because there is no direct relation between the structure and the measured spectrum. Theoretical calculations can bridge this gap by explaining the structural origins of the spectral features. Reliable methods for modeling inelastic x-ray scattering require accurate electronic structure calculations. This work presents the development and implementation of new schemes for modeling the inelastic scattering of x-rays from non-periodic systems. The methods are based on density functional theory and are applicable for a wide variety of molecular materials. Applications are presented in this work for amorphous silicon monoxide and several gas phase systems. Valuable new information on their structure and properties could be extracted with the combination of experimental and computational methods.

Computational Methods for Studies of Quantum Dots and Boson Condensates

This thesis presents ab initio studies of two kinds of physical systems, quantum dots and bosons, using two program packages of which the bosonic one has mainly been developed by the author. The implemented models, \emph{i.e.}, configuration interaction (CI) and coupled cluster (CC) take the correlated motion of the particles into account, and provide a hierarchy of computational schemes, on top of which the exact solution, within the limit of the single-particle basis set, is obtained. The theory underlying the models is presented in some detail, in order to provide insight into the approximations made and the circumstances under which they hold. Some of the computational methods are also highlighted. In the final sections the results are summarized. The CI and CC calculations on multiexciton complexes in self-assembled semiconductor quantum dots are presented and compared, along with radiative and non-radiative transition rates. Full CI calculations on quantum rings and double quantum rings are also presented. In the latter case, experimental and theoretical results from the literature are re-examined and an alternative explanation for the reported photoluminescence spectra is found. The boson program is first applied on a fictitious model system consisting of bosonic electrons in a central Coulomb field for which CI at the singles and doubles level is found to account for almost all of the correlation energy. Finally, the boson program is employed to study Bose-Einstein condensates confined in different anisotropic trap potentials. The effects of the anisotropy on the relative correlation energy is examined, as well as the effect of varying the interaction potential.}

Cross-field gas breakdown in a coaxial cylinder geometry

Sparking potentials in a coaxial cylinder geometry in oxygen and dry air were measured in crossed electric and magnetic fields. From the data effective collision frequencies were calculated using the equivalent pressure concept. It is shown that the equivalent pressure concept holds good for deriving the effective collision frequencies in non-uniform electric fields.

Computational Fluid Dynamics in Combustion

Computational fluid dynamics has reached a stage where flow field in practical situation can be predicted to aid the design and to probe into the fundamental flow physics to understand and resolve the issues in fundamental fluid mechanics The study examines the computation of reacting flows After exploring the conservation equations for species and energy, the methods of closing the reaction rate terms in turbulent flow have been examined briefly Two cases of computation where combustion-flow interaction plays important role, have been discussed to illustrate the computational aspects and the physical insight that can be gained by the reacting flow computation

HELANAL: A program to characterise helix geometry in proteins,

A detailed analysis of structural and position dependent characteristic features of helices will give a better understanding of the secondary structure formation in globular proteins. Here we describe an algorithm that quantifies the geometry of helices in proteins on the basis of their C-alpha atoms alone. The Fortran program HELANAL can extract the helices from the PDB files and then characterises the overall geometry of each helix as being linear, curved or kinked, in terms of its local structural features, viz. local helical twist and rise, virtual torsion angle, local helix origins and bending angles between successive local helix axes. Even helices with large radius of curvature are unambiguously identified as being linear or curved. The program can also be used to differentiate a kinked helix and other motifs, such as helix-loop-helix or a helix-turn-helix (with a single residue linker) with the help of local bending angles. In addition to these, the program can also be used to characterise the helix start and end as well as other types of secondary structures.