Probing the flavor violating scalar top quark signal at the LHC

The Large Hadron Collider (LHC) has completed its run at 8 TeV with the experiments ATLAS and CMS having collected about 25 fb(-1) of data each. Discovery of a light Higgs boson coupled with lack of evidence for supersymmetry at the LHC so far, has motivated studies of supersymmetry in the context of naturalness with the principal focus being the third generation squarks. In this work, we analyze the prospects of the flavor violating decay mode (t) over tilde (1) -> c chi(0)(1) at 8 and 13 TeV center-of-mass energy at the LHC. This channel is also relevant in the dark matter context for the stop-coannihilation scenario, where the relic density depends on the mass difference between the lighter stop quark ((t) over tilde (1)) and the lightest neutralino (chi(0)(1)) states. This channel is extremely challenging to probe, especially for situations when the mass difference between the lighter stop quark and the lightest neutralino is small. Using certain kinematical properties of signal events we find that the level of backgrounds can be reduced substantially. We find that the prospect for this channel is limited due to the low production cross section for top squarks and limited luminosity at 8 TeV, but at the 13 TeV LHC with 100 fb(-1) luminosity, it is possible to probe top squarks with masses up to similar to 450 GeV. We also discuss how the sensitivity could be significantly improved by tagging charm jets.

Looking for BSM physics using top-quark polarization and decay-lepton kinematic asymmetries

We explore beyond-standard-model (BSM) physics signatures in the l + jets channel of the t (t) over bar pair production process at the Tevatron and the LHC. We study the effects of BSM physics scenarios on the top-quark polarization and on the kinematics of the decay leptons. To this end, we construct asymmetries using the lepton energy and angular distributions. Further, we find their correlations with the top polarization, net charge asymmetry and top forward-backward asymmetry. We show that when used together, these observables can help discriminate effectively between SM and different BSM scenarios, which can lead to varying degrees of top polarization at the Tevatron as well as the LHC. We use two types of colored mediator models to demonstrate the effectiveness of proposed observables, an s-channel axigluon and a u-channel diquark.

A Monte Carlo study of crystal structure transformations

The Metropolis algorithm has been generalized to allow for the variation of shape and size of the MC cell. A calculation using different potentials illustrates how the generalized method can be used for the study of crystal structure transformations. A restricted MC integration in the nine dimensional space of the cell components also leads to the stable structure for the Lennard-Jones potential.

Monte-Carlo estimation of time-dependent statistical characteristics of random dynamical systems

The problem of estimating the time-dependent statistical characteristics of a random dynamical system is studied under two different settings. In the first, the system dynamics is governed by a differential equation parameterized by a random parameter, while in the second, this is governed by a differential equation with an underlying parameter sequence characterized by a continuous time Markov chain. We propose, for the first time in the literature, stochastic approximation algorithms for estimating various time-dependent process characteristics of the system. In particular, we provide efficient estimators for quantities such as the mean, variance and distribution of the process at any given time as well as the joint distribution and the autocorrelation coefficient at different times. A novel aspect of our approach is that we assume that information on the parameter model (i.e., its distribution in the first case and transition probabilities of the Markov chain in the second) is not available in either case. This is unlike most other work in the literature that assumes availability of such information. Also, most of the prior work in the literature is geared towards analyzing the steady-state system behavior of the random dynamical system while our focus is on analyzing the time-dependent statistical characteristics which are in general difficult to obtain. We prove the almost sure convergence of our stochastic approximation scheme in each case to the true value of the quantity being estimated. We provide a general class of strongly consistent estimators for the aforementioned statistical quantities with regular sample average estimators being a specific instance of these. We also present an application of the proposed scheme on a widely used model in population biology. Numerical experiments in this framework show that the time-dependent process characteristics as obtained using our algorithm in each case exhibit excellent agreement with exact results. (C) 2010 Elsevier Inc. All rights reserved.

Monte Carlo simulation of network formation based on structural fragments in epoxy-anhydride systems

A method combining the Monte Carlo technique and the simple fragment approach has been developed for simulating network formation in amine-catalysed epoxy-anhydride systems. The method affords a detailed insight into the nature and composition of the network, showing the distribution of various fragments. It has been used to characterize the network formation in the reaction of the diglycidyl ester of isophthalic acid with hexahydrophthalic anhydride, catalysed by benzyldimethylamine. Pre-gel properties like number and weight distributions and average molecular weights have been calculated as a function of epoxy conversion, leading to a prediction of the gel-point conversion. Analysis of the simulated network further yields other characteristic properties such as concentration of crosslink points, distribution and concentration of elastically active chains, average molecular weight between crosslinks, sol content and mass fraction of pendent chains. A comparison has been made of the properties obtained through simulation with those predicted by the fragment approach alone, which, however, gives only average properties. The Monte Carlo simulation results clearly show that loops and other cyclic structures occur in the gel. This may account for the differences observed between the results of the simulation and the fragment model in the post-gel phase. Copyright (C) 1996 Elsevier Science Ltd.

A Kushner-Stratonovich Monte Carlo filter applied to nonlinear dynamical system identification

A Monte Carlo filter, based on the idea of averaging over characteristics and fashioned after a particle-based time-discretized approximation to the Kushner-Stratonovich (KS) nonlinear filtering equation, is proposed. A key aspect of the new filter is the gain-like additive update, designed to approximate the innovation integral in the KS equation and implemented through an annealing-type iterative procedure, which is aimed at rendering the innovation (observation prediction mismatch) for a given time-step to a zero-mean Brownian increment corresponding to the measurement noise. This may be contrasted with the weight-based multiplicative updates in most particle filters that are known to precipitate the numerical problem of weight collapse within a finite-ensemble setting. A study to estimate the a-priori error bounds in the proposed scheme is undertaken. The numerical evidence, presently gathered from the assessed performance of the proposed and a few other competing filters on a class of nonlinear dynamic system identification and target tracking problems, is suggestive of the remarkably improved convergence and accuracy of the new filter. (C) 2013 Elsevier B.V. All rights reserved.

Calculation of three-phase methane-ethane binary clathrate hydrate phase equilibrium from Monte Carlo molecular simulations

Methane and ethane are the simplest hydrocarbon molecules that can form clathrate hydrates. Previous studies have reported methods for calculating the three-phase equilibrium using Monte Carlo simulation methods in systems with a single component in the gas phase. Here we extend those methods to a binary gas mixture of methane and ethane. Methane-ethane system is an interesting one in that the pure components form sII clathrate hydrate whereas a binary mixture of the two can form the sII clathrate. The phase equilibria computed from Monte Carlo simulations show a good agreement with experimental data and are also able to predict the sI-sII structural transition in the clathrate hydrate. This is attributed to the quality of the TIP4P/Ice and TRaPPE models used in the simulations. (C) 2014 Elsevier B.V. All rights reserved.

Probabilistic stress variation studies on composite single lap joint using Monte Carlo simulation

The work presented in this paper involves the stochastic finite element analysis of composite-epoxy adhesive lap joints using Monte Carlo simulation. A set of composite adhesive lap joints were prepared and loaded till failure to obtain their strength. The peel and shear strain in the bond line region at different levels of load were obtained using digital image correlation (DIC). The corresponding stresses were computed assuming a plane strain condition. The finite element model was verified by comparing the numerical and experimental stresses. The stresses exhibited a similar behavior and a good correlation was obtained. Further, the finite element model was used to perform the stochastic analysis using Monte Carlo simulation. The parameters influencing stress distribution were provided as a random input variable and the resulting probabilistic variation of maximum peel and shear stresses were studied. It was found that the adhesive modulus and bond line thickness had significant influence on the maximum stress variation. While the adherend thickness had a major influence, the effect of variation in longitudinal and shear modulus on the stresses was found to be little. (C) 2014 Elsevier Ltd. All rights reserved.

Experimentally validated Raman Monte Carlo simulation for a cuboid object to obtain Raman spectroscopic signatures for hidden material

In conventional Raman spectroscopic measurements of liquids or surfaces the preferred geometry for detection of the Raman signal is the backscattering (or reflection) mode. For non-transparent layered materials, sub-surface Raman signals have been retrieved using spatially offset Raman spectroscopy (SORS), usually with light collection in the same plane as the point of excitation. However, as a result of multiple scattering in a turbid medium, Raman photons will be emitted in all directions. In this study, Monte Carlo simulations for a three-dimensional layered sample with finite geometry have been performed to confirm the detectability of Raman signals at all angles and at all sides of the object. We considered a non-transparent cuboid container (high density polyethylene) with explosive material (ammonium nitrate) inside. The simulation results were validated with experimental Raman intensities. Monte Carlo simulation results reveal that the ratio of sub-surface to surface signals improves at geometries other than backscattering. In addition, we demonstrate through simulations the effects of the absorption and scattering coefficients of the layers, and that of the diameter of the excitation beam. The advantage of collecting light from all possible 4 angles, over other collection modes, is that this technique is not geometry specific and molecular identification of layers underneath non-transparent surfaces can be obtained with minimal interference from the surface layer. To what extent all sides of the object will contribute to the total signal will depend on the absorption and scattering coefficients and the physical dimensions. Copyright (c) 2015 John Wiley & Sons, Ltd.

An investigation of solid adamantane by a modified isothermal-isobaric ensemble Monte Carlo simulation

Isothermal-isobaric ensemble Monte Carlo simulation studies of adamantane have been carried out at different temperatures. Thermodynamic properties and radial distribution functions calculated by employing a simple potential model based on sitesite interactions show good agreement with experiment and suggest that the solid is orientationally disordered at high temperatures.

Glass-Transition In Isopentane - A Monte-Carlo Study

Monte Carlo simulations with realistic interaction potentials have been carried out on isopentane to investigate the glass transition. Intermolecular pair-correlation functions of the glass show distinct differences from those of the liquid, the CH-CH pair-correlation function being uniquely different from the other pair-correlation functions. The coordination number of the glass is higher than that of the liquid, and the packing in the glass seems to be mainly governed by the geometrical constraints of the molecule. Annealing affects the properties of the glass significantly.

A conditionally linearized Monte Carlo filter in non-linear structural dynamics

State and parameter estimations of non-linear dynamical systems, based on incomplete and noisy measurements, are considered using Monte Carlo simulations. Given the measurements. the proposed method obtains the marginalized posterior distribution of an appropriately chosen (ideally small) subset of the state vector using a particle filter. Samples (particles) of the marginalized states are then used to construct a family of conditionally linearized system of equations and thus obtain the posterior distribution of the states using a bank of Kalman filters. Discrete process equations for the marginalized states are derived through truncated Ito-Taylor expansions. Increased analyticity and reduced dispersion of weights computed over a smaller sample space of marginalized states are the key features of the filter that help achieve smaller sample variance of the estimates. Numerical illustrations are provided for state/parameter estimations of a Duffing oscillator and a 3-DOF non-linear oscillator. Performance of the filter in parameter estimation is also assessed using measurements obtained through experiments on simple models in the laboratory. Despite an added computational cost, the results verify that the proposed filter generally produces estimates with lower sample variance over the standard sequential importance sampling (SIS) filter.

Monte Carlo simulation of mictomagnets

Monte Carlo simulations of a binary alloy with impurity concentrations between 20 and 45 at.% have been carried out. The proportion of large clusters relative to that of small clusters increases with the number of MC diffusion steps as well as impurity concentration. Magnetic susceptibility peaks become more prominent and occur at higher temperatures with increasing impurity concentration. The different peaks in the susceptibility and specific heat curves seem to correspond to different sized clusters. A freezing model would explain the observed behaviour with the large clusters freezing first and the small clusters contributing to susceptibility (specific heat) peaks at lower temperatures.

Structure, energetics and diffusion properties of isomers of trimethyl benzene in beta zeolite: Uptake and Monte Carlo simulation study

A Monte Carlo study along with experimental uptake measurements of 1,2,3-trimethyl benzene, 1,2,4-trimethyl benzene and 1,3,5-trimethyl benzene (TMB) in beta zeolite is reported. The TraPPE potential has been employed for hydrocarbon interaction and harmonic potential of Demontis for modeling framework of the zeolite. Structure, energetics and dynamics of TMB in zeolite beta from Monte Carlo runs reveal interesting information about the diameter, properties of these isomers on confinement. Of the three isomers, 135TMB is supposed to have the largest diameter. It is seen TraPPE with Demontis potential predicts a restricted motion of 135TMB in the channels of zeolite beta.Experimentally, 135TMB has the highest transport diffusivity whereas MID results suggest this has the lowest self diffusivity. (C) 2009 Elsevier Inc. Ail rights reserved.

Monte Carlo simulation of Ising chains showing several phase transitions

A Monte Carlo simulation of Ising chains with competing short-range and infiniterange interactions has been carried out. Results show that whenever the system does not enter a metastable state, variation of temperature brings about phase transitions in the Ising chain. These phase transitions, except for two sets of interaction strengths, are generally of higher order and involve changes in the long-range order while the short-range order remains unaffected.