Perfectly Matched Layers in a Divergence Preserving ADI Scheme for Electromagnetics

For numerical simulations of highly relativistic and transversely accelerated charged particles including radiation fast algorithms are needed. While the radiation in particle accelerators has wavelengths in the order of 100 um the computational domain has dimensions roughly 5 orders of magnitude larger resulting in very large mesh sizes. The particles are confined to a small area of this domain only. To resolve the smallest scales close to the particles subgrids are envisioned. For reasons of stability the alternating direction implicit (ADI) scheme by D. N. Smithe et al. (J. Comput. Phys. 228 (2009) pp.7289-7299) for Maxwell equations has been adopted. At the boundary of the domain absorbing boundary conditions have to be employed to prevent reflection of the radiation. In this paper we show how the divergence preserving ADI scheme has to be formulated in perfectly matched layers (PML) and compare the performance in several scenarios.Comment: 8 pages, 6 figure

Electrical Detection of Spin Accumulation at a Ferromagnet-Semiconductor Interface

We show that the accumulation of spin-polarized electrons at a forward-biased Schottky tunnel barrier between Fe and n-GaAs can be detected electrically. The spin accumulation leads to an additional voltage drop across the barrier that is suppressed by a small transverse magnetic field, which depolarizes the spins in the semiconductor. The dependence of the electrical accumulation signal on magnetic field, bias current, and temperature is in good agreement with the predictions of a drift-diffusion model for spin-polarized transport.Comment: Submitted to Phys. Rev. Let

Spin wave emission by spin-orbit torque antennas

We study the generation of propagating spin waves in Ta/CoFeB waveguides by spin-orbit torque antennas and compare them to conventional inductive antennas. The spin-orbit torque was generated by a transverse microwave current across the magnetic waveguide. The detected spin wave signals for an in-plane magnetization across the waveguide (Damon-Eshbach configuration) exhibited the expected phase rotation and amplitude decay upon propagation when the current spreading was taken into account. Wavevectors up to about 6 rad/$\mu$m could be excited by the spin-orbit torque antennas despite the current spreading, presumably due to the non-uniformity of the microwave current. The relative magnitude of generated anti-damping spin-Hall and Oersted fields was calculated within an analytic model and it was found that they contribute approximately equally to the total effective field generated by the spin-orbit torque antenna. Due to the ellipticity of the precession in the ultrathin waveguide and the different orientation of the anti-damping spin-Hall and Oersted fields, the torque was however still dominated by the Oersted field. The prospects for obtaining a pure spin-orbit torque response are discussed, as are the energy efficiency and the scaling properties of spin-orbit torque antennas.Comment: 20 pages, 5 figure

Beam Dynamics in High Intensity Cyclotrons Including Neighboring Bunch Effects: Model, Implementation and Application

Space charge effects, being one of the most significant collective effects, play an important role in high intensity cyclotrons. However, for cyclotrons with small turn separation, other existing effects are of equal importance. Interactions of radially neighboring bunches are also present, but their combined effects has not yet been investigated in any great detail. In this paper, a new particle in cell based self-consistent numerical simulation model is presented for the first time. The model covers neighboring bunch effects and is implemented in the three-dimensional object-oriented parallel code OPAL-cycl, a flavor of the OPAL framework. We discuss this model together with its implementation and validation. Simulation results are presented from the PSI 590 MeV Ring Cyclotron in the context of the ongoing high intensity upgrade program, which aims to provide a beam power of 1.8 MW (CW) at the target destination

Magnetisation distribution in the tetragonal phase of BaFe2As2

We have determined the spatial distribution of the magnetisation induced by a field of 9 T in the tetragonal phase of BaFe2As2 using polarised neutron diffraction. Magnetic structure factors derived from the polarisation dependence of the intensities of Bragg reflections were used to make a maximum entropy reconstruction of the distribution projected on the 110 plane. The reconstruction shows clearly that the magnetisation is confined to the region around the iron atoms and that there is no significant magnetisation associated with either the As or Ba atoms. The distribution of magnetisation around the Fe atom is significantly non-spherical with a shape which is extended in the directions in the projection. These results show that the electrons which give rise to the paramagnetic susceptibility are confined to the Fe atoms their distribution suggests that they occupy 3d t_2g type orbitals with about 60% in those of xy symmetry

Calorimetric Evidence of Multiband Superconductivity in Ba(Fe0.925Co0.075)2As2

We report on the determination of the electronic heat capacity of a slightly overdoped (x = 0.075) Ba(Fe1-xCox)2As2 single crystal with a Tc of 21.4 K. Our analysis of the temperature dependence of the superconducting-state specific heat provides strong evidence for a two-band s-wave order parameter with gap amplitudes 2D1(0)/kBTc=1.9 and 2D2(0)/kBTc=4.4. Our result is consistent with the recently predicted s+- order parameter [I. I. Mazin et al., Phys. Rev. Lett. 101, 057003 (2008)].Comment: 4 pages, 3 figure

A Parallel General Purpose Multi-Objective Optimization Framework, with Application to Beam Dynamics

Particle accelerators are invaluable tools for research in the basic and applied sciences, in fields such as materials science, chemistry, the biosciences, particle physics, nuclear physics and medicine. The design, commissioning, and operation of accelerator facilities is a non-trivial task, due to the large number of control parameters and the complex interplay of several conflicting design goals. We propose to tackle this problem by means of multi-objective optimization algorithms which also facilitate a parallel deployment. In order to compute solutions in a meaningful time frame a fast and scalable software framework is required. In this paper, we present the implementation of such a general-purpose framework for simulation-based multi-objective optimization methods that allows the automatic investigation of optimal sets of machine parameters. The implementation is based on a master/slave paradigm, employing several masters that govern a set of slaves executing simulations and performing optimization tasks. Using evolutionary algorithms as the optimizer and OPAL as the forward solver, validation experiments and results of multi-objective optimization problems in the domain of beam dynamics are presented. The high charge beam line at the Argonne Wakefield Accelerator Facility was used as the beam dynamics model. The 3D beam size, transverse momentum, and energy spread were optimized