A method for atomistic spin dynamics simulations: implementation and examples

We present a method for performing atomistic spin dynamic simulations. A comprehensive summary of all pertinent details for performing the simulations such as equations of motions, models for including temperature, methods of extracting data and numerical schemes for performing the simulations is given. The method can be applied in a first principles mode, where all interatomic exchange is calculated self-consistently, or it can be applied with frozen parameters estimated from experiments or calculated for a fixed spin-configuration. Areas of potential applications to different magnetic questions are also discussed. The method is finally applied to one situation where the macrospin model breaks down; magnetic switching in ultra strong fields.Comment: 14 pages, 19 figure

Atomistic spin dynamics of the CuMn spin glass alloy

We demonstrate the use of Langevin spin dynamics for studying dynamical properties of an archetypical spin glass system. Simulations are performed on CuMn (20% Mn) where we study the relaxation that follows a sudden quench of the system to the low temperature phase. The system is modeled by a Heisenberg Hamiltonian where the Heisenberg interaction parameters are calculated by means of first-principles density functional theory. Simulations are performed by numerically solving the Langevin equations of motion for the atomic spins. It is shown that dynamics is governed, to a large degree, by the damping parameter in the equations of motion and the system size. For large damping and large system sizes we observe the typical aging regime.Comment: 18 pages, 9 figure

Dynamics of diluted magnetic semiconductors from atomistic spin dynamics simulations: Mn doped GaAs as a case study

The dynamical behavior of the magnetism of diluted magnetic semiconductors (DMS) has been investigated by means of atomistic spin dynamics simulations. The conclusions drawn from the study are argued to be general for DMS systems in the low concentration limit, although all simulations are done for 5% Mn-doped GaAs with various concentrations of As antisite defects. The magnetization curve, $M(T)$, and the Curie temperature $T_C$ have been calculated, and are found to be in good correspondence to results from Monte Carlo simulations and experiments. Furthermore, equilibrium and non-equilibrium behavior of the magnetic pair correlation function have been extracted. The dynamics of DMS systems reveals a substantial short ranged magnetic order even at temperatures at or above the ordering temperature, with a non-vanishing pair correlation function extending up to several atomic shells. For the high As antisite concentrations the simulations show a short ranged anti-ferromagnetic coupling, and a weakened long ranged ferromagnetic coupling. For sufficiently large concentrations we do not observe any long ranged ferromagnetic correlation. A typical dynamical response shows that starting from a random orientation of moments, the spin-correlation develops very fast ($\sim$ 1ps) extending up to 15 atomic shells. Above $\sim$ 10 ps in the simulations, the pair correlation is observed to extend over some 40 atomic shells. The autocorrelation function has been calculated and compared with ferromagnets like bcc Fe and spin-glass materials. We find no evidence in our simulations for a spin-glass behaviour, for any concentration of As antisites. Instead the magnetic response is better described as slow dynamics, at least when compared to that of a regular ferromagnet like bcc Fe.Comment: 24 pages, 15 figure

Performance of an Operating High Energy Physics Data Grid: D0SAR-Grid

The D0 experiment at Fermilab's Tevatron will record several petabytes of data over the next five years in pursuing the goals of understanding nature and searching for the origin of mass. Computing resources required to analyze these data far exceed capabilities of any one institution. Moreover, the widely scattered geographical distribution of D0 collaborators poses further serious difficulties for optimal use of human and computing resources. These difficulties will exacerbate in future high energy physics experiments, like the LHC. The computing grid has long been recognized as a solution to these problems. This technology is being made a more immediate reality to end users in D0 by developing a grid in the D0 Southern Analysis Region (D0SAR), D0SAR-Grid, using all available resources within it and a home-grown local task manager, McFarm. We will present the architecture in which the D0SAR-Grid is implemented, the use of technology and the functionality of the grid, and the experience from operating the grid in simulation, reprocessing and data analyses for a currently running HEP experiment.Comment: 3 pages, no figures, conference proceedings of DPF04 tal

Thermally activated magnetization reversal in monoatomic magnetic chains on surfaces studied by classical atomistic spin-dynamics simulations

We analyze the spontaneous magnetization reversal of supported monoatomic chains of finite length due to thermal fluctuations via atomistic spin-dynamics simulations. Our approach is based on the integration of the Landau-Lifshitz equation of motion of a classical spin Hamiltonian at the presence of stochastic forces. The associated magnetization lifetime is found to obey an Arrhenius law with an activation barrier equal to the domain wall energy in the chain. For chains longer than one domain-wall width, the reversal is initiated by nucleation of a reversed magnetization domain primarily at the chain edge followed by a subsequent propagation of the domain wall to the other edge in a random-walk fashion. This results in a linear dependence of the lifetime on the chain length, if the magnetization correlation length is not exceeded. We studied chains of uniaxial and tri-axial anisotropy and found that a tri-axial anisotropy leads to a reduction of the magnetization lifetime due to a higher reversal attempt rate, even though the activation barrier is not changed.Comment: 2nd version contains some improvements and new Appendi

Study of the Radiation-Hardness of VCSEL and PIN

The silicon trackers of the ATLAS experiment at the Large Hadron Collider (LHC) at CERN (Geneva) use optical links for data transmission. An upgrade of the trackers is planned for the Super LHC (SLHC), an upgraded LHC with ten times higher luminosity. We study the radiation-hardness of VCSELs (Vertical-Cavity Surface-Emitting Laser) and GaAs and silicon PINs using 24 GeV/c protons at CERN for possible application in the data transmission upgrade. The optical power of VCSEL arrays decreases significantly after the irradiation but can be partially annealed with high drive currents. The responsivities of the PIN diodes also decrease significantly after irradiation, but can be recovered by operating at higher bias voltage. This provides a simple mechanism to recover from the radiation damage

Observation of the Hadronic Transitions Chi_{b 1,2}(2P) -> omega Upsilon(1S)

The CLEO Collaboration has observed the first hadronic transition among bottomonium (b bbar) states other than the dipion transitions among vector states, Upsilon(nS) -> pi pi Upsilon(mS). In our study of Upsilon(3S) decays, we find a significant signal for Upsilon(3S) -> gamma omega Upsilon(1S) that is consistent with radiative decays Upsilon(3S) -> gamma chi_{b 1,2}(2P), followed by chi_{b 1,2} -> omega Upsilon(1S). The branching ratios we obtain are Br(chi_{b1} -> omega Upsilon(1S) = 1.63 (+0.35 -0.31) (+0.16 -0.15) % and Br(chi_{b2} -> omega Upsilon(1S) = 1.10 (+0.32 -0.28) (+0.11 - 0.10)%, in which the first error is statistical and the second is systematic.Comment: submitted to XXI Intern'l Symp on Lepton and Photon Interact'ns at High Energies, August 2003, Fermila

Improved Measurement of the Form Factors in the Decay Lambda_c^+ --> Lambda e^+ nu_e

Using the CLEO detector at the Cornell Electron Storage Ring, we have studied the distribution of kinematic variables in the decay Lambda_c^+ -> Lambda e^+ nu_e. By performing a four-dimensional maximum likelihood fit, we determine the form factor ratio, R = f_2/f_1 = -0.31 +/- 0.05(stat) +/- 0.04(syst), the pole mass, M_{pole} = (2.21 +/- 0.08(stat) +/- 0.14(syst)) GeV/c^2, and the decay asymmetry parameter of the Lambda_c, alpha_{Lambda_c} = -0.86 +/- 0.03(stat) +/- 0.02(syst), for = 0.67 (GeV/c^2)^2. We compare the angular distributions of the Lambda_c^+ and Lambda_c^- and find no evidence for CP-violation: A_{Lambda_c} = (alpha_{Lambda_c^+} + alpha_{Lambda_c^-})/ (alpha_{Lambda_c^+} - alpha_{Lambda_c^-}) = 0.00 +/- 0.03(stat) +/- 0.01(syst) +/- 0.02, where the third error is from the uncertainty in the world average of the CP-violating parameter, A_{Lambda}, for Lambda -> p pi^-.Comment: 8 pages postscript,also available through http://www.lns.cornell.edu/public/CLNS/2004/, submitted to PR

Measurement of the B-Meson Inclusive Semileptonic Branching Fraction and Electron-Energy Moments

We report a new measurement of the B-meson semileptonic decay momentum spectrum that has been made with a sample of 9.4/fb of electron-positron annihilation data collected with the CLEO II detector at the Y(4S) resonance. Electrons from primary semileptonic decays and secondary charm decays were separated by using charge and angular correlations in Y(4S) events with a high-momentum lepton and an additional electron. We determined the semileptonic branching fraction to be (10.91 +- 0.09 +- 0.24)% from the normalization of the electron-energy spectrum. We also measured the moments of the electron energy spectrum with minimum energies from 0.6 GeV to 1.5 GeV.Comment: 36 pages postscript, als available through http://w4.lns.cornell.edu/public/CLNS/, Submitted to PRD (back-to-back with preceding preprint hep-ex/0403052