Computationally-efficient stochastic cluster dynamics method for modeling damage accumulation in irradiated materials

An improved version of a recently developed stochastic cluster dynamics (SCD) method {[}Marian, J. and Bulatov, V. V., {\it J. Nucl. Mater.} \textbf{415} (2014) 84-95{]} is introduced as an alternative to rate theory (RT) methods for solving coupled ordinary differential equation (ODE) systems for irradiation damage simulations. SCD circumvents by design the curse of dimensionality of the variable space that renders traditional ODE-based RT approaches inefficient when handling complex defect population comprised of multiple (more than two) defect species. Several improvements introduced here enable efficient and accurate simulations of irradiated materials up to realistic (high) damage doses characteristic of next-generation nuclear systems. The first improvement is a procedure for efficiently updating the defect reaction-network and event selection in the context of a dynamically expanding reaction-network. Next is a novel implementation of the $\tau$-leaping method that speeds up SCD simulations by advancing the state of the reaction network in large time increments when appropriate. Lastly, a volume rescaling procedure is introduced to control the computational complexity of the expanding reaction-network through occasional reductions of the defect population while maintaining accurate statistics. The enhanced SCD method is then applied to model defect cluster accumulation in iron thin films subjected to triple ion-beam ($\text{Fe}^{3+}$, $\text{He}^{+}$ and $$\text{H\ensuremath{{}^{+}}}$$) irradiations, for which standard RT or spatially-resolved kinetic Monte Carlo simulations are prohibitively expensive

Temperature dependent photoluminescence of single CdS nanowires

Temperature dependent photoluminescence (PL) is used to study the electronic properties of single CdS nanowires. At low temperatures, both near-band edge (NBE) photoluminescence (PL) and spatially-localized defect-related PL are observed in many nanowires. The intensity of the defect states is a sensitive tool to judge the character and structural uniformity of nanowires. As the temperature is raised, the defect states rapidly quench at varying rates leaving the NBE PL which dominates up to room temperature. All PL lines from nanowires follow closely the temperature-dependent band edge, similar to that observed in bulk CdS.Comment: 11 pages, 4 figure

Type 2 solar radio events observed in the interplanetary medium. Part 1: General characteristics

Twelve type 2 solar radio events were observed in the 2 MHz to 30 kHz frequency range by the radio astronomy experiment on the ISEE-3 satellite over the period from September 1978 to December 1979. These data provide the most comprehensive sample of type 2 radio bursts observed at kilometer wavelengths. Dynamic spectra of a number of events are presented. Where possible, the 12 events were associated with an initiating flare, ground based radio data, the passage of a shock at the spacecraft, and the sudden commencement of a geomagnetic storm. The general characteristics of kilometric type 2 bursts are discussed

Low temperature photoluminescence imaging and time-resolved spectroscopy of single CdS nanowires

Time-resolved photoluminescence (PL) and micro-PL imaging were used to study single CdS nanowires at 10 K. The low-temperature PL of all CdS nanowires exhibit spectral features near energies associated with free and bound exciton transitions, with the transition energies and emission intensities varying along the length of the nanowire. In addition, several nanowires show spatially localized PL at lower energies which are associated with morphological irregularities in the nanowires. Time-resolved PL measurements indicate that exciton recombination in all CdS nanowires is dominated by non-radiative recombination at the surface of the nanowires.Comment: 9 pages, 3 figures, to be published in Applied Physics Letter

Selection and reconstruction of the top quarks in the all-hadronic decays at a Linear Collider

A method of reconstruction of the top quarks produced in the process E+E- -> t\bar{t} -> 6 jets at a Linear Collider (LC) is proposed. The approach does not involve a kinematic fit, as well as assumptions on the invariant masses of the dijets originating from the decays of W bosons and, therefore, the method is expected to be less sensitive to theoretical and experimental uncertainties on the top-mass measurement than traditional reconstruction methods. For the first time, the reconstruction of the top quarks was investigated using the full LC detector simulation after taking into account the background arising from QCD multi-jet production.Comment: 22 pages, including 13 figures and 3 table

SU(3)_LxU(1)_N Model for Right-Handed Neutrino Neutral Currents

A model based on the \mbox{SU(3)}_L\otimes \mbox{U(1)}_N gauge group, in which neutrinos have right-handed neutral currents is considered. We argue that in order to have a result consistent with low-energy one, the right-handed neutrino component must be treated as correction instead of an equivalent spin state.Comment: 6 pages, Latex, no figures, Accepted for publication in Phys. Rev.

Bayesian inverse problems for recovering coefficients of two scale elliptic equations

We consider the Bayesian inverse homogenization problem of recovering the locally periodic two scale coefficient of a two scale elliptic equation, given limited noisy information on the solution. We consider both the uniform and the Gaussian prior probability measures. We use the two scale homogenized equation whose solution contains the solution of the homogenized equation which describes the macroscopic behaviour, and the corrector which encodes the microscopic behaviour. We approximate the posterior probability by a probability measure determined by the solution of the two scale homogenized equation. We show that the Hellinger distance of these measures converges to zero when the microscale converges to zero, and establish an explicit convergence rate when the solution of the two scale homogenized equation is sufficiently regular. Sampling the posterior measure by Markov Chain Monte Carlo (MCMC) method, instead of solving the two scale equation using fine mesh for each proposal with extremely high cost, we can solve the macroscopic two scale homogenized equation. Although this equation is posed in a high dimensional tensorized domain, it can be solved with essentially optimal complexity by the sparse tensor product finite element method, which reduces the computational complexity of the MCMC sampling method substantially. We show numerically that observations on the macrosopic behaviour alone are not sufficient to infer the microstructure. We need also observations on the corrector. Solving the two scale homogenized equation, we get both the solution to the homogenized equation and the corrector. Thus our method is particularly suitable for sampling the posterior measure of two scale coefficients

Top Quark Pair Production close to Threshold: Top Mass, Width and Momentum Distribution

The complete NNLO QCD corrections to the total cross section $\sigma(e^+e^- \to Z*,\gamma*\to t\bar t)$ in the kinematic region close to the top-antitop threshold are calculated by solving the corresponding Schroedinger equations exactly in momentum space in a consistent momentum cutoff regularization scheme. The corrections coming from the same NNLO QCD effects to the top quark three-momentum distribution $d\sigma/d |\vec k_t|$ are determined. We discuss the origin of the large NNLO corrections to the peak position and the normalization of the total cross section observed in previous works and propose a new top mass definition, the 1S mass M_1S, which stabilizes the peak in the total cross section. If the influence of beamstrahlung and initial state radiation on the mass determination is small, a theoretical uncertainty on the 1S top mass measurement of 200 MeV from the total cross section at the linear collider seems possible. We discuss how well the 1S mass can be related to the $\bar{MS}$ mass. We propose a consistent way to implement the top quark width at NNLO by including electroweak effects into the NRQCD matching coefficients, which then can become complex.Comment: 53 pages, latex; minor changes, a number of typos correcte