Modified conjugated gradient method for diagonalising large matrices

We present an iterative method to diagonalise large matrices. The basic idea is the same as the conjugated gradient (CG) method, i.e, minimizing the Rayleigh quotient via its gradient and avoiding reintroduce errors to the directions of previous gradients. Each iteration step is to find lowest eigenvector of the matrix in a subspace spanned by the current trial vector and the corresponding gradient of the Rayleigh quotient, as well as some previous trial vectors. The gradient, together with the previous trail vectors, play a similar role of the conjugated gradient of the original CG algorithm. Our numeric tests indicate that this method converges significantly faster than the original CG method. And the computational cost of one iteration step is about the same as the original CG method. It is suitably for first principle calculations.Comment: 6 Pages, 2EPS figures. (To appear in Phys. Rev. E

Growth of (110) Diamond using pure Dicarbon

We use a density-functional based tight-binding method to study diamond growth steps by depositing dicarbon species onto a hydrogen-free diamond (110) surface. Subsequent C_2 molecules are deposited on an initially clean surface, in the vicinity of a growing adsorbate cluster, and finally, near vacancies just before completion of a full new monolayer. The preferred growth stages arise from C_2n clusters in near ideal lattice positions forming zigzag chains running along the [-110] direction parallel to the surface. The adsorption energies are consistently exothermic by 8--10 eV per C_2, depending on the size of the cluster. The deposition barriers for these processes are in the range of 0.0--0.6 eV. For deposition sites above C_2n clusters the adsorption energies are smaller by 3 eV, but diffusion to more stable positions is feasible. We also perform simulations of the diffusion of C_2 molecules on the surface in the vicinity of existing adsorbate clusters using an augmented Lagrangian penalty method. We find migration barriers in excess of 3 eV on the clean surface, and 0.6--1.0 eV on top of graphene-like adsorbates. The barrier heights and pathways indicate that the growth from gaseous dicarbons proceeds either by direct adsorption onto clean sites or after migration on top of the existing C_2n chains.Comment: 8 Pages, 7 figure

Physical restraint in residential child care : the experiences of young people and residential workers

There have long been concerns about the use of physical restraint in residential care. This paper presents the findings of a qualitative study which explores the experiences of children, young people and residential workers about physical restraint. The research identifies the dilemmas and ambiguities for both staff and young people, and participants discuss the situations where they feel physical restraint is appropriate as well as their concerns about unjustified or painful restraints. They describe the negative emotions involved in restraint but also those situations where, through positive relationships and trust, restraint can help young people through unsafe situations

Efficient photoionization for barium ion trapping using a dipole-allowed resonant two-photon transition

Two efficient and isotope-selective resonant two-photon ionization techniques for loading barium ions into radio-frequency (RF)-traps are demonstrated. The scheme of using a strong dipole-allowed transition at \lambda=553 nm as a first step towards ionization is compared to the established technique of using a weak intercombination line (\lambda=413 nm). An increase of two orders of magnitude in the ionization efficiency is found favoring the transition at 553 nm. This technique can be implemented using commercial all-solid-state laser systems and is expected to be advantageous compared to other narrowband photoionization schemes of barium in cases where highest efficiency and isotope-selectivity are required.Comment: 8 pages, 5 figure

Gas-liquid transition in the model of particles interacting at high energy

An application of the ideas of the inertial confinement fusion process in the case of particles interacting at high energy is investigated. A possibility of the gas-liquid transition in the gas is considered using different approaches. In particular, a shock wave description of interactions between particles is studied and a self-similar solution of Euler's equation is discussed. Additionally, Boltzmann equation is solved for self-consistent field (Vlasov's equation) in linear approximation for the case of a gas under external pressure and the corresponding change of Knudsen number of the system is calculated.Comment: 24 pages, 2 figur

Vector meson production and nucleon resonance analysis in a coupled-channel approach for energies m_N < sqrt(s) < 2 GeV I: pion-induced results and hadronic parameters

We present a nucleon resonance analysis by simultaneously considering all pion- and photon-induced experimental data on the final states gamma N, pi N, 2 pi N, eta N, K Lambda, K Sigma, and omega N for energies from the nucleon mass up to sqrt(s) = 2 GeV. In this analysis we find strong evidence for the resonances P_{31}(1750), P_{13}(1900), P_{33}(1920), and D_{13}(1950). The omega N production mechanism is dominated by large P_{11}(1710) and P_{13}(1900) contributions. In this first part, we present the results of the pion-induced reactions and the extracted resonance and background properties with emphasis on the difference between global and purely hadronic fits.Comment: 54 pages, 26 figures, discussion extended, typos corrected, references updated, to appear in Phys. Rev.

Gauge and Scheme Dependence of Mixing Matrix Renormalization

We revisit the issue of mixing matrix renormalization in theories that include Dirac or Majorana fermions. We show how a gauge-variant on-shell renormalized mixing matrix can be related to a manifestly gauge-independent one within a generalized ${\bar {\rm MS}}$ scheme of renormalization. This scheme-dependent relation is a consequence of the fact that in any scheme of renormalization, the gauge-dependent part of the mixing-matrix counterterm is ultra-violet safe and has a pure dispersive form. Employing the unitarity properties of the theory, we can successfully utilize the afore-mentioned scheme-dependent relation to preserve basic global or local symmetries of the bare Lagrangian through the entire process of renormalization. As an immediate application of our study, we derive the gauge-independent renormalization-group equations of mixing matrices in a minimal extension of the Standard Model with isosinglet neutrinos.Comment: 31 pages, LaTeX, uses axodraw.st

SO(10) unified models and soft leptogenesis

Motivated by the fact that, in some realistic models combining SO(10) GUTs and flavour symmetries, it is not possible to achieve the required baryon asymmetry through the CP asymmetry generated in the decay of right-handed neutrinos, we take a fresh look on how deep this connection is in SO(10). The common characteristics of these models are that they use the see-saw with right-handed neutrinos, predict a normal hierarchy of masses for the neutrinos observed in oscillating experiments and in the basis where the right-handed Majorana mass is diagonal, the charged lepton mixings are tiny. In addition these models link the up-quark Yukawa matrix to the neutrino Yukawa matrix Y^\nu with the special feature of Y^\nu_{11}-> 0 Using this condition, we find that the required baryon asymmetry of the Universe can be explained by the soft leptogenesis using the soft B parameter of the second lightest right-handed neutrino whose mass turns out to be around 10^8 GeV. It is pointed out that a natural way to do so is to use no-scale supergravity where the value of B ~1 GeV is set through gauge-loop corrections.Comment: 26 pages, 2 figures. Added references, new appendix of a relevant fit and improved comment

Vector meson production and nucleon resonance analysis in a coupled-channel approach for energies m_N < sqrt(s) < 2 GeV II: photon-induced results

We present a nucleon resonance analysis by simultaneously considering all pion- and photon-induced experimental data on the final states gamma N, pi N, 2 pi N, eta N, K Lambda, K Sigma, and omega N for energies from the nucleon mass up to sqrt(s) = 2 GeV. In this analysis we find strong evidence for the resonances P_{31}(1750), P_{13}(1900), P_{33}(1920), and D_{13}(1950). The omega N production mechanism is dominated by large P_{11}(1710) and P_{13}(1900) contributions. In this second part we present the results on the photoproduction reactions and the electromagnetic properties of the resonances. The inclusion of all important final states up to sqrt(s) = 2 GeV allows for estimates on the importance of the individual states for the GDH sum rule.Comment: 41 pages, 26 figures, discussion extended, typos corrected, references updated, to appear in Phys. Rev.

Effective Lagrangian Approach to the Theory of Eta Photoproduction in the N∗(1535)N^{*}(1535) Region

We investigate eta photoproduction in the $N^{*}(1535)$ resonance region within the effective Lagrangian approach (ELA), wherein leading contributions to the amplitude at the tree level are taken into account. These include the nucleon Born terms and the leading $t$-channel vector meson exchanges as the non-resonant pieces. In addition, we consider five resonance contributions in the $s$- and $u$- channel; besides the dominant $N^{*}(1535)$, these are: $N^{*}(1440),N^{*}(1520),N^{*}(1650)$ and $N^{*}(1710)$. The amplitudes for the $\pi^\circ$ and the $\eta$ photoproduction near threshold have significant differences, even as they share common contributions, such as those of the nucleon Born terms. Among these differences, the contribution to the $\eta$ photoproduction of the $s$-channel excitation of the $N^{*}(1535)$ is the most significant. We find the off-shell properties of the spin-3/2 resonances to be important in determining the background contributions. Fitting our effective amplitude to the available data base allows us to extract the quantity $\sqrt{\chi \Gamma_\eta} A_{1/2}/\Gamma_T$, characteristic of the photoexcitation of the $N^{*}(1535)$ resonance and its decay into the $\eta$-nucleon channel, of interest to precise tests of hadron models. At the photon point, we determine it to be $(2.2\pm 0.2)\times 10^{-1} GeV^{-1}$ from the old data base, and $(2.2\pm 0.1) \times 10^{-1} GeV^{-1}$ from a combination of old data base and new Bates data. We obtain the helicity amplitude for $N^{*}(1535)\rightarrow \gamma p$ to be $A_{1/2}=(97\pm 7)\times 10^{-3} GeV^{-1/2}$ from the old data base, and $A_{1/2}=(97\pm 6)\times 10^{-3} GeV^{-1/2}$ from the combination of the old data base and new Bates data, compared with the results of the analysis of pion photoproduction yielding $74\pm 11$, in the same units.Comment: 43 pages, RevTeX, 9 figures available upon request, to appear in Phys. Rev.