Revision of FMM-Yukawa: An adaptive fast multipole method for screened Coulomb interactions
Zhang, Bo; Huang, Jingfang; Pitsianis, Nikos P.; Sun, Xiaobai
2010-12-01
FMM-YUKAWA is a mathematical software package primarily for rapid evaluation of the screened Coulomb interactions of N particles in three dimensional space. Since its release, we have revised and re-organized the data structure, software architecture, and user interface, for the purpose of enabling more flexible, broader and easier use of the package. The package and its documentation are available at http://www.fastmultipole.org/, along with a few other closely related mathematical software packages. New version program summaryProgram title: FMM-Yukawa Catalogue identifier: AEEQ_v2_0 Program summary URL:http://cpc.cs.qub.ac.uk/summaries/AEEQ_v2_0.html Program obtainable from: CPC Program Library, Queen's University, Belfast, N. Ireland Licensing provisions: GNU GPL 2.0 No. of lines in distributed program, including test data, etc.: 78 704 No. of bytes in distributed program, including test data, etc.: 854 265 Distribution format: tar.gz Programming language: FORTRAN 77, FORTRAN 90, and C. Requires gcc and gfortran version 4.4.3 or later Computer: All Operating system: Any Classification: 4.8, 4.12 Catalogue identifier of previous version: AEEQ_v1_0 Journal reference of previous version: Comput. Phys. Comm. 180 (2009) 2331 Does the new version supersede the previous version?: Yes Nature of problem: To evaluate the screened Coulomb potential and force field of N charged particles, and to evaluate a convolution type integral where the Green's function is the fundamental solution of the modified Helmholtz equation. Solution method: The new version of fast multipole method (FMM) that diagonalizes the multipole-to-local translation operator is applied with the tree structure adaptive to sample particle locations. Reasons for new version: To handle much larger particle ensembles, to enable the iterative use of the subroutines in a solver, and to remove potential contention in assignments for parallelization. Summary of revisions: The software package FMM-Yukawa has been
An implementation of the fast multiple method (FMM) is performed for magnetic systems with long-ranged dipolar interactions. Expansion in spherical harmonics of the original FMM is replaced by expansion of polynomials in Cartesian coordinates, which is considerably simpler. Under open boundary conditions, an expression for multipole moments of point dipoles in a cell is derived. These make the program appropriate for nanomagnetic simulations, including magnetic nanoparticles and ferrofluids. The performance is optimized in terms of cell size and parameter set (expansion order and opening angle) and the trade off between computing time and accuracy is quantitatively studied. A rule of thumb is proposed to decide the appropriate average number of dipoles in the smallest cells, and an optimal choice of parameter set is suggested. Finally, the superiority of Cartesian coordinate FMM is demonstrated by comparison to spherical harmonics FMM and FFT.
Zhang, Wen
2008-01-01
An implementation of the fast multiple method (FMM) is performed for magnetic systems with long-ranged dipolar interactions. Expansion in spherical harmonics of the original FMM is replaced by expansion of polynomials in cartesian coordinates, which is considerably simpler. Under open boundary conditions, an expression for multipole moments of point dipoles in a cell is derived. These make the program appropriate for nanomagnetic simulations, including magnetic nanoparticles and ferrofluids. The performance is optimized in terms of cell size and parameter set (expansion order and opening angle) and trade off between computing time and accuracy is quantitatively studied. A rule of thumb is proposed to decide the average number of dipoles in the smallest cells, and an optimal choice of parameter set is suggested.
Simple recursive implementation of fast multipole method
In this paper we present an implementation of the well known 'fast multipole' method (FMM) for the efficient calculation of dipole fields. The main advantage of the present implementation is simplicity-we believe that a major reason for the lack of use of FMMs is their complexity. One of the simplifications is the use of polynomials in the Cartesian coordinates rather than spherical harmonics. We have implemented it in the context of an arbitrary hierarchical system of cells-no periodic mesh is required, as it is for FFT (fast Fourier transform) methods. The implementation is in terms of recursive functions. Results are given for application to micromagnetic simulation. Complete source code is provided for an open-source implementation of this method, as well as an installer for the resulting program.
Scalable force directed graph layout algorithms using fast multipole methods
Yunis, Enas
2012-06-01
We present an extension to ExaFMM, a Fast Multipole Method library, as a generalized approach for fast and scalable execution of the Force-Directed Graph Layout algorithm. The Force-Directed Graph Layout algorithm is a physics-based approach to graph layout that treats the vertices V as repelling charged particles with the edges E connecting them acting as springs. Traditionally, the amount of work required in applying the Force-Directed Graph Layout algorithm is O(|V|2 + |E|) using direct calculations and O(|V| log |V| + |E|) using truncation, filtering, and/or multi-level techniques. Correct application of the Fast Multipole Method allows us to maintain a lower complexity of O(|V| + |E|) while regaining most of the precision lost in other techniques. Solving layout problems for truly large graphs with millions of vertices still requires a scalable algorithm and implementation. We have been able to leverage the scalability and architectural adaptability of the ExaFMM library to create a Force-Directed Graph Layout implementation that runs efficiently on distributed multicore and multi-GPU architectures. © 2012 IEEE.
Scalable fast multipole accelerated vortex methods
Hu, Qi
2014-05-01
The fast multipole method (FMM) is often used to accelerate the calculation of particle interactions in particle-based methods to simulate incompressible flows. To evaluate the most time-consuming kernels - the Biot-Savart equation and stretching term of the vorticity equation, we mathematically reformulated it so that only two Laplace scalar potentials are used instead of six. This automatically ensuring divergence-free far-field computation. Based on this formulation, we developed a new FMM-based vortex method on heterogeneous architectures, which distributed the work between multicore CPUs and GPUs to best utilize the hardware resources and achieve excellent scalability. The algorithm uses new data structures which can dynamically manage inter-node communication and load balance efficiently, with only a small parallel construction overhead. This algorithm can scale to large-sized clusters showing both strong and weak scalability. Careful error and timing trade-off analysis are also performed for the cutoff functions induced by the vortex particle method. Our implementation can perform one time step of the velocity+stretching calculation for one billion particles on 32 nodes in 55.9 seconds, which yields 49.12 Tflop/s.
Performance Benchmarking of Fast Multipole Methods
Al-Harthi, Noha A.
2013-06-01
The current trends in computer architecture are shifting towards smaller byte/flop ratios, while available parallelism is increasing at all levels of granularity – vector length, core count, and MPI process. Intel’s Xeon Phi coprocessor, NVIDIA’s Kepler GPU, and IBM’s BlueGene/Q all have a Byte/flop ratio close to 0.2, which makes it very difficult for most algorithms to extract a high percentage of the theoretical peak flop/s from these architectures. Popular algorithms in scientific computing such as FFT are continuously evolving to keep up with this trend in hardware. In the meantime it is also necessary to invest in novel algorithms that are more suitable for computer architectures of the future. The fast multipole method (FMM) was originally developed as a fast algorithm for ap- proximating the N-body interactions that appear in astrophysics, molecular dynamics, and vortex based fluid dynamics simulations. The FMM possesses have a unique combination of being an efficient O(N) algorithm, while having an operational intensity that is higher than a matrix-matrix multiplication. In fact, the FMM can reduce the requirement of Byte/flop to around 0.01, which means that it will remain compute bound until 2020 even if the cur- rent trend in microprocessors continues. Despite these advantages, there have not been any benchmarks of FMM codes on modern architectures such as Xeon Phi, Kepler, and Blue- Gene/Q. This study aims to provide a comprehensive benchmark of a state of the art FMM code “exaFMM” on the latest architectures, in hopes of providing a useful reference for deciding when the FMM will become useful as the computational engine in a given application code. It may also serve as a warning to certain problem size domains areas where the FMM will exhibit insignificant performance improvements. Such issues depend strongly on the asymptotic constants rather than the asymptotics themselves, and therefore are strongly implementation and hardware
A new simple multidomain fast multipole boundary element method
Huang, S.; Liu, Y. J.
2016-06-01
A simple multidomain fast multipole boundary element method (BEM) for solving potential problems is presented in this paper, which can be applied to solve a true multidomain problem or a large-scale single domain problem using the domain decomposition technique. In this multidomain BEM, the coefficient matrix is formed simply by assembling the coefficient matrices of each subdomain and the interface conditions between subdomains without eliminating any unknown variables on the interfaces. Compared with other conventional multidomain BEM approaches, this new approach is more efficient with the fast multipole method, regardless how the subdomains are connected. Instead of solving the linear system of equations directly, the entire coefficient matrix is partitioned and decomposed using Schur complement in this new approach. Numerical results show that the new multidomain fast multipole BEM uses fewer iterations in most cases with the iterative equation solver and less CPU time than the traditional fast multipole BEM in solving large-scale BEM models. A large-scale fuel cell model with more than 6 million elements was solved successfully on a cluster within 3 h using the new multidomain fast multipole BEM.
A new simple multidomain fast multipole boundary element method
Huang, S.; Liu, Y. J.
2016-09-01
A simple multidomain fast multipole boundary element method (BEM) for solving potential problems is presented in this paper, which can be applied to solve a true multidomain problem or a large-scale single domain problem using the domain decomposition technique. In this multidomain BEM, the coefficient matrix is formed simply by assembling the coefficient matrices of each subdomain and the interface conditions between subdomains without eliminating any unknown variables on the interfaces. Compared with other conventional multidomain BEM approaches, this new approach is more efficient with the fast multipole method, regardless how the subdomains are connected. Instead of solving the linear system of equations directly, the entire coefficient matrix is partitioned and decomposed using Schur complement in this new approach. Numerical results show that the new multidomain fast multipole BEM uses fewer iterations in most cases with the iterative equation solver and less CPU time than the traditional fast multipole BEM in solving large-scale BEM models. A large-scale fuel cell model with more than 6 million elements was solved successfully on a cluster within 3 h using the new multidomain fast multipole BEM.
Gaussian translation operator for Multi-Level Fast Multipole Method
Borries, Oscar Peter; Hansen, Per Christian; Sorensen, Stig B.; Meincke, Peter; Jorgensen, Erik
Results using a new translation operator for the Multi-Level Fast Multipole Method are presented. Based on Gaussian beams, the translation operator allows a significant portion of the plane-wave directions to be neglected, resulting in a much faster translation step.......Results using a new translation operator for the Multi-Level Fast Multipole Method are presented. Based on Gaussian beams, the translation operator allows a significant portion of the plane-wave directions to be neglected, resulting in a much faster translation step....
Chang-Jun Zheng; Hai-Bo Chen; Lei-Lei Chen
2013-01-01
This paper presents a novel wideband fast multipole boundary element approach to 3D half-space/planesymmetric acoustic wave problems.The half-space fundamental solution is employed in the boundary integral equations so that the tree structure required in the fast multipole algorithm is constructed for the boundary elements in the real domain only.Moreover,a set of symmetric relations between the multipole expansion coefficients of the real and image domains are derived,and the half-space fundamental solution is modified for the purpose of applying such relations to avoid calculating,translating and saving the multipole/local expansion coefficients of the image domain.The wideband adaptive multilevel fast multipole algorithm associated with the iterative solver GMRES is employed so that the present method is accurate and efficient for both lowand high-frequency acoustic wave problems.As for exterior acoustic problems,the Burton-Miller method is adopted to tackle the fictitious eigenfrequency problem involved in the conventional boundary integral equation method.Details on the implementation of the present method are described,and numerical examples are given to demonstrate its accuracy and efficiency.
Multilevel Fast Multipole Method for Higher Order Discretizations
Borries, Oscar Peter; Meincke, Peter; Jorgensen, Erik; Hansen, Per Christian
2014-01-01
The multi-level fast multipole method (MLFMM) for a higher order (HO) discretization is demonstrated on high-frequency (HF) problems, illustrating for the first time how an efficient MLFMM for HO can be achieved even for very large groups. Applying several novel ideas, beneficial to both lower...
Data-driven execution of fast multipole methods
Ltaief, Hatem
2013-09-17
Fast multipole methods (FMMs) have O (N) complexity, are compute bound, and require very little synchronization, which makes them a favorable algorithm on next-generation supercomputers. Their most common application is to accelerate N-body problems, but they can also be used to solve boundary integral equations. When the particle distribution is irregular and the tree structure is adaptive, load balancing becomes a non-trivial question. A common strategy for load balancing FMMs is to use the work load from the previous step as weights to statically repartition the next step. The authors discuss in the paper another approach based on data-driven execution to efficiently tackle this challenging load balancing problem. The core idea consists of breaking the most time-consuming stages of the FMMs into smaller tasks. The algorithm can then be represented as a directed acyclic graph where nodes represent tasks and edges represent dependencies among them. The execution of the algorithm is performed by asynchronously scheduling the tasks using the queueing and runtime for kernels runtime environment, in a way such that data dependencies are not violated for numerical correctness purposes. This asynchronous scheduling results in an out-of-order execution. The performance results of the data-driven FMM execution outperform the previous strategy and show linear speedup on a quad-socket quad-core Intel Xeon system.Copyright © 2013 John Wiley & Sons, Ltd. Copyright © 2013 John Wiley & Sons, Ltd.
Electrical impedance tomography and the fast multipole method
Bikowski, Jutta; Mueller, Jennifer L.
2004-10-01
A 3-D linearization-based reconstruction algorithm for Electrical Impedance Tomography suitable for breast cancer detection using data collected on a rectangular array was introduced by Mueller et al. [IEEE Biomed. Eng., 46(11), 1999]. By considering the scenario as an electrostatic problem, it is possible to model the electrodes with various charges, facilitating the use of the Fast Multipole Method (FMM) for calculating particle interactions and also supporting the use of different electrode models. In this paper the use of FMM is explained and results in form of reconstructed images from experimental data show that this method is an improvement.
Fourier-Based Fast Multipole Method for the Helmholtz Equation
Cecka, Cris
2013-01-01
The fast multipole method (FMM) has had great success in reducing the computational complexity of solving the boundary integral form of the Helmholtz equation. We present a formulation of the Helmholtz FMM that uses Fourier basis functions rather than spherical harmonics. By modifying the transfer function in the precomputation stage of the FMM, time-critical stages of the algorithm are accelerated by causing the interpolation operators to become straightforward applications of fast Fourier transforms, retaining the diagonality of the transfer function, and providing a simplified error analysis. Using Fourier analysis, constructive algorithms are derived to a priori determine an integration quadrature for a given error tolerance. Sharp error bounds are derived and verified numerically. Various optimizations are considered to reduce the number of quadrature points and reduce the cost of computing the transfer function. © 2013 Society for Industrial and Applied Mathematics.
Pipelining the Fast Multipole Method over a Runtime System
Agullo, Emmanuel; Coulaud, Olivier; Darve, Eric; Messner, Matthias; Toru, Takahashi
2012-01-01
Fast Multipole Methods (FMM) are a fundamental operation for the simulation of many physical problems. The high performance design of such methods usually requires to carefully tune the algorithm for both the targeted physics and the hardware. In this paper, we propose a new approach that achieves high performance across architectures. Our method consists of expressing the FMM algorithm as a task flow and employing a state-of-the-art runtime system, StarPU, in order to process the tasks on the different processing units. We carefully design the task flow, the mathematical operators, their Central Processing Unit (CPU) and Graphics Processing Unit (GPU) implementations, as well as scheduling schemes. We compute potentials and forces of 200 million particles in 48.7 seconds on a homogeneous 160 cores SGI Altix UV 100 and of 38 million particles in 13.34 seconds on a heterogeneous 12 cores Intel Nehalem processor enhanced with 3 Nvidia M2090 Fermi GPUs.
Scalable fast multipole methods for vortex element methods
Hu, Qi
2012-11-01
We use a particle-based method to simulate incompressible flows, where the Fast Multipole Method (FMM) is used to accelerate the calculation of particle interactions. The most time-consuming kernelsâ\\'the Biot-Savart equation and stretching term of the vorticity equationâ\\'are mathematically reformulated so that only two Laplace scalar potentials are used instead of six, while automatically ensuring divergence-free far-field computation. Based on this formulation, and on our previous work for a scalar heterogeneous FMM algorithm, we develop a new FMM-based vortex method capable of simulating general flows including turbulence on heterogeneous architectures, which distributes the work between multi-core CPUs and GPUs to best utilize the hardware resources and achieve excellent scalability. The algorithm also uses new data structures which can dynamically manage inter-node communication and load balance efficiently but with only a small parallel construction overhead. This algorithm can scale to large-sized clusters showing both strong and weak scalability. Careful error and timing trade-off analysis are also performed for the cutoff functions induced by the vortex particle method. Our implementation can perform one time step of the velocity+stretching for one billion particles on 32 nodes in 55.9 seconds, which yields 49.12 Tflop/s. © 2012 IEEE.
A task parallel implementation of fast multipole methods
Taura, Kenjiro
2012-11-01
This paper describes a task parallel implementation of ExaFMM, an open source implementation of fast multipole methods (FMM), using a lightweight task parallel library MassiveThreads. Although there have been many attempts on parallelizing FMM, experiences have almost exclusively been limited to formulation based on flat homogeneous parallel loops. FMM in fact contains operations that cannot be readily expressed in such conventional but restrictive models. We show that task parallelism, or parallel recursions in particular, allows us to parallelize all operations of FMM naturally and scalably. Moreover it allows us to parallelize a \\'\\'mutual interaction\\'\\' for force/potential evaluation, which is roughly twice as efficient as a more conventional, unidirectional force/potential evaluation. The net result is an open source FMM that is clearly among the fastest single node implementations, including those on GPUs; with a million particles on a 32 cores Sandy Bridge 2.20GHz node, it completes a single time step including tree construction and force/potential evaluation in 65 milliseconds. The study clearly showcases both programmability and performance benefits of flexible parallel constructs over more monolithic parallel loops. © 2012 IEEE.
Lorenzen, Konstantin; Mathias, Gerald; Tavan, Paul
2015-11-01
Hamiltonian Dielectric Solvent (HADES) is a recent method [S. Bauer et al., J. Chem. Phys. 140, 104103 (2014)] which enables atomistic Hamiltonian molecular dynamics (MD) simulations of peptides and proteins in dielectric solvent continua. Such simulations become rapidly impractical for large proteins, because the computational effort of HADES scales quadratically with the number N of atoms. If one tries to achieve linear scaling by applying a fast multipole method (FMM) to the computation of the HADES electrostatics, the Hamiltonian character (conservation of total energy, linear, and angular momenta) may get lost. Here, we show that the Hamiltonian character of HADES can be almost completely preserved, if the structure-adapted fast multipole method (SAMM) as recently redesigned by Lorenzen et al. [J. Chem. Theory Comput. 10, 3244-3259 (2014)] is suitably extended and is chosen as the FMM module. By this extension, the HADES/SAMM forces become exact gradients of the HADES/SAMM energy. Their translational and rotational invariance then guarantees (within the limits of numerical accuracy) the exact conservation of the linear and angular momenta. Also, the total energy is essentially conserved—up to residual algorithmic noise, which is caused by the periodically repeated SAMM interaction list updates. These updates entail very small temporal discontinuities of the force description, because the employed SAMM approximations represent deliberately balanced compromises between accuracy and efficiency. The energy-gradient corrected version of SAMM can also be applied, of course, to MD simulations of all-atom solvent-solute systems enclosed by periodic boundary conditions. However, as we demonstrate in passing, this choice does not offer any serious advantages.
Hamiltonian Dielectric Solvent (HADES) is a recent method [S. Bauer et al., J. Chem. Phys. 140, 104103 (2014)] which enables atomistic Hamiltonian molecular dynamics (MD) simulations of peptides and proteins in dielectric solvent continua. Such simulations become rapidly impractical for large proteins, because the computational effort of HADES scales quadratically with the number N of atoms. If one tries to achieve linear scaling by applying a fast multipole method (FMM) to the computation of the HADES electrostatics, the Hamiltonian character (conservation of total energy, linear, and angular momenta) may get lost. Here, we show that the Hamiltonian character of HADES can be almost completely preserved, if the structure-adapted fast multipole method (SAMM) as recently redesigned by Lorenzen et al. [J. Chem. Theory Comput. 10, 3244-3259 (2014)] is suitably extended and is chosen as the FMM module. By this extension, the HADES/SAMM forces become exact gradients of the HADES/SAMM energy. Their translational and rotational invariance then guarantees (within the limits of numerical accuracy) the exact conservation of the linear and angular momenta. Also, the total energy is essentially conserved—up to residual algorithmic noise, which is caused by the periodically repeated SAMM interaction list updates. These updates entail very small temporal discontinuities of the force description, because the employed SAMM approximations represent deliberately balanced compromises between accuracy and efficiency. The energy-gradient corrected version of SAMM can also be applied, of course, to MD simulations of all-atom solvent-solute systems enclosed by periodic boundary conditions. However, as we demonstrate in passing, this choice does not offer any serious advantages
Lorenzen, Konstantin; Mathias, Gerald; Tavan, Paul, E-mail: tavan@physik.uni-muenchen.de [Lehrstuhl für BioMolekulare Optik, Ludig–Maximilians Universität München, Oettingenstr. 67, 80538 München (Germany)
2015-11-14
Hamiltonian Dielectric Solvent (HADES) is a recent method [S. Bauer et al., J. Chem. Phys. 140, 104103 (2014)] which enables atomistic Hamiltonian molecular dynamics (MD) simulations of peptides and proteins in dielectric solvent continua. Such simulations become rapidly impractical for large proteins, because the computational effort of HADES scales quadratically with the number N of atoms. If one tries to achieve linear scaling by applying a fast multipole method (FMM) to the computation of the HADES electrostatics, the Hamiltonian character (conservation of total energy, linear, and angular momenta) may get lost. Here, we show that the Hamiltonian character of HADES can be almost completely preserved, if the structure-adapted fast multipole method (SAMM) as recently redesigned by Lorenzen et al. [J. Chem. Theory Comput. 10, 3244-3259 (2014)] is suitably extended and is chosen as the FMM module. By this extension, the HADES/SAMM forces become exact gradients of the HADES/SAMM energy. Their translational and rotational invariance then guarantees (within the limits of numerical accuracy) the exact conservation of the linear and angular momenta. Also, the total energy is essentially conserved—up to residual algorithmic noise, which is caused by the periodically repeated SAMM interaction list updates. These updates entail very small temporal discontinuities of the force description, because the employed SAMM approximations represent deliberately balanced compromises between accuracy and efficiency. The energy-gradient corrected version of SAMM can also be applied, of course, to MD simulations of all-atom solvent-solute systems enclosed by periodic boundary conditions. However, as we demonstrate in passing, this choice does not offer any serious advantages.
Palmesi, Pietro; Bruckner, Florian; Abert, Claas; Suess, Dieter
2016-01-01
The long-range magnetic field is the most time-consuming part in micromagnetic simulations. Improvements both on a numerical and computational basis can relief problems related to this bottleneck. This work presents an efficient implementation of the Fast Multipole Method [FMM] for the magnetic scalar potential as used in micromagnetics. We assume linearly magnetized tetrahedral sources, treat the near field directly and use analytical integration on the multipole expansion in the far field. This approach tackles important issues like the vectorial and continuous nature of the magnetic field. By using FMM the calculations scale linearly in time and memory.
A fast multipole boundary element method for three dimensional potential flow problems
TENG Bin; NING Dezhi; GOU Ying
2004-01-01
A fast multipole methodology (FMM) is developed as a numerical approach to reduce the computational cost and memory requirements in solving large-scale problems. It is applied to the boundary element method (BEM) for threedimensional potential flow problems. The algorithm based on mixed multipole expansion and numerical integration is implemented in combination with an iterative solver. Numerical examinations, on Dirichlet and Neumann problems,are carried out to demonstrate the capability and accuracy of the present method. It has been shown that the method has evident advantages in saving memory and computing time when used to solve huge-scale problems which may be prohibitive for the traditional BEM implementation.
Application of Fast Multipole Methods to the NASA Fast Scattering Code
Dunn, Mark H.; Tinetti, Ana F.
2008-01-01
The NASA Fast Scattering Code (FSC) is a versatile noise prediction program designed to conduct aeroacoustic noise reduction studies. The equivalent source method is used to solve an exterior Helmholtz boundary value problem with an impedance type boundary condition. The solution process in FSC v2.0 requires direct manipulation of a large, dense system of linear equations, limiting the applicability of the code to small scales and/or moderate excitation frequencies. Recent advances in the use of Fast Multipole Methods (FMM) for solving scattering problems, coupled with sparse linear algebra techniques, suggest that a substantial reduction in computer resource utilization over conventional solution approaches can be obtained. Implementation of the single level FMM (SLFMM) and a variant of the Conjugate Gradient Method (CGM) into the FSC is discussed in this paper. The culmination of this effort, FSC v3.0, was used to generate solutions for three configurations of interest. Benchmarking against previously obtained simulations indicate that a twenty-fold reduction in computational memory and up to a four-fold reduction in computer time have been achieved on a single processor.
Schwörer, Magnus; Lorenzen, Konstantin; Mathias, Gerald; Tavan, Paul
2015-03-14
Recently, a novel approach to hybrid quantum mechanics/molecular mechanics (QM/MM) molecular dynamics (MD) simulations has been suggested [Schwörer et al., J. Chem. Phys. 138, 244103 (2013)]. Here, the forces acting on the atoms are calculated by grid-based density functional theory (DFT) for a solute molecule and by a polarizable molecular mechanics (PMM) force field for a large solvent environment composed of several 10(3)-10(5) molecules as negative gradients of a DFT/PMM hybrid Hamiltonian. The electrostatic interactions are efficiently described by a hierarchical fast multipole method (FMM). Adopting recent progress of this FMM technique [Lorenzen et al., J. Chem. Theory Comput. 10, 3244 (2014)], which particularly entails a strictly linear scaling of the computational effort with the system size, and adapting this revised FMM approach to the computation of the interactions between the DFT and PMM fragments of a simulation system, here, we show how one can further enhance the efficiency and accuracy of such DFT/PMM-MD simulations. The resulting gain of total performance, as measured for alanine dipeptide (DFT) embedded in water (PMM) by the product of the gains in efficiency and accuracy, amounts to about one order of magnitude. We also demonstrate that the jointly parallelized implementation of the DFT and PMM-MD parts of the computation enables the efficient use of high-performance computing systems. The associated software is available online. PMID:25770527
The fast multipole method and point dipole moment polarizable force fields
Coles, Jonathan P.; Masella, Michel
2015-01-01
We present an implementation of the fast multipole method for computing Coulombic electrostatic and polarization forces from polarizable force-fields based on induced point dipole moments. We demonstrate the expected O(N) scaling of that approach by performing single energy point calculations on hexamer protein subunits of the mature HIV-1 capsid. We also show the long time energy conservation in molecular dynamics at the nanosecond scale by performing simulations of a protein complex embedded in a coarse-grained solvent using a standard integrator and a multiple time step integrator. Our tests show the applicability of fast multipole method combined with state-of-the-art chemical models in molecular dynamical systems.
The Fast Multipole Method and Point Dipole Moment Polarizable Force Fields
Coles, Jonathan P.; Masella, Michel
2014-01-01
We present an implementation of the fast multipole method for computing coulombic electrostatic and polarization forces from polarizable force-fields based on induced point dipole moments. We demonstrate the expected $O(N)$ scaling of that approach by performing single energy point calculations on hexamer protein subunits of the mature HIV-1 capsid. We also show the long time energy conservation in molecular dynamics at the nanosecond scale by performing simulations of a protein complex embed...
Error control of the vectorial nondirective stable plane wave multilevel fast multipole algorithm
Bogaert, Ignace; Peeters, Joris; De Zutter, Daniël
2011-01-01
Novel formulas are presented that allow the rapid estimation of the number of terms L that needs to be taken into account in the translation operator of the vectorial Nondirective Stable Plane Wave Multilevel Fast Multipole Algorithm (NSPWMLFMA). This is especially important for low frequencies, since the L needed for error-controllability can be substantially higher than the L required in the scalar case. Although these formulas were originally derived for use in the NSPWMLFMA, they are equa...
Ergul, Ozgur
2014-01-01
The Multilevel Fast Multipole Algorithm (MLFMA) for Solving Large-Scale Computational Electromagnetic Problems provides a detailed and instructional overview of implementing MLFMA. The book: Presents a comprehensive treatment of the MLFMA algorithm, including basic linear algebra concepts, recent developments on the parallel computation, and a number of application examplesCovers solutions of electromagnetic problems involving dielectric objects and perfectly-conducting objectsDiscusses applications including scattering from airborne targets, scattering from red
White, J.; Phillips, J.R.; Korsmeyer, T. [Massachusetts Institute of Technology, Cambridge, MA (United States)
1994-12-31
Mixed first- and second-kind surface integral equations with (1/r) and {partial_derivative}/{partial_derivative} (1/r) kernels are generated by a variety of three-dimensional engineering problems. For such problems, Nystroem type algorithms can not be used directly, but an expansion for the unknown, rather than for the entire integrand, can be assumed and the product of the singular kernal and the unknown integrated analytically. Combining such an approach with a Galerkin or collocation scheme for computing the expansion coefficients is a general approach, but generates dense matrix problems. Recently developed fast algorithms for solving these dense matrix problems have been based on multipole-accelerated iterative methods, in which the fast multipole algorithm is used to rapidly compute the matrix-vector products in a Krylov-subspace based iterative method. Another approach to rapidly computing the dense matrix-vector products associated with discretized integral equations follows more along the lines of a multigrid algorithm, and involves projecting the surface unknowns onto a regular grid, then computing using the grid, and finally interpolating the results from the regular grid back to the surfaces. Here, the authors describe a precorrectted-FFT approach which can replace the fast multipole algorithm for accelerating the dense matrix-vector product associated with discretized potential integral equations. The precorrected-FFT method, described below, is an order n log(n) algorithm, and is asymptotically slower than the order n fast multipole algorithm. However, initial experimental results indicate the method may have a significant constant factor advantage for a variety of engineering problems.
An efficient blocking M2L translation for low-frequency fast multipole method in three dimensions
Takahashi, Toru; Shimba, Yuta; Isakari, Hiroshi; Matsumoto, Toshiro
2016-05-01
We propose an efficient scheme to perform the multipole-to-local (M2L) translation in the three-dimensional low-frequency fast multipole method (LFFMM). Our strategy is to combine a group of matrix-vector products associated with M2L translation into a matrix-matrix product in order to diminish the memory traffic. For this purpose, we first developed a grouping method (termed as internal blocking) based on the congruent transformations (rotational and reflectional symmetries) of M2L-translators for each target box in the FMM hierarchy (adaptive octree). Next, we considered another method of grouping (termed as external blocking) that was able to handle M2L translations for multiple target boxes collectively by using the translational invariance of the M2L translation. By combining these internal and external blockings, the M2L translation can be performed efficiently whilst preservingthe numerical accuracy exactly. We assessed the proposed blocking scheme numerically and applied it to the boundary integral equation method to solve electromagnetic scattering problems for perfectly electrical conductor. From the numerical results, it was found that the proposed M2L scheme achieved a few times speedup compared to the non-blocking scheme.
Fast multipole acceleration of the MEG/EEG boundary element method
The accurate solution of the forward electrostatic problem is an essential first step before solving the inverse problem of magneto- and electroencephalography (MEG/EEG). The symmetric Galerkin boundary element method is accurate but cannot be used for very large problems because of its computational complexity and memory requirements. We describe a fast multipole-based acceleration for the symmetric boundary element method (BEM). It creates a hierarchical structure of the elements and approximates far interactions using spherical harmonics expansions. The accelerated method is shown to be as accurate as the direct method, yet for large problems it is both faster and more economical in terms of memory consumption
Periodic Boundary Conditions and the Error-Controlled Fast Multipole Method
Kabadshow, I.
2012-01-01
Die Simulation von paarweisen Wechselwirkungen für sehr große Teilchen-Ensembles ist eine wesentliche Aufgabe in der wissenschaftlichen Forschung. Besonders die Berechnung von langreichweitigen Wechselwirkungen legt der Systemgröße Beschränkungen auf, da die Anzahl der zu berechnenden Wechselwirkungen quadratisch mit der Anzahl der Teilchen skaliert. Schnelle Summationsverfahren wie die Fast Multipole Method (FMM) können dazu beitragen die Komplexität auf O(N) zu reduzieren. Diese Arbeit erwe...
Application of A Fast Multipole BIEM for Flow Diffraction from A 3D Body
滕斌; 宁德志
2004-01-01
A Fast Multipole Method (FMM) is developed as a numerical approach to the reduction of the computational cost and requirement memory capacity for a large in solving large-scale problems. In this paper it is applied to the boundary integral equation method (BIEM) for current diffraction from arbitrary 3D bodies. The boundary integral equation is discretized by higher order elements, the FMM is applied to avoid the matrix/vector product, and the resulting algebraic equation is solved by the Generalized Conjugate Residual method (GCR). Numerical examination shows that the FMM is more efficient than the direct evaluation method in computational cost and storage of computers.
Fast multipole boundary element analysis of 2D viscoelastic composites with imperfect interfaces
无
2010-01-01
A fast multipole boundary element method(FMBEM)is developed for the analysis of 2D linear viscoelastic composites with imperfect viscoelastic interfaces.The transformed fast multipole formulations are established using the time domain method. To simulate the viscoelastic behavior of imperfect interfaces that are frequently encountered in practice,the Kelvin type model is introduced.The FMBEM is further improved by incorporating naturally the interaction among inclusions as well as eliminating the phenomenon of material penetration.Since all the integrals are evaluated analytically,high accuracy and fast convergence of the numerical scheme are obtained.Several numerical examples,including planar viscoelastic composites with a single inclusion or randomly distributed multi-inclusions are presented.The numerical results are compared with the developed analytical solutions,which illustrates that the proposed FMBEM is very efficient in determining the macroscopic viscoelastic behavior of the particle-reinforced composites with the presence of imperfect interfaces.The laboratory measurements of the mixture creep compliance of asphalt concrete are also compared with the prediction by the developed model.
Nitadori, Keigo
2014-01-01
We propose an efficient algorithm for the evaluation of the potential and its gradient of gravitational/electrostatic $N$-body systems, which we call particle mesh multipole method (PMMM or PM$^3$). PMMM can be understood both as an extension of the particle mesh (PM) method and as an optimization of the fast multipole method (FMM).In the former viewpoint, the scalar density and potential held by a grid point are extended to multipole moments and local expansions in $(p+1)^2$ real numbers, where $p$ is the order of expansion. In the latter viewpoint, a hierarchical octree structure which brings its $\\mathcal O(N)$ nature, is replaced with a uniform mesh structure, and we exploit the convolution theorem with fast Fourier transform (FFT) to speed up the calculations. Hence, independent $(p+1)^2$ FFTs with the size equal to the number of grid points are performed. The fundamental idea is common to PPPM/MPE by Shimada et al. (1993) and FFTM by Ong et al. (2003). PMMM differs from them in supporting both the open ...
A Fast Multipole Algorithm with Virtual Cube Partitioning for 3-D Capacitance Extraction
YANGZhaozhi; WANGZeyi
2004-01-01
In this paper a fast indirect boundaryelement method based on the multipole algorithm for capacitance extraction of three-dimensional (3-D) geometries, virtual cube multipole algorithm, is described. First,each 2-D boundary element is regarded as a set of particles with charge rather than a single particle, so the relations between the positions of elements themselves are considered instead of the relations between the center-points of the elements, and a new strategy for cube partitioning is introduced. This strategy overcomes the inadequacy of the methods that associating panels to particles, does not need to break up every panel contained in more than one cube, and has higher speed and precision. Next, a new method is proposed to accelerate the potential integration between the panels that are near to each other. Making good use of the similarity in the 2-D boundary integration,the fast potential integral approach decreases the burden of direct potential computing. Experiments confirm that the algorithm is accurate and has nearly linear computational growth as O(nm), where n is the number of panels and rn is the number of conductors. The new algorithm is implemented and the performance is compared with previous algorithms, such as Fastcap2 of MIT, for k×k bus examples.
Periodic boundary conditions and the error-controlled fast multipole method
Kabadshow, Ivo
2012-08-22
The simulation of pairwise interactions in huge particle ensembles is a vital issue in scientific research. Especially the calculation of long-range interactions poses limitations to the system size, since these interactions scale quadratically with the number of particles. Fast summation techniques like the Fast Multipole Method (FMM) can help to reduce the complexity to O(N). This work extends the possible range of applications of the FMM to periodic systems in one, two and three dimensions with one unique approach. Together with a tight error control, this contribution enables the simulation of periodic particle systems for different applications without the need to know and tune the FMM specific parameters. The implemented error control scheme automatically optimizes the parameters to obtain an approximation for the minimal runtime for a given energy error bound.
Efficient analysis of dielectric radomes using multilevel fast multipole algorithm with CRWG basis
无
2008-01-01
A full-wave analysis of the electromagnetic problem of a three-dimensional (3-D) antenna radiating through a 3-D dielectric radome is preserued.The problem is formulated using the Poggio-Miller-Chang-Harrington-Wu(PMCHW) approach for homogeneous dielectric objects and the electric field integral equation for conducting objects.The integral equations are discretized by the method of moment (MoM),in which the conducting and dielectric surface/interfaces are represented by curvilinear triangular patches and the unknown equivalent electric and magnetic currents are expanded using curvilinear RWG basis functions.The resultant matrix equation is then solved by the multilevel fast multipole algorithm (MLFMA) and fast far-field approximation (FAFFA) is used to further accelerate the computation.The radiation patterns of dipole arrays in the presence of radomes are presented.The numerical results demonstrate the accuracy and versatility of this method.
Computing dislocation stress fields in anisotropic elastic media using fast multipole expansions
The calculation of stress fields due to dislocations and hence the forces they exert on each other is the most time consuming step in dislocation dynamics (DD) simulations. The fast multipole method (FMM) can reduce the computational cost at each simulation step from O(N2) to O(N) for an ensemble of N dislocation segments. However, FMM has not yet been applied to three-dimensional DD simulations which take into account anisotropic elasticity. We demonstrate a systematic procedure to establish this capability by first obtaining the derivatives of the elastic Green's function to arbitrary order for a medium of general anisotropy. We then compute the stress field of a dislocation loop using multipole expansions based on these derivatives, and analyze the dependence of numerical errors on the expansion order. This method can be implemented in large scale DD simulations when the consideration of elastic anisotropy is necessary, for example the technologically important cases of iron and ferritic steels at high temperatures. (paper)
Another preprocessing algorithm for generalized one-dimensional fast multipole method
The fast multipole method (FMM), which is originally an algorithm for fast evaluation of particle interactions, is also effective for accelerating several numerical computations. Yarvin and Rokhlin proposed 'generalized' FMM using the singular value decomposition (SVD), which gives the optimum low-rank approximation. Their algorithm reduces the computational costs of the FMM evaluation and frees the FMM from analytical approximation formulae. However, the computational complexity of the preprocessing for an NxN matrix is O(N3) because of the SVD, and it requires orthogonal matrices of the low-rank approximations. In this paper we propose another preprocessing algorithm for the generalized FMM. Our algorithm runs in time O(N2) even with the SVD and releases the low-rank approximations from orthogonal matrices. The triangularization by the QR decomposition with sparsification, which reduces the costs of the FMM more than the diagonalization, is enabled. Although the algorithm by Yarvin and Rokhlin can be accelerated to O(N2) using the QR decomposition, our preprocessing algorithm outperforms it in fast spherical filter, fast polynomial interpolation and fast Legendre transform
A Tuned and Scalable Fast Multipole Method as a Preeminent Algorithm for Exascale Systems
Yokota, Rio
2011-01-01
Achieving computing at the exascale means accelerating today's applications by one thousand times. Clearly, this cannot be accomplished by hardware alone, at least not in the short time frame expected for reaching this performance milestone. Thus, a lively discussion has begun in the last couple of years about programming models, software components and tools, and algorithms that will facilitate exascale computing. Among the algorithms that are likely to play a preeminent role in the new world of computing, the fast multipole method (F MM) appears as a rising star. Due to its hierarchical nature and the techniques used to access the data via a tree structure, it is not a locality-sensitive application. It also enjoys favorable synchronization patterns, again, of a hierarchical nature, where many operations can happen simultaneously at each level of the hierarchy. In this paper, we present a discussion of the features of the F MM that make it a particularly favorable algorithm for the emerging heterogeneous, m...
The Fast Multipole Method and Point Dipole Moment Polarizable Force Fields
Coles, Jonathan P
2014-01-01
We present a momentum conserving implementation of the fast multipole method for computing coulombic electrostatic and polarization forces from polarizable force-fields based on induced point dipole moments. We demonstrate the expected $O(N)$ scaling of that approach by performing single energy point calculations on hexamer protein subunits of the mature HIV-1 capsid. We also show the long time energy conservation in molecular dynamics at the nanosecond scale by performing simulations of a protein complex embedded in a coarse-grained solvent using both a standard integrator and a multiple time step one. Our tests show the applicability of FMM combined with state-of-the-art chemical models in molecular dynamical systems.
Pan Xiaomin; Sheng Xinqing
2008-01-01
A general and efficient parallel approach is proposed for the first time to parallelize the hybrid finite-element-boundary-integral-multi-level fast multipole algorithm (FE-BI-MLFMA). Among many algorithms of FE-BI-MLFMA, the decomposition algorithm (DA) is chosen as a basis for the parallelization of FE-BI-MLFMA because of its distinct numerical characteristics suitable for parallelization. On the basis of the DA, the parallelization of FE-BI-MLFMA is carried out by employing the parallelized multi-frontal method for the matrix from the finite-element method and the parallelized MLFMA for the matrix from the boundary integral method respectively. The programming and numerical experiments of the proposed parallel approach are carried out in the high perfor-mance computing platform CEMS-Liuhui. Numerical experiments demonstrate that FE-BI-MLFMA is efficiently parallelized and its computational capacity is greatly improved without losing accuracy, efficiency, and generality.
Toivanen, Elias A; Losilla, Sergio A; Sundholm, Dage
2015-12-21
Algorithms and working expressions for a grid-based fast multipole method (GB-FMM) have been developed and implemented. The computational domain is divided into cubic subdomains, organized in a hierarchical tree. The contribution to the electrostatic interaction energies from pairs of neighboring subdomains is computed using numerical integration, whereas the contributions from further apart subdomains are obtained using multipole expansions. The multipole moments of the subdomains are obtained by numerical integration. Linear scaling is achieved by translating and summing the multipoles according to the tree structure, such that each subdomain interacts with a number of subdomains that are almost independent of the size of the system. To compute electrostatic interaction energies of neighboring subdomains, we employ an algorithm which performs efficiently on general purpose graphics processing units (GPGPU). Calculations using one CPU for the FMM part and 20 GPGPUs consisting of tens of thousands of execution threads for the numerical integration algorithm show the scalability and parallel performance of the scheme. For calculations on systems consisting of Gaussian functions (α = 1) distributed as fullerenes from C20 to C720, the total computation time and relative accuracy (ppb) are independent of the system size. PMID:26006111
Jakobus, U.; J. Van Tonder
2005-01-01
An existing method of moments (MoM) code for the solution of complex scattering bodies has been accelerated by means of a multilevel fast multipole method (MLFMM). We demonstrate the usage of this technique both for metallic structures (wires and surfaces) and for dielectric bodies (volume and surface equivalence principle). Aspects like the effect of the type of integral equation, preconditioning schemes, or iterative solution techniques are discussed. But also limitations are addressed, whi...
Ergül, Özgür; Gürel, Levent
2013-03-01
Accurate electromagnetic modeling of complicated optical structures poses several challenges. Optical metamaterial and plasmonic structures are composed of multiple coexisting dielectric and/or conducting parts. Such composite structures may possess diverse values of conductivities and dielectric constants, including negative permittivity and permeability. Further challenges are the large sizes of the structures with respect to wavelength and the complexities of the geometries. In order to overcome these challenges and to achieve rigorous and efficient electromagnetic modeling of three-dimensional optical composite structures, we have developed a parallel implementation of the multilevel fast multipole algorithm (MLFMA). Precise formulation of composite structures is achieved with the so-called "electric and magnetic current combined-field integral equation." Surface integral equations are carefully discretized with piecewise linear basis functions, and the ensuing dense matrix equations are solved iteratively with parallel MLFMA. The hierarchical strategy is used for the efficient parallelization of MLFMA on distributed-memory architectures. In this paper, fast and accurate solutions of large-scale canonical and complicated real-life problems, such as optical metamaterials, discretized with tens of millions of unknowns are presented in order to demonstrate the capabilities of the proposed electromagnetic solver. PMID:23456127
Gumerov, Nail A; Duraiswami, Ramani
2009-01-01
The development of a fast multipole method (FMM) accelerated iterative solution of the boundary element method (BEM) for the Helmholtz equations in three dimensions is described. The FMM for the Helmholtz equation is significantly different for problems with low and high kD (where k is the wavenumber and D the domain size), and for large problems the method must be switched between levels of the hierarchy. The BEM requires several approximate computations (numerical quadrature, approximations of the boundary shapes using elements), and these errors must be balanced against approximations introduced by the FMM and the convergence criterion for iterative solution. These different errors must all be chosen in a way that, on the one hand, excess work is not done and, on the other, that the error achieved by the overall computation is acceptable. Details of translation operators for low and high kD, choice of representations, and BEM quadrature schemes, all consistent with these approximations, are described. A novel preconditioner using a low accuracy FMM accelerated solver as a right preconditioner is also described. Results of the developed solvers for large boundary value problems with 0.0001 less, similarkD less, similar500 are presented and shown to perform close to theoretical expectations. PMID:19173406
Petascale molecular dynamics simulation using the fast multipole method on K computer
Ohno, Yousuke
2014-10-01
In this paper, we report all-atom simulations of molecular crowding - a result from the full node simulation on the "K computer", which is a 10-PFLOPS supercomputer in Japan. The capability of this machine enables us to perform simulation of crowded cellular environments, which are more realistic compared to conventional MD simulations where proteins are simulated in isolation. Living cells are "crowded" because macromolecules comprise ∼30% of their molecular weight. Recently, the effects of crowded cellular environments on protein stability have been revealed through in-cell NMR spectroscopy. To measure the performance of the "K computer", we performed all-atom classical molecular dynamics simulations of two systems: target proteins in a solvent, and target proteins in an environment of molecular crowders that mimic the conditions of a living cell. Using the full system, we achieved 4.4 PFLOPS during a 520 million-atom simulation with cutoff of 28 Å. Furthermore, we discuss the performance and scaling of fast multipole methods for molecular dynamics simulations on the "K computer", as well as comparisons with Ewald summation methods. © 2014 Elsevier B.V. All rights reserved.
Combining the multilevel fast multipole method with the uniform geometrical theory of diffraction
A. Tzoulis
2005-01-01
Full Text Available The presence of arbitrarily shaped and electrically large objects in the same environment leads to hybridization of the Method of Moments (MoM with the Uniform Geometrical Theory of Diffraction (UTD. The computation and memory complexity of the MoM solution is improved with the Multilevel Fast Multipole Method (MLFMM. By expanding the k-space integrals in spherical harmonics, further considerable amount of memory can be saved without compromising accuracy and numerical speed. However, until now MoM-UTD hybrid methods are restricted to conventional MoM formulations only with Electric Field Integral Equation (EFIE. In this contribution, a MLFMM-UTD hybridization for Combined Field Integral Equation (CFIE is proposed and applied within a hybrid Finite Element - Boundary Integral (FEBI technique. The MLFMM-UTD hybridization is performed at the translation procedure on the various levels of the MLFMM, using a far-field approximation of the corresponding translation operator. The formulation of this new hybrid technique is presented, as well as numerical results.
Zhang, He
2013-01-01
The space charge effect is one of the most important collective effects in beam dynamic studies. In many cases, numerical simulations are inevitable in order to get a clear understanding of this effect. The particle-particle interaction algorithms and the article-in-cell algorithms are widely used in space charge effect simulations. But they both have difficulties in dealing with highly correlated beams with abnormal distributions or complicated geometries. We developed a new algorithm to calculate the three dimensional self-field between charged particles by combining the differential algebra (DA) techniques with the fast multi-pole method (FMM). The FMM hierarchically decomposes the whole charged domain into many small regions. For each region it uses multipole expansions to represent the potential/field contributions from the particles far away from the region and then converts the multipole expansions into a local expansion inside the region. The potential/field due to the far away particles is calculated from the expansions and the potential/field due to the nearby particles is calculated from the Coulomb force law. The DA techniques are used in the calculation, translation and converting of the expansions. The new algorithm scales linearly with the total number of particles and it is suitable for any arbitrary charge distribution. Using the DA techniques, we can calculate both the potential/field and its high order derivatives, which will be useful for the purpose of including the space charge effect into transfer maps in the future. We first present the single level FMM, which decomposes the whole domain into boxes of the same size. It works best for charge distributions that are not overly non-uniform. Then we present the multilevel fast multipole algorithm (MLFMA), which decomposes the whole domain into different sized boxes according to the charge density. Finer boxes are generated where the higher charge density exists; thus the algorithm works for any
Inexact Krylov iterations and relaxation strategies with fast-multipole boundary element method
Layton, Simon K.; Barba, Lorena A.
2015-01-01
Boundary element methods produce dense linear systems that can be accelerated via multipole expansions. Solved with Krylov methods, this implies computing the matrix-vector products within each iteration with some error, at an accuracy controlled by the order of the expansion, $p$. We take advantage of a unique property of Krylov iterations that allow lower accuracy of the matrix-vector products as convergence proceeds, and propose a relaxation strategy based on progressively decreasing $p$. ...
Fajardo, Kristel C Meza; Chaillat, Stéphanie; Lenti, Luca
2016-01-01
In this work, we study seismic wave amplification in alluvial basins having 3D standard geometries through the Fast Multipole Boundary Element Method in the frequency domain. We investigate how much 3D amplification differs from the 1D (horizontal layering) case. Considering incident fields of plane harmonic waves, we examine the relationships between the amplification level and the most relevant physical parameters of the problem (impedance contrast, 3D aspect ratio, vertical and oblique incidence of plane waves). The FMBEM results show that the most important parameters for wave amplification are the impedance contrast and the so-called equivalent shape ratio. Using these two parameters, we derive simple rules to compute the fundamental frequency for various 3D basin shapes and the corresponding 3D/1D amplification factor for 5% damping. Effects on amplification due to 3D basin asymmetry are also studied and incorporated in the derived rules.
Poursina, Mohammad; Anderson, Kurt S.
2014-08-01
This paper presents a novel algorithm to approximate the long-range electrostatic potential field in the Cartesian coordinates applicable to 3D coarse-grained simulations of biopolymers. In such models, coarse-grained clusters are formed via treating groups of atoms as rigid and/or flexible bodies connected together via kinematic joints. Therefore, multibody dynamic techniques are used to form and solve the equations of motion of such coarse-grained systems. In this article, the approximations for the potential fields due to the interaction between a highly negatively/positively charged pseudo-atom and charged particles, as well as the interaction between clusters of charged particles, are presented. These approximations are expressed in terms of physical and geometrical properties of the bodies such as the entire charge, the location of the center of charge, and the pseudo-inertia tensor about the center of charge of the clusters. Further, a novel substructuring scheme is introduced to implement the presented far-field potential evaluations in a binary tree framework as opposed to the existing quadtree and octree strategies of implementing fast multipole method. Using the presented Lagrangian grids, the electrostatic potential is recursively calculated via sweeping two passes: assembly and disassembly. In the assembly pass, adjacent charged bodies are combined together to form new clusters. Then, the potential field of each cluster due to its interaction with faraway resulting clusters is recursively calculated in the disassembly pass. The method is highly compatible with multibody dynamic schemes to model coarse-grained biopolymers. Since the proposed method takes advantage of constant physical and geometrical properties of rigid clusters, improvement in the overall computational cost is observed comparing to the tradition application of fast multipole method.
This paper presents a novel algorithm to approximate the long-range electrostatic potential field in the Cartesian coordinates applicable to 3D coarse-grained simulations of biopolymers. In such models, coarse-grained clusters are formed via treating groups of atoms as rigid and/or flexible bodies connected together via kinematic joints. Therefore, multibody dynamic techniques are used to form and solve the equations of motion of such coarse-grained systems. In this article, the approximations for the potential fields due to the interaction between a highly negatively/positively charged pseudo-atom and charged particles, as well as the interaction between clusters of charged particles, are presented. These approximations are expressed in terms of physical and geometrical properties of the bodies such as the entire charge, the location of the center of charge, and the pseudo-inertia tensor about the center of charge of the clusters. Further, a novel substructuring scheme is introduced to implement the presented far-field potential evaluations in a binary tree framework as opposed to the existing quadtree and octree strategies of implementing fast multipole method. Using the presented Lagrangian grids, the electrostatic potential is recursively calculated via sweeping two passes: assembly and disassembly. In the assembly pass, adjacent charged bodies are combined together to form new clusters. Then, the potential field of each cluster due to its interaction with faraway resulting clusters is recursively calculated in the disassembly pass. The method is highly compatible with multibody dynamic schemes to model coarse-grained biopolymers. Since the proposed method takes advantage of constant physical and geometrical properties of rigid clusters, improvement in the overall computational cost is observed comparing to the tradition application of fast multipole method
Inexact Krylov iterations and relaxation strategies with fast-multipole boundary element method
Layton, Simon K
2015-01-01
Boundary element methods produce dense linear systems that can be accelerated via multipole expansions. Solved with Krylov methods, this implies computing the matrix-vector products within each iteration with some error, at an accuracy controlled by the order of the expansion, $p$. We take advantage of a unique property of Krylov iterations that allow lower accuracy of the matrix-vector products as convergence proceeds, and propose a relaxation strategy based on progressively decreasing $p$. Via extensive numerical tests, we show that the relaxed Krylov iterations converge with speed-ups of between 2x and 4x for Laplace problems and between 3.5x and 4.5x for Stokes problems. We include an application to Stokes flow around red blood cells, computing with up to 64 cells and problem size up to 131k boundary elements and nearly 400k unknowns. The study was done with an in-house multi-threaded C++ code, on a quad-core CPU.
Zheng, Chang-Jun; Gao, Hai-Feng; Du, Lei; Chen, Hai-Bo; Zhang, Chuanzeng
2016-01-01
An accurate numerical solver is developed in this paper for eigenproblems governed by the Helmholtz equation and formulated through the boundary element method. A contour integral method is used to convert the nonlinear eigenproblem into an ordinary eigenproblem, so that eigenvalues can be extracted accurately by solving a set of standard boundary element systems of equations. In order to accelerate the solution procedure, the parameters affecting the accuracy and efficiency of the method are studied and two contour paths are compared. Moreover, a wideband fast multipole method is implemented with a block IDR (s) solver to reduce the overall solution cost of the boundary element systems of equations with multiple right-hand sides. The Burton-Miller formulation is employed to identify the fictitious eigenfrequencies of the interior acoustic problems with multiply connected domains. The actual effect of the Burton-Miller formulation on tackling the fictitious eigenfrequency problem is investigated and the optimal choice of the coupling parameter as α = i / k is confirmed through exterior sphere examples. Furthermore, the numerical eigenvalues obtained by the developed method are compared with the results obtained by the finite element method to show the accuracy and efficiency of the developed method.
Quantum fast-roll initial conditions for the inflaton which are different from the classical fast-roll conditions and from the quantum slow-roll conditions can lead to inflation that lasts long enough. These quantum fast-roll initial conditions for the inflaton allow for kinetic energies of the same order of the potential energies and nonperturbative inflaton modes with nonzero wave numbers. Their evolution starts with a transitory epoch where the redshift due to the expansion succeeds to assemble the quantum excited modes of the inflaton in a homogeneous (zero mode) condensate, and the large value of the Hubble parameter succeeds to overdamp the fast roll of the redshifted inflaton modes. After this transitory stage the effective classical slow-roll epoch is reached. Most of the e-folds are produced during the slow-roll epoch, and we recover the classical slow-roll results for the scalar and tensor metric perturbations plus corrections. These corrections are important if scales which are horizon size today exited the horizon by the end of the transitory stage and, as a consequence, the lower cosmic microwave background (CMB) multipoles get suppressed or enhanced. Both for scalar and tensor metric perturbations, fast roll leads to a suppression of the amplitude of the perturbations (and of the low CMB multipoles), while the quantum precondensate epoch gives an enhancement of the amplitude of the perturbations (and of the low CMB multipoles). These two types of corrections can compete and combine in a scale dependent manner. They turn out to be smaller in new inflation than in chaotic inflation. These corrections arise as natural consequences of the quantum nonperturbative inflaton dynamics, and can allow a further improvement of the fitting of inflation plus the ΛCMB model to the observed CMB spectra. In addition, the corrections to the tensor metric perturbations will provide an independent test of this model. Thus, the effects of quantum inflaton fast
To simulate numerically a non-destructive by eddy current testing (NDT-CF), the sensor response can be modeled through a semi-analytical approach by volume integral equations. Faster than the finite element method, this approach is however restricted to the study of plane or cylindrical parts (without taking into account the edge effects) because of the complexity of the expression of the dyadic Green function for more general configurations. However, there is an industrial demand to extend the capabilities of the CF model in complex configurations (deformed plates, edges effects...). We were thus brought to formulate the electromagnetic problem differently, by setting ourselves the goal of maintaining a semi-analytical approach. The surface integral equation (SIE) expresses the volume problem by an equivalent transmission one at the interfaces (2D) between homogeneous sub-domains. This problem is approached by a linear system (by the method of moments), whose number of unknowns is reduced due to the nature of the surfacic mesh. Therefore, this system can be solved by a direct solver for small configurations. That enabled us to treat several various positions of the sensor for only one inversion of the impedance matrix. The numerical results obtained using this formulation involve plates with consideration of edge effects such as edge and corner. They are consistent with results obtained by the finite element method. For larger configurations, we conducted a preliminary study for the adaptation of an acceleration method of the matrix vector product involved in an iterative solver (fast multipole method or FMM) to define the conditions under which the FMM calculation works correctly (accuracy, convergence...) in the NDT's domain. A special attention has been given to the choice of basis functions (which have to satisfy an Hdiv conforming property) and on the evaluation of near interactions (which are weakly singular). (author)
Motion Tracking with Fast Adaptive Background Subtraction
Xiao De-gui; Yu Sheng-sheng; Zhou Jing-li
2003-01-01
To extract and track moving objects is usually one of the most important tasks of intelligent video surveillance systems. This paper presents a fast and adaptive background subtraction algorithm and the motion tracking process using this algorithm. The algorithm uses only luminance components of sampled image sequence pixels and models every pixel in a statistical model. The algorithm is characterized by its ability of real time detecting sudden lighting changes, and extracting and tracking motion objects faster. It is shown that our algorithm can be realized with lower time and space complexity and adjustable object detection error rate with comparison to other background subtraction algorithms. Making use of the algorithm, an indoor monitoring system is also worked out and the motion tracking process is presented in this paper. Experimental results testify the algorithm's good performances when used in an indoor monitoring system.
Jiang, Xikai; Zhao, Xujun; Qin, Jian; Karpeev, Dmitry; Hernandez-Ortiz, Juan; de Pablo, Juan; Heinonen, Olle
2016-01-01
Large classes of materials systems in physics and engineering are governed by magnetic and electrostatic interactions. Continuum or mesoscale descriptions of such systems can be cast in terms of integral equations, whose direct computational evaluation requires O(N^2) operations, where N is the number of unknowns. Such a scaling, which arises from the many-body nature of the relevant Green's function, has precluded wide-spread adoption of integral methods for solution of large-scale scientific and engineering problems. In this work, a parallel computational approach is presented that relies on using scalable open source libraries and utilizes a kernel-independent Fast Multipole Method to evaluate the integrals in O(N) operations, with O(N) memory cost, thereby substantially improving the scalability and efficiency of computational integral methods. We demonstrate the accuracy, efficiency, and scalability of our approach in the contest of two examples. In the first, we solve a boundary value problem for a ferr...
Adaptive grouping for the higher-order multilevel fast multipole method
Borries, Oscar Peter; Jørgensen, Erik; Meincke, Peter; Hansen, Per Christian
2014-01-01
functions. Results from both a uniformly and nonuniformly meshed scatterer are presented, showing how the technique is worthwhile even for regular meshes, and demonstrating that there is no loss of accuracy in spite of the large reduction in memory requirements and the relatively low computational cost....
Lee, Dongryeol; Gray, Alexander G
2011-01-01
A three-body potential function can account for interactions among triples of particles which are uncaptured by pairwise interaction functions such as Coulombic or Lennard-Jones potentials. Likewise, a multibody potential of order $n$ can account for interactions among $n$-tuples of particles uncaptured by interaction functions of lower orders. To date, the computation of multibody potential functions for a large number of particles has not been possible due to its $O(N^n)$ scaling cost. In this paper we describe a fast tree-code for efficiently approximating multibody potentials. For the first time, we show how to extend the series-expansion-based approach of Fast Multipole Method-like algorithms to handle interactions among more than two particles. Our approach guarantees a user-specified bound on the absolute or relative error in the computed potential. We provide speedup results on a three-body dispersion potential, the Axilrod-Teller potential.
Jiang, Xikai; Li, Jiyuan; Zhao, Xujun; Qin, Jian; Karpeev, Dmitry; Hernandez-Ortiz, Juan; de Pablo, Juan J.; Heinonen, Olle
2016-08-01
Large classes of materials systems in physics and engineering are governed by magnetic and electrostatic interactions. Continuum or mesoscale descriptions of such systems can be cast in terms of integral equations, whose direct computational evaluation requires O(N2) operations, where N is the number of unknowns. Such a scaling, which arises from the many-body nature of the relevant Green's function, has precluded wide-spread adoption of integral methods for solution of large-scale scientific and engineering problems. In this work, a parallel computational approach is presented that relies on using scalable open source libraries and utilizes a kernel-independent Fast Multipole Method (FMM) to evaluate the integrals in O(N) operations, with O(N) memory cost, thereby substantially improving the scalability and efficiency of computational integral methods. We demonstrate the accuracy, efficiency, and scalability of our approach in the context of two examples. In the first, we solve a boundary value problem for a ferroelectric/ferromagnetic volume in free space. In the second, we solve an electrostatic problem involving polarizable dielectric bodies in an unbounded dielectric medium. The results from these test cases show that our proposed parallel approach, which is built on a kernel-independent FMM, can enable highly efficient and accurate simulations and allow for considerable flexibility in a broad range of applications.
Fast-start adaptation funding: keeping promises from Copenhagen
Ciplet, David; Roberts, J. Timmons [Brown University (United States); Chandani, Achala; Huq, Saleemul
2010-11-15
The most concrete commitment to come out of the international climate negotiations in Copenhagen was US$30 billion dollars in 'fast-start climate finance' to developing countries, with balanced support for both mitigation and adaptation. Fast-start adaptation finance, in particular, is crusial for poor countries facing rapid climate change. But so far, pledges for adaptation from developed countries have been inadequate and unclear. This briefing outlines ways for the Cancun negotiations to address five crucial issues for adaptation finance: (1) the amount and type of funding being offered, (2) the definition of adaptation, (3) global oversight and accounting, (4) a clear baseline and transparent spending and (5) the channel for delivering funds.
Sox17 regulates liver lipid metabolism and adaptation to fasting.
Rommelaere, Samuel; Millet, Virginie; Vu Manh, Thien-Phong; Gensollen, Thomas; Andreoletti, Pierre; Cherkaoui-Malki, Mustapha; Bourges, Christophe; Escalière, Bertrand; Du, Xin; Xia, Yu; Imbert, Jean; Beutler, Bruce; Kanai, Yoshiakira; Malissen, Bernard; Malissen, Marie; Tailleux, Anne; Staels, Bart; Galland, Franck; Naquet, Philippe
2014-01-01
Liver is a major regulator of lipid metabolism and adaptation to fasting, a process involving PPARalpha activation. We recently showed that the Vnn1 gene is a PPARalpha target gene in liver and that release of the Vanin-1 pantetheinase in serum is a biomarker of PPARalpha activation. Here we set up a screen to identify new regulators of adaptation to fasting using the serum Vanin-1 as a marker of PPARalpha activation. Mutagenized mice were screened for low serum Vanin-1 expression. Functional interactions with PPARalpha were investigated by combining transcriptomic, biochemical and metabolic approaches. We characterized a new mutant mouse in which hepatic and serum expression of Vanin-1 is depressed. This mouse carries a mutation in the HMG domain of the Sox17 transcription factor. Mutant mice display a metabolic phenotype featuring lipid abnormalities and inefficient adaptation to fasting. Upon fasting, a fraction of the PPARα-driven transcriptional program is no longer induced and associated with impaired fatty acid oxidation. The transcriptional phenotype is partially observed in heterozygous Sox17+/- mice. In mutant mice, the fasting phenotype but not all transcriptomic signature is rescued by the administration of the PPARalpha agonist fenofibrate. These results identify a novel role for Sox17 in adult liver as a modulator of the metabolic adaptation to fasting. PMID:25141153
Fast Histograms using Adaptive CUDA Streams
Koppaka, Sisir; Narasimhan, Srihari; Narayanan, Babu
2010-01-01
Histograms are widely used in medical imaging, network intrusion detection, packet analysis and other stream-based high throughput applications. However, while porting such software stacks to the GPU, the computation of the histogram is a typical bottleneck primarily due to the large impact on kernel speed by atomic operations. In this work, we propose a stream-based model implemented in CUDA, using a new adaptive kernel that can be optimized based on latency hidden CPU compute. We also explore the tradeoffs of using the new kernel vis-\\`a-vis the stock NVIDIA SDK kernel, and discuss an intelligent kernel switching method for the stream based on a degeneracy criterion that is adaptively computed from the input stream.
Information criterion based fast PCA adaptive algorithm
Li Jiawen; Li Congxin
2007-01-01
The novel information criterion (NIC) algorithm can find the principal subspace quickly, but it is not an actual principal component analysis (PCA) algorithm and hence it cannot find the orthonormal eigen-space which corresponds to the principal component of input vector.This defect limits its application in practice.By weighting the neural network's output of NIC, a modified novel information criterion (MNIC) algorithm is presented.MNIC extractes the principal components and corresponding eigenvectors in a parallel online learning program, and overcomes the NIC's defect.It is proved to have a single global optimum and nonquadratic convergence rate, which is superior to the conventional PCA online algorithms such as Oja and LMSER.The relationship among Oja, LMSER and MNIC is exhibited.Simulations show that MNIC could converge to the optimum fast.The validity of MNIC is proved.
Sox17 Regulates Liver Lipid Metabolism and Adaptation to Fasting
Rommelaere, Samuel; Millet, Virginie; Vu Manh, Thien-Phong; Gensollen, Thomas; Andreoletti, Pierre; Cherkaoui-Malki, Mustapha; Bourges, Christophe; Escalière, Bertrand; Du, Xin; Xia, Yu; Imbert, Jean; Beutler, Bruce; Kanai, Yoshiakira; Malissen, Bernard; Malissen, Marie
2014-01-01
Liver is a major regulator of lipid metabolism and adaptation to fasting, a process involving PPARalpha activation. We recently showed that the Vnn1 gene is a PPARalpha target gene in liver and that release of the Vanin-1 pantetheinase in serum is a biomarker of PPARalpha activation. Here we set up a screen to identify new regulators of adaptation to fasting using the serum Vanin-1 as a marker of PPARalpha activation. Mutagenized mice were screened for low serum Vanin-1 expression. Functional...
An Optimal Control Modification to Model-Reference Adaptive Control for Fast Adaptation
Nguyen, Nhan T.; Krishnakumar, Kalmanje; Boskovic, Jovan
2008-01-01
This paper presents a method that can achieve fast adaptation for a class of model-reference adaptive control. It is well-known that standard model-reference adaptive control exhibits high-gain control behaviors when a large adaptive gain is used to achieve fast adaptation in order to reduce tracking error rapidly. High gain control creates high-frequency oscillations that can excite unmodeled dynamics and can lead to instability. The fast adaptation approach is based on the minimization of the squares of the tracking error, which is formulated as an optimal control problem. The necessary condition of optimality is used to derive an adaptive law using the gradient method. This adaptive law is shown to result in uniform boundedness of the tracking error by means of the Lyapunov s direct method. Furthermore, this adaptive law allows a large adaptive gain to be used without causing undesired high-gain control effects. The method is shown to be more robust than standard model-reference adaptive control. Simulations demonstrate the effectiveness of the proposed method.
Multipole expansions in magnetostatics
Multipole expansions of the magnetic field of a spatially restricted system of stationary currents and those for the potential function of such currents in an external magnetic field are studied using angular momentum algebraic techniques. It is found that the expansion for the magnetic induction vector is made identical to that for the electric field strength of a neutral system of charges by substituting electric for magnetic multipole moments. The toroidal part of the multipole expansion for the magnetic field vector potential can, due to its potential nature, be omitted in the static case. Also, the potential function of a system of currents in an external magnetic field and the potential energy of a neutral system of charges in an external electric field have identical multipole expansions. For axisymmetric systems, the expressions for the field and those for the potential energy of electric and magnetic multipoles are reduced to simple forms, with symmetry axis orientation dependence separated out. (methodological notes)
Multipole expansions in magnetostatics
Agre, Mark Ya [National University of ' Kyiv-Mohyla Academy' , Kyiv (Ukraine)
2011-02-28
Multipole expansions of the magnetic field of a spatially restricted system of stationary currents and those for the potential function of such currents in an external magnetic field are studied using angular momentum algebraic techniques. It is found that the expansion for the magnetic induction vector is made identical to that for the electric field strength of a neutral system of charges by substituting electric for magnetic multipole moments. The toroidal part of the multipole expansion for the magnetic field vector potential can, due to its potential nature, be omitted in the static case. Also, the potential function of a system of currents in an external magnetic field and the potential energy of a neutral system of charges in an external electric field have identical multipole expansions. For axisymmetric systems, the expressions for the field and those for the potential energy of electric and magnetic multipoles are reduced to simple forms, with symmetry axis orientation dependence separated out. (methodological notes)
左风丽; 刘旭; 张宝印; 胡晓燕
2013-01-01
A parallel software module of FMM ( fast multipole method) for three-dimensional Laplace kernel functions, JASMIN-3DLapFMM, is designed and implemented. The module is based on two phases parallel strategy of both processes and threads. A parallel software module is successfully used to solve the far field potential of electrostatic fields. With a fixed problem size of single processor,almost linear weak parallel scalability is obtained for a grand scale problem with 10 particles on more than ten thousand processors. With a fixed size of total problem and 1 024 processors, about three times speedup is obtained on four threads.%在JASMIN上,基于进程/线程两级并行实现策略,研制三维Laplace核函数FMM (fast multipole method)的解法器模块“JASMIN-3DLapFMM”.该解法器已成功应用于三维静电场远场势的并行计算.固定单机问题规模,在上万个处理器核上运行百亿粒子的大规模问题,获得进程级几乎线性的并行可扩展性.固定总的问题规模和1 024个进程,4个线程时,获得大约3倍的加速.
Lim, T.
2011-04-28
To simulate numerically a non-destructive by eddy current testing (NDT-CF), the sensor response can be modeled through a semi-analytical approach by volume integral equations. Faster than the finite element method, this approach is however restricted to the study of plane or cylindrical parts (without taking into account the edge effects) because of the complexity of the expression of the dyadic Green function for more general configurations. However, there is an industrial demand to extend the capabilities of the CF model in complex configurations (deformed plates, edges effects...). We were thus brought to formulate the electromagnetic problem differently, by setting ourselves the goal of maintaining a semi-analytical approach. The surface integral equation (SIE) expresses the volume problem by an equivalent transmission one at the interfaces (2D) between homogeneous sub-domains. This problem is approached by a linear system (by the method of moments), whose number of unknowns is reduced due to the nature of the surfacic mesh. Therefore, this system can be solved by a direct solver for small configurations. That enabled us to treat several various positions of the sensor for only one inversion of the impedance matrix. The numerical results obtained using this formulation involve plates with consideration of edge effects such as edge and corner. They are consistent with results obtained by the finite element method. For larger configurations, we conducted a preliminary study for the adaptation of an acceleration method of the matrix vector product involved in an iterative solver (fast multipole method or FMM) to define the conditions under which the FMM calculation works correctly (accuracy, convergence...) in the NDT's domain. A special attention has been given to the choice of basis functions (which have to satisfy an Hdiv conforming property) and on the evaluation of near interactions (which are weakly singular). (author) [French] Pour simuler
Brain source localization based on fast fully adaptive approach.
Ravan, Maryam; Reilly, James P
2012-01-01
In the electroencephalogram (EEG) or magnetoencephalogram (MEG) context, brain source localization (beamforming) methods often fail when the number of observations is small. This is particularly true when measuring evoked potentials, especially when the number of electrodes is large. Due to the nonstationarity of the EEG/MEG, an adaptive capability is desirable. Previous work has addressed these issues by reducing the adaptive degrees of freedom (DoFs). This paper develops and tests a new multistage adaptive processing for brain source localization that has been previously used for radar statistical signal processing application with uniform linear antenna array. This processing, referred to as the fast fully adaptive (FFA) approach, could significantly reduce the required sample support and computational complexity, while still processing all available DoFs. The performance improvement offered by the FFA approach in comparison to the fully adaptive minimum variance beamforming (MVB) with limited data is demonstrated by bootstrapping simulated data to evaluate the variability of the source location. PMID:23367106
A Fast Adaptive Receive Antenna Selection Method in MIMO System
Chaowei Wang
2013-01-01
Full Text Available Antenna selection has been regarded as an effective method to acquire the diversity benefits of multiple antennas while potentially reduce hardware costs. This paper focuses on receive antenna selection. According to the proportion between the numbers of total receive antennas and selected antennas and the influence of each antenna on system capacity, we propose a fast adaptive antenna selection algorithm for wireless multiple-input multiple-output (MIMO systems. Mathematical analysis and numerical results show that our algorithm significantly reduces the computational complexity and memory requirement and achieves considerable system capacity gain compared with the optimal selection technique in the same time.
Tang, Guanglin; Yang, Ping; Sun, Bingqiang; Panetta, R. Lee; Kattawar, George W.
2016-06-01
The finite-difference time-domain (FDTD) and ray-by-ray (RBR) methods are techniques used to calculate the optical properties of nonspherical particles for small-to-moderate and large size parameters, respectively. The former is a rigorous method, and the latter is an approximate geometric-physical optics-hybrid method that takes advantage of both high efficiency of the geometric optics approach and high accuracy of the physical optics approach. In these two methods, the far field is calculated by mapping the near field to the far field with consideration of the phase interference. The mapping computation is more time-consuming than the near-field simulation when multiple scattering directions are involved, particularly in the case of the RBR implementation. To overcome the difficulty, in this study the fast multi-pole method is applied to both FDTD and RBR towards accelerating the far-field calculation, without degrading the accuracy of the simulation results.
Fast and Adaptive Lossless Onboard Hyperspectral Data Compression System
Aranki, Nazeeh I.; Keymeulen, Didier; Kimesh, Matthew A.
2012-01-01
Modern hyperspectral imaging systems are able to acquire far more data than can be downlinked from a spacecraft. Onboard data compression helps to alleviate this problem, but requires a system capable of power efficiency and high throughput. Software solutions have limited throughput performance and are power-hungry. Dedicated hardware solutions can provide both high throughput and power efficiency, while taking the load off of the main processor. Thus a hardware compression system was developed. The implementation uses a field-programmable gate array (FPGA). The implementation is based on the fast lossless (FL) compression algorithm reported in Fast Lossless Compression of Multispectral-Image Data (NPO-42517), NASA Tech Briefs, Vol. 30, No. 8 (August 2006), page 26, which achieves excellent compression performance and has low complexity. This algorithm performs predictive compression using an adaptive filtering method, and uses adaptive Golomb coding. The implementation also packetizes the coded data. The FL algorithm is well suited for implementation in hardware. In the FPGA implementation, one sample is compressed every clock cycle, which makes for a fast and practical realtime solution for space applications. Benefits of this implementation are: 1) The underlying algorithm achieves a combination of low complexity and compression effectiveness that exceeds that of techniques currently in use. 2) The algorithm requires no training data or other specific information about the nature of the spectral bands for a fixed instrument dynamic range. 3) Hardware acceleration provides a throughput improvement of 10 to 100 times vs. the software implementation. A prototype of the compressor is available in software, but it runs at a speed that does not meet spacecraft requirements. The hardware implementation targets the Xilinx Virtex IV FPGAs, and makes the use of this compressor practical for Earth satellites as well as beyond-Earth missions with hyperspectral instruments.
Holographic Adaptive Laser Optics System (HALOS): Fast, Autonomous Aberration Correction
Andersen, G.; MacDonald, K.; Gelsinger-Austin, P.
2013-09-01
We present an adaptive optics system which uses a multiplexed hologram to deconvolve the phase aberrations in an input beam. This wavefront characterization is extremely fast as it is based on simple measurements of the intensity of focal spots and does not require any computations. Furthermore, the system does not require a computer in the loop and is thus much cheaper, less complex and more robust as well. A fully functional, closed-loop prototype incorporating a 32-element MEMS mirror has been constructed. The unit has a footprint no larger than a laptop but runs at a bandwidth of 100kHz over an order of magnitude faster than comparable, conventional systems occupying a significantly larger volume. Additionally, since the sensing is based on parallel, all-optical processing, the speed is independent of actuator number running at the same bandwidth for one actuator as for a million. We are developing the HALOS technology with a view towards next-generation surveillance systems for extreme adaptive optics applications. These include imaging, lidar and free-space optical communications for unmanned aerial vehicles and SSA. The small volume is ideal for UAVs, while the high speed and high resolution will be of great benefit to the ground-based observation of space-based objects.
Fast, multiphase volume adaptation to hyperosmotic shock by Escherichia coli.
Teuta Pilizota
Full Text Available All living cells employ an array of different mechanisms to help them survive changes in extra cellular osmotic pressure. The difference in the concentration of chemicals in a bacterium's cytoplasm and the external environment generates an osmotic pressure that inflates the cell. It is thought that the bacterium Escherichia coli use a number of interconnected systems to adapt to changes in external pressure, allowing them to maintain turgor and live in surroundings that range more than two-hundred-fold in external osmolality. Here, we use fluorescence imaging to make the first measurements of cell volume changes over time during hyperosmotic shock and subsequent adaptation on a single cell level in vivo with a time resolution on the order of seconds. We directly observe two previously unseen phases of the cytoplasmic water efflux upon hyperosmotic shock. Furthermore, we monitor cell volume changes during the post-shock recovery and observe a two-phase response that depends on the shock magnitude. The initial phase of recovery is fast, on the order of 15-20 min and shows little cell-to-cell variation. For large sucrose shocks, a secondary phase that lasts several hours adds to the recovery. We find that cells are able to recover fully from shocks as high as 1 Osmol/kg using existing systems, but that for larger shocks, protein synthesis is required for full recovery.
Choi, Cheol Ho
2004-02-22
A new way of generating the multipole moments of Cartesian Gaussian functions in spherical polar coordinates has been established, bypassing the intermediary of Cartesian moment tensors. A new set of recurrence relations have also been derived for the resulting analytic integral values. The new method furnishes a conceptually simple and numerically efficient evaluation procedure for the multipole moments. The advantages over existing methods are documented. The results are relevant for the linear scaling quantum theories based on the fast multipole method. PMID:15268515
Multipoles are being employed as devices to study fusion issues and plasma phenomena at high values of beta (plasma pressure/magnetic pressure) in a controlled manner. Due to their large volume, low magnetic field (low synchrotron radiation) region, they are also under consideration as potential steady state advanced fuel (low neutron yield) reactors. Present experiments are investigating neoclassical (bootstrap and Pfirsch-Schlueter) currents and plasma stability at extremely high beta
Fast Adaptive Blind MMSE Equalizer for Multichannel FIR Systems
Kacha, Ibrahim; Abed-Meraim, Karim; Belouchrani, Adel
2006-12-01
We propose a new blind minimum mean square error (MMSE) equalization algorithm of noisy multichannel finite impulse response (FIR) systems, that relies only on second-order statistics. The proposed algorithm offers two important advantages: a low computational complexity and a relative robustness against channel order overestimation errors. Exploiting the fact that the columns of the equalizer matrix filter belong both to the signal subspace and to the kernel of truncated data covariance matrix, the proposed algorithm achieves blindly a direct estimation of the zero-delay MMSE equalizer parameters. We develop a two-step procedure to further improve the performance gain and control the equalization delay. An efficient fast adaptive implementation of our equalizer, based on the projection approximation and the shift invariance property of temporal data covariance matrix, is proposed for reducing the computational complexity from[InlineEquation not available: see fulltext.] to[InlineEquation not available: see fulltext.], where[InlineEquation not available: see fulltext.] is the number of emitted signals,[InlineEquation not available: see fulltext.] the data vector length, and[InlineEquation not available: see fulltext.] the dimension of the signal subspace. We then derive a statistical performance analysis to compare the equalization performance with that of the optimal MMSE equalizer. Finally, simulation results are provided to illustrate the effectiveness of the proposed blind equalization algorithm.
Fast Adaptive Blind MMSE Equalizer for Multichannel FIR Systems
Abed-Meraim Karim
2006-01-01
Full Text Available We propose a new blind minimum mean square error (MMSE equalization algorithm of noisy multichannel finite impulse response (FIR systems, that relies only on second-order statistics. The proposed algorithm offers two important advantages: a low computational complexity and a relative robustness against channel order overestimation errors. Exploiting the fact that the columns of the equalizer matrix filter belong both to the signal subspace and to the kernel of truncated data covariance matrix, the proposed algorithm achieves blindly a direct estimation of the zero-delay MMSE equalizer parameters. We develop a two-step procedure to further improve the performance gain and control the equalization delay. An efficient fast adaptive implementation of our equalizer, based on the projection approximation and the shift invariance property of temporal data covariance matrix, is proposed for reducing the computational complexity from to , where is the number of emitted signals, the data vector length, and the dimension of the signal subspace. We then derive a statistical performance analysis to compare the equalization performance with that of the optimal MMSE equalizer. Finally, simulation results are provided to illustrate the effectiveness of the proposed blind equalization algorithm.
Multipole moments for embedding potentials
Nørby, Morten Steen; Olsen, Jógvan Magnus Haugaard; Kongsted, Jacob;
2016-01-01
and RNA/DNA, and of molecules in solution. Our aim is to develop a computational methodology for distributed multipole moments and their associated multipole polarizabilities which is accurate, computationally efficient, and with smooth convergence with respect to multipole order. As the first step...... toward this goal, we herein investigate different ways of obtaining distributed atom-centered multipole moments that are used in the construction of the electrostatic part of the embedding potential. Our objective is methods that not only are accurate and computationally efficient, but which can be...
El-Shenawee, M; Rappaport, C; Silevitch, M
2001-12-01
We present a statistical study of the electric field scattered from a three-dimensional penetrable object buried under a two-dimensional random rough surface. Monte Carlo simulations using the steepest-descent fast multipole method (SDFMM) are conducted to calculate the average and the standard deviation of the near-zone scattered fields. The SDFMM, originally developed at the University of Illinois at Urbana-Champaign, has been modified to calculate the unknown surface currents both on the rough ground and on the buried object that are due to excitation by a tapered Gaussian beam. The rough ground medium used is an experimentally measured typical dry Bosnian soil with 3.8% moisture, while the buried object represents a plastic land mine modeled as an oblate spheroid with dimensions and burial depth smaller than the free-space wavelength. Both vertical and horizontal polarizations for the incident waves are studied. The numerical results show that the TNT mine signature is almost 5% of the total field scattered from the ground. Moreover, relatively recognizable object signatures are observed even when the object is buried under the tail of the incident beam. Interestingly, even for the small surface roughness parameters considered, the standard deviation of the object signature is almost 30% of the signal itself, indicating significant clutter distortion that is due to the roughness of the ground. PMID:11760205
The origin of soft multipole states in halo nuclei is discussed by using a double square well potential. It is found that the huge enhancement of the transition strength near particle threshold energy is caused by the coherent contribution of particle and hole wave functions from an extremely large radial distance R=10∼45 fm. The soft excitation is thus characterized as an independent particle-hole excitation rather than a coherent superposition of particle-hole states. The characteristic feature of the soft mode is demonstrated to manifest itself in the response to the momentum dependent transition operator. (author) 12 refs., 4 figs., 1 tab
Multipole Expansion in Generalized Electrodynamics
Bonin, C A; Ortega, P H
2016-01-01
In this article we study some classical aspects of Podolsky Electrodynamics in the static regime. We develop the multipole expansion for the theory in both the electrostatic and the magnetostatic cases. We also address the problem of consistently truncating the infinite series associated with the several kinds of multipoles, yielding approximations for the static Podolskian electromagnetic field to any degree of precision required. Moreover, we apply the general theory of multipole expansion to some specific physical problems. In those problems we identify the first terms of the series with the monopole, dipole and quadrupole terms in the generalized theory. We also propose a situation in which Podolsky theory can be experimentally tested.
FFA-CONTEXT AWARE ENERGY EFFICIENT ROUTING USING FAST REACTIVE AND ADAPTIVE ALGORITHM
P. Kavitha Rani; E. Kannan
2014-01-01
A vital issue of routing is how to efficiently strengthen energy consumption of the whole network and to avoid the energy-hole which leads to node failure or node death; this scenario makes the network to work improperly and leads to network death, Here we proves our hypothesis with linearity between routing and linear principles. Our proposed work effectively utilizes the adaption of Optimized particle swarm algorithm (Fast Reactive and fast adaptive algorithm) in order to solve the routing ...
Multipole moments of isolated horizons
To every axi-symmetric isolated horizon we associate two sets of numbers, Mn and Jn with n = 0, 1, 2, ..., representing its mass and angular momentum multipoles. They provide a diffeomorphism invariant characterization of the horizon geometry. Physically, they can be thought of as the 'source multipoles' of black holes in equilibrium. These structures have a variety of potential applications ranging from equations of motion of black holes and numerical relativity to quantum gravity
Multipole structure of compact objects
Quevedo, Hernando
2016-01-01
We analyze the applications of general relativity in relativistic astrophysics in order to solve the problem of describing the geometric and physical properties of the interior and exterior gravitational and electromagnetic fields of compact objects. We focus on the interpretation of exact solutions of Einstein's equations in terms of their multipole moments structure. In view of the lack of physical interior solutions, we propose an alternative approach in which higher multipoles should be taken into account.
Measured multipole moments of continuum electron transfer angular distributions
The velocity space distribution of electrons emitted near the forward direction from collisions involving fast, highly stripped oxygen ions with gaseous and solid targets is presented and described in terms of multipole moments of the ejected charge distribution, which permits direct comparison with recent theory. The measurements are produced by employing position-sensitive electron detection to combine emission angle definition with conventional electrostatic spectrometry. Agreement obtained between theory and distributions observed for binary continuum electron loss processes coupled with a similar multipole content observed with solid targets suggests a model of convoy electron production dominated by electron loss from the projectile within the bulk of the target. Further, the connection between multipoles of the projectile electron emission distribution in single collisions and the state of excitation of that projectile excited states may provide the basis for a probe of the state of ions traversing bulk solid matter. 14 refs., 3 figs., 1 tab
Fast calibration of high-order adaptive optics systems
Kasper, Markus; Fedrigo, Enrico; Looze, Douglas P.; Bonnet, Henri; Ivanescu, Liviu; Oberti, Sylvain
2004-06-01
We present a new method of calibrating adaptive optics systems that greatly reduces the required calibration time or, equivalently, improves the signal-to-noise ratio. The method uses an optimized actuation scheme with Hadamard patterns and does not scale with the number of actuators for a given noise level in the wave-front sensor channels. It is therefore highly desirable for high-order systems and/or adaptive secondary systems on a telescope without a Gregorian focal plane. In the latter case, the measurement noise is increased by the effects of the turbulent atmosphere when one is calibrating on a natural guide star.
Centralized and Distributed Solutions for Fast Muting Adaptation in LTE-Advanced HetNets
Soret, Beatriz; Pedersen, Klaus I.
2015-01-01
though the eICIC configuration (RE and ABS) ideally should be instantaneously adapted to follow the fluctuations of the traffic and the channel conditions over time, previous studies have focused on slow intercell coordination. In this paper, we investigate fast dynamic eICIC solutions for centralized...... and distributed Radio Resource Management (RRM) architectures. The centralized RRM architecture assumes macro and Remote Radio Heads (RRHs) inter-connected via high-speed fronthaul connections, while the distributed architecture is based on traditional macro and pico cells deployments with X2 backhaul...... interface. Two different fast muting adaptation algorithms are derived, and it is shown how those can be appplied to both the centralized and the distributed architecture. Performance results with bursty traffic show that the fast dynamic adaptation provides significant gains, both in 5%-ile and 50%-ile...
Adaptation of fibers in fast-twitch muscles of rats to spaceflight and hindlimb suspension
Jiang, Bian; Ohira, Yoshi; Roy, Roland R.; Nguyen, Quyet; Il'ina-Kakueva, E. I.; Oganov, V.; Edgerton, V. R.
1992-01-01
The adaptation of single fibers in medial gastrocnemius (MG), a fast-twitch extensor, and in tibialis anterior (TA), a fast-twitch flexor, was studied after 14 days of spaceflight onboard Cosmos 2044 or hindlimb suspension. Quantitative myosin ATPase activities of single fibers were measured in flight and suspended rats. Each of the enzyme and size measurements were directly correlated within each fiber with respect to its qualitative myosin ATPase staining properties and its expression of fast, slow, or both myosin heavy chains (MHC). The percentage of slow- and fast-twitch fibers of the MG and TA were found to be unchanged. Mean fiber size of all fibers was unaffected after flight or suspension. The ATPase activity in the MG was higher in flight than in control or suspended rats. In comparison to Cosmos 1887 spaceflight, the adaptations in the muscle fibers of the MG were more moderate.
Radiation reaction for multipole moments
Kazinski, P O
2006-01-01
We propose a Poincare-invariant description for the effective dynamics of systems of charged particles by means of intrinsic multipole moments. To achieve this goal we study the effective dynamics of such systems within two frameworks -- the particle itself and hydrodynamical one. We give a relativistic-invariant definition for the intrinsic multipole moments both pointlike and extended relativistic objects. Within the hydrodynamical framework we suggest a covariant action functional for a perfect fluid with pressure. In the case of a relativistic charged dust we prove the equivalence of the particle approach to the hydrodynamical one to the problem of radiation reaction for multipoles. As the particular example of a general procedure we obtain the effective model for a neutral system of charged particles with dipole moment.
Radiation reaction of multipole moments
A Poincare-invariant description is proposed for the effective dynamics of a localized system of charged particles in classical electrodynamics in terms of the intrinsic multipole moments of the system. A relativistic-invariant definition for the intrinsic multipole moments of a system of charged particles is given. A new generally covariant action functional for a relativistic perfect fluid is proposed. In the case of relativistic charged dust, it is proven that the description of the problem of radiation reaction of multipole moments by the model of particles is equivalent to the description of this problem by a hydrodynamic model. An effective model is obtained for a pointlike neutral system of charged particles that possesses an intrinsic dipole moment, and the free dynamics of this system is described. The bound momentum of a point dipole is found
Radiation reaction of multipole moments
Kazinski, P. O.
2007-08-01
A Poincaré-invariant description is proposed for the effective dynamics of a localized system of charged particles in classical electrodynamics in terms of the intrinsic multipole moments of the system. A relativistic-invariant definition for the intrinsic multipole moments of a system of charged particles is given. A new generally covariant action functional for a relativistic perfect fluid is proposed. In the case of relativistic charged dust, it is proven that the description of the problem of radiation reaction of multipole moments by the model of particles is equivalent to the description of this problem by a hydrodynamic model. An effective model is obtained for a pointlike neutral system of charged particles that possesses an intrinsic dipole moment, and the free dynamics of this system is described. The bound momentum of a point dipole is found.
Improving GPU-accelerated Adaptive IDW Interpolation Algorithm Using Fast kNN Search
Mei, Gang; Xu, Nengxiong; Xu, Liangliang
2016-01-01
This paper presents an efficient parallel Adaptive Inverse Distance Weighting (AIDW) interpolation algorithm on modern Graphics Processing Unit (GPU). The presented algorithm is an improvement of our previous GPU-accelerated AIDW algorithm by adopting fast k-Nearest Neighbors (kNN) search. In AIDW, it needs to find several nearest neighboring data points for each interpolated point to adaptively determine the power parameter; and then the desired prediction value of the interpolated point is ...
Identification of fast-changing signals by means of adaptive chaotic transformations
Berezowski, Marek; Lawnik, Marcin
2016-01-01
The adaptive approach of strongly non-linear fast-changing signals identification is discussed. The approach is devised by adaptive sampling based on chaotic mapping in yourself of a signal. Presented sampling way may be utilized online in the automatic control of chemical reactor (throughout identification of concentrations and temperature oscillations in real-time), in medicine (throughout identification of ECG and EEG signals in real-time), etc. In this paper, we presented it to identify t...
Adaptive approach to global synchronization of directed networks with fast switching topologies
Global synchronization of directed networks with switching topologies is investigated. It is found that if there exists at least one directed spanning tree in the network with the fixed time-average topology and the time-average topology is achieved sufficiently fast, the network will reach global synchronization for appreciate coupling strength. Furthermore, this appreciate coupling strength may be obtained by local adaptive approach. A sufficient condition about the global synchronization is given. Numerical simulations verify the effectiveness of the adaptive strategy.
Faust, Thomas W.; Assous, Maxime; Shah, Fulva; Tepper, James M.; Koós, Tibor
2015-01-01
Previous work suggests that neostriatal cholinergic interneurons control the activity of several classes of GABAergic interneurons through fast nicotinic receptor mediated synaptic inputs. Although indirect evidence has suggested the existence of several classes of interneurons controlled by this mechanism only one such cell type, the neuropeptide-Y expressing neurogliaform neuron, has been identified to date. Here we tested the hypothesis that in addition to the neurogliaform neurons that el...
Polarizable atomic multipole X-ray refinement: application to peptide crystals
Schnieders, Michael J.; Fenn, Timothy D.; Pande, Vijay S.; Brunger, Axel T.
2009-01-01
Recent advances in computational chemistry have produced force fields based on a polarizable atomic multipole description of biomolecular electrostatics. In this work, the Atomic Multipole Optimized Energetics for Biomolecular Applications (AMOEBA) force field is applied to restrained refinement of molecular models against X-ray diffraction data from peptide crystals. A new formalism is also developed to compute anisotropic and aspherical structure factors using fast Fourier transformation (F...
Static multipole deformations in nuclei
The physics of static multipole deformations in nuclei is reviewed. Nuclear static moments result from the delicate balance between the vibronic Jahn-Teller interaction (particle-vibration coupling) and the residual interaction (pairing force). Examples of various permanent nuclear deformations are discussed
Fast Block-match Motion Estimation Based on Multilevel Adaptive Diamond Search
Li, Shan; Yi, Qing-Ming; Shi, Min
In this paper, a novel fast block-match algorithm called MADS based on multilevel adaptive diamond search is proposed. The algorithm adaptively estimates the frame-level motion complexity with the reference frame texture information and the macro-block residual value at first, and then estimates the block-level motion complexity according to the spatial-temporal correlation of the vector field. The threshold is applied to stop the stationary block from searching. The initial search point and different diamond search modes are adaptively selected based on motion type for non-stationary block. Experimental results show that MADS algorithm has better performance than other popular fast algorithms for a wide range of video sequences.
Mathematical Programming Solution for the Frictional Contact Multipole BEM
YU Chunxiao; SHEN Guangxian; LIU Deyi
2005-01-01
This paper presents a new mathematical model for the highly nonlinear problem of frictional contact. A programming model, multipole boundary element method (BEM), was developed for 3-D elastic contact with friction to replace the Monte Carlo method. A numerical example shows that the optimization programming model for the point-to-surface contact with friction and the fast optimization generalized minimal residual algorithm (GMRES(m)) significantly improve the analysis of such problems relative to the conventional BEM.
Fast non-local means with size-adaptive search window
Hancheng Yu
2016-03-01
Full Text Available In this study, the authors present a fast non-local means (NLM image denoising algorithm with size-adaptive search window. On the basis of the edge gradient and direction of the noisy image, the proposed fast scheme divides all the pixels into significant edge, moderate edge or non-edge region, and the size of the search window for most pixels which belong to non-edge region can be reduced. The proposed fast scheme also adopts different strategies to pre-select similar patches in the search window for efficient NLM denoising. Experimental results show that compared with the standard NLM method, the proposed fast scheme achieves a substantial reduction in computational cost and improvement in the denoising performance, both in terms of visual quality and numerical results.
Polarizable atomic multipole X-ray refinement: application to peptide crystals
A method to accelerate the computation of structure factors from an electron density described by anisotropic and aspherical atomic form factors via fast Fourier transformation is described for the first time. Recent advances in computational chemistry have produced force fields based on a polarizable atomic multipole description of biomolecular electrostatics. In this work, the Atomic Multipole Optimized Energetics for Biomolecular Applications (AMOEBA) force field is applied to restrained refinement of molecular models against X-ray diffraction data from peptide crystals. A new formalism is also developed to compute anisotropic and aspherical structure factors using fast Fourier transformation (FFT) of Cartesian Gaussian multipoles. Relative to direct summation, the FFT approach can give a speedup of more than an order of magnitude for aspherical refinement of ultrahigh-resolution data sets. Use of a sublattice formalism makes the method highly parallelizable. Application of the Cartesian Gaussian multipole scattering model to a series of four peptide crystals using multipole coefficients from the AMOEBA force field demonstrates that AMOEBA systematically underestimates electron density at bond centers. For the trigonal and tetrahedral bonding geometries common in organic chemistry, an atomic multipole expansion through hexadecapole order is required to explain bond electron density. Alternatively, the addition of interatomic scattering (IAS) sites to the AMOEBA-based density captured bonding effects with fewer parameters. For a series of four peptide crystals, the AMOEBA–IAS model lowered Rfree by 20–40% relative to the original spherically symmetric scattering model
Polarizable atomic multipole X-ray refinement: application to peptide crystals
Schnieders, Michael J. [Department of Chemistry, Stanford, CA 94305 (United States); Fenn, Timothy D. [Department of Molecular and Cellular Physiology, Stanford, CA 94305 (United States); Howard Hughes Medical Institute (United States); Pande, Vijay S., E-mail: pande@stanford.edu [Department of Chemistry, Stanford, CA 94305 (United States); Brunger, Axel T., E-mail: pande@stanford.edu [Department of Molecular and Cellular Physiology, Stanford, CA 94305 (United States); Howard Hughes Medical Institute (United States); Department of Chemistry, Stanford, CA 94305 (United States)
2009-09-01
A method to accelerate the computation of structure factors from an electron density described by anisotropic and aspherical atomic form factors via fast Fourier transformation is described for the first time. Recent advances in computational chemistry have produced force fields based on a polarizable atomic multipole description of biomolecular electrostatics. In this work, the Atomic Multipole Optimized Energetics for Biomolecular Applications (AMOEBA) force field is applied to restrained refinement of molecular models against X-ray diffraction data from peptide crystals. A new formalism is also developed to compute anisotropic and aspherical structure factors using fast Fourier transformation (FFT) of Cartesian Gaussian multipoles. Relative to direct summation, the FFT approach can give a speedup of more than an order of magnitude for aspherical refinement of ultrahigh-resolution data sets. Use of a sublattice formalism makes the method highly parallelizable. Application of the Cartesian Gaussian multipole scattering model to a series of four peptide crystals using multipole coefficients from the AMOEBA force field demonstrates that AMOEBA systematically underestimates electron density at bond centers. For the trigonal and tetrahedral bonding geometries common in organic chemistry, an atomic multipole expansion through hexadecapole order is required to explain bond electron density. Alternatively, the addition of interatomic scattering (IAS) sites to the AMOEBA-based density captured bonding effects with fewer parameters. For a series of four peptide crystals, the AMOEBA–IAS model lowered R{sub free} by 20–40% relative to the original spherically symmetric scattering model.
Multipole correction in large synchrotrons
A new method of correcting dynamic nonlinearities due to the multipole content of a synchrotron such as the Superconducting Super Collider is discussed. The method uses lumped multipole elements placed at the center (C) of the accelerator half-cells as well as elements near the focusing (F) and defocusing (D) quads. In a first approximation, the corrector strengths follow Simpson's Rule. Correction of second-order sextupole nonlinearities may also be obtained with the F, C, and D octupoles. Correction of nonlinearities by about three orders of magnitude are obtained, and simple solutions to a fundamental problem in synchrotrons are demonstrated. Applications to the CERN Large Hadron Collider and lower energy machines, as well as extensions for quadrupole correction, are also discussed
Harmonic projection and multipole Vectors
Lachièze-Rey, M
2004-01-01
We show that the multipole vector decomposition, recently introduced by Copi et al., is a consequence of Sylvester s theorem, and corresponds to the Maxwell representation. Analyzing it in terms of harmonic polynomials, we show that this decomposition, as well as its extension proposed by Katz and Weeks, result in fact from the application of the harmonic projection operator and its inverse. We derive the coefficients of the usual harmonic decomposition from the multipole vectors. We answer to '' an open question '' asked by Katz and Weeks, by showing that the decomposition resulting from their corollary is unstable. We propose however a new decomposition which is stable. We generalize these results to complex functions and polynomials.
Fast implementation of length-adaptive privacy amplification in quantum key distribution
Post-processing is indispensable in quantum key distribution (QKD), which is aimed at sharing secret keys between two distant parties. It mainly consists of key reconciliation and privacy amplification, which is used for sharing the same keys and for distilling unconditional secret keys. In this paper, we focus on speeding up the privacy amplification process by choosing a simple multiplicative universal class of hash functions. By constructing an optimal multiplication algorithm based on four basic multiplication algorithms, we give a fast software implementation of length-adaptive privacy amplification. “Length-adaptive” indicates that the implementation of privacy amplification automatically adapts to different lengths of input blocks. When the lengths of the input blocks are 1 Mbit and 10 Mbit, the speed of privacy amplification can be as fast as 14.86 Mbps and 10.88 Mbps, respectively. Thus, it is practical for GHz or even higher repetition frequency QKD systems. (general)
Luithardt, Annette Franziska geb. Schmid
2016-01-01
Purpose To validate the commercially available “EyeSuite” version of the new fast thresholding algorithm GATE (German Adaptive Thresholding Estimation) for automated static perimetry. Methods Thirty patients suffering from visual pathway lesions of various origin (anterior ischemic optic neuropathy [n=3], glaucoma [n=15], (post-) chiasmal visual pathway lesion [n=6], retinitis pigmentosa (RP) [n=6]) were tested on three Octopus 900 perimeters (Haag-Streit AG, Köniz, Switzerland) wi...
Fast Linear Adaptive Skipping Training Algorithm for Training Artificial Neural Network
Manjula Devi, R.; R. C. Suganthe; S. KUPPUSWAMI
2013-01-01
Artificial neural network has been extensively consumed training model for solving pattern recognition tasks. However, training a very huge training data set using complex neural network necessitates excessively high training time. In this correspondence, a new fast Linear Adaptive Skipping Training (LAST) algorithm for training artificial neural network (ANN) is instituted. The core essence of this paper is to ameliorate the training speed of ANN by exhibiting only the input samples that do ...
Strength of a multipole residual interaction
The consequences of assuming different radial dependences for the multipole residual interaction are explored, using macroscopic and microscopic arguments. Calculations are done both for schematic and realistic cases
Improvement of the Simplified Fast Transversal Filter Type Algorithm for Adaptive Filtering
Madjid Arezki
2009-01-01
Full Text Available Problem statement: In this study, we proposed a new algorithm M-SMFTF for adaptive filtering with fast convergence and low complexity. Approach: It was the result of a simplified FTF type algorithm, where the adaptation gain was obtained only from the forward prediction variables and using a new recursive method to compute the likelihood variable. Results: The computational complexity was reduced from 7L-6L, where L is the finite impulse response filter length. Furthermore, this computational complexity can be significantly reduced to (2L+4P when used with a reduced P-size forward predictor. Conclusion: This algorithm presented a certain interest, for the adaptation of very long filters, like those used in the problems of echo acoustic cancellation, due to its reduced complexity, its numerical stability and its convergence in the presence of the speech signal.
Dorian S Houser
2013-11-01
Full Text Available Insulin resistance in modern society is perceived as a pathological consequence of excess energy consumption and reduced physical activity. Its presence in relation to the development of cardiovascular risk factors has been termed the metabolic syndrome, which produces increased mortality and morbidity and which is rapidly increasing in human populations. Ironically, insulin resistance likely evolved to assist animals during food shortages by increasing the availability of endogenous lipid for catabolism while protecting protein from use in gluconeogenesis and eventual oxidation. Some species that incorporate fasting as a predictable component of their life history demonstrate physiological traits similar to the metabolic syndrome during prolonged fasts. One such species is the northern elephant seal (Mirounga angustirostris, which fasts from food and water for periods of up to three months. During this time, ~90% of the seals metabolic demands are met through fat oxidation and circulating non-esterified fatty acids are high (0.7-3.2 mM. All life history stages of elephant seal studied to date demonstrate insulin resistance and fasting hyperglycemia as well as variations in hormones and adipocytokines that reflect the metabolic syndrome to some degree. Elephant seals demonstrate some intriguing adaptations with the potential for medical advancement; for example, ketosis is negligible despite significant and prolonged fatty acid oxidation and investigation of this feature might provide insight into the treatment of diabetic ketoacidosis. The parallels to the metabolic syndrome are likely reflected to varying degrees in other marine mammals, most of which evolved on diets high in lipid and protein content but essentially devoid of carbohydrate. Utilization of these natural models of insulin resistance may further our understanding of the pathophysiology of the metabolic syndrome in humans and better assist the development of preventative measures
Polarizable atomic multipole X-ray refinement: application to peptide crystals.
Schnieders, Michael J; Fenn, Timothy D; Pande, Vijay S; Brunger, Axel T
2009-09-01
Recent advances in computational chemistry have produced force fields based on a polarizable atomic multipole description of biomolecular electrostatics. In this work, the Atomic Multipole Optimized Energetics for Biomolecular Applications (AMOEBA) force field is applied to restrained refinement of molecular models against X-ray diffraction data from peptide crystals. A new formalism is also developed to compute anisotropic and aspherical structure factors using fast Fourier transformation (FFT) of Cartesian Gaussian multipoles. Relative to direct summation, the FFT approach can give a speedup of more than an order of magnitude for aspherical refinement of ultrahigh-resolution data sets. Use of a sublattice formalism makes the method highly parallelizable. Application of the Cartesian Gaussian multipole scattering model to a series of four peptide crystals using multipole coefficients from the AMOEBA force field demonstrates that AMOEBA systematically underestimates electron density at bond centers. For the trigonal and tetrahedral bonding geometries common in organic chemistry, an atomic multipole expansion through hexadecapole order is required to explain bond electron density. Alternatively, the addition of interatomic scattering (IAS) sites to the AMOEBA-based density captured bonding effects with fewer parameters. For a series of four peptide crystals, the AMOEBA-IAS model lowered R(free) by 20-40% relative to the original spherically symmetric scattering model. PMID:19690373
Adaptive Hybrid Visual Servo Regulation of Mobile Robots Based on Fast Homography Decomposition
Chunfu Wu
2015-01-01
Full Text Available For the monocular camera-based mobile robot system, an adaptive hybrid visual servo regulation algorithm which is based on a fast homography decomposition method is proposed to drive the mobile robot to its desired position and orientation, even when object’s imaging depth and camera’s position extrinsic parameters are unknown. Firstly, the homography’s particular properties caused by mobile robot’s 2-DOF motion are taken into account to induce a fast homography decomposition method. Secondly, the homography matrix and the extracted orientation error, incorporated with the desired view’s single feature point, are utilized to form an error vector and its open-loop error function. Finally, Lyapunov-based techniques are exploited to construct an adaptive regulation control law, followed by the experimental verification. The experimental results show that the proposed fast homography decomposition method is not only simple and efficient, but also highly precise. Meanwhile, the designed control law can well enable mobile robot position and orientation regulation despite the lack of depth information and camera’s position extrinsic parameters.
FFA-CONTEXT AWARE ENERGY EFFICIENT ROUTING USING FAST REACTIVE AND ADAPTIVE ALGORITHM
P. Kavitha Rani
2014-01-01
Full Text Available A vital issue of routing is how to efficiently strengthen energy consumption of the whole network and to avoid the energy-hole which leads to node failure or node death; this scenario makes the network to work improperly and leads to network death, Here we proves our hypothesis with linearity between routing and linear principles. Our proposed work effectively utilizes the adaption of Optimized particle swarm algorithm (Fast Reactive and fast adaptive algorithm in order to solve the routing problem in wireless sensor networks to avoid the energy hole. Our approach effectively defines the node leader based on energy level and path optimization for node traversing using particle swarm which effectively utilizes the solution for problem in linearity into routing problem. Our algorithm FFA redefines the particle swarm rules which are reliable in routing strategies and adapts to the working principle. The optimization in our algorithm is to balance energy level and to withstand for its context. This study proves the accuracy using our algorithm respectively by theoretical survey and analysis with simulated results.
The use of the spectral method within the fast adaptive composite grid method
McKay, S.M.
1994-12-31
The use of efficient algorithms for the solution of partial differential equations has been sought for many years. The fast adaptive composite grid (FAC) method combines an efficient algorithm with high accuracy to obtain low cost solutions to partial differential equations. The FAC method achieves fast solution by combining solutions on different grids with varying discretizations and using multigrid like techniques to find fast solution. Recently, the continuous FAC (CFAC) method has been developed which utilizes an analytic solution within a subdomain to iterate to a solution of the problem. This has been shown to achieve excellent results when the analytic solution can be found. The CFAC method will be extended to allow solvers which construct a function for the solution, e.g., spectral and finite element methods. In this discussion, the spectral methods will be used to provide a fast, accurate solution to the partial differential equation. As spectral methods are more accurate than finite difference methods, the ensuing accuracy from this hybrid method outside of the subdomain will be investigated.
A single-site multipole model for liquid water
Tran, Kelly N.; Tan, Ming-Liang; Ichiye, Toshiko
2016-07-01
Accurate and efficient empirical potential energy models that describe the atomistic interactions between water molecules in the liquid phase are essential for computer simulations of many problems in physics, chemistry, and biology, especially when long length or time scales are important. However, while models with non-polarizable partial charges at four or five sites in a water molecule give remarkably good values for certain properties, deficiencies have been noted in other properties and increasing the number of sites decreases computational efficiency. An alternate approach is to utilize a multipole expansion of the electrostatic potential due to the molecular charge distribution, which is exact outside the charge distribution in the limits of infinite distances or infinite orders of multipoles while partial charges are a qualitative representation of electron density as point charges. Here, a single-site multipole model of water is presented, which is as fast computationally as three-site models but is also more accurate than four- and five-site models. The dipole, quadrupole, and octupole moments are from quantum mechanical-molecular mechanical calculations so that they account for the average polarization in the liquid phase, and represent both the in-plane and out-of-plane electrostatic potentials of a water molecule in the liquid phase. This model gives accurate thermodynamic, dynamic, and dielectric properties at 298 K and 1 atm, as well as good temperature and pressure dependence of these properties.
Adaptive control of piezoelectric fast steering mirror for high precision tracking application
Wang, Geng; Rao, Changhui
2015-03-01
A piezoelectric fast steering mirror (PFSM) is a complex, strong coupling nonlinear system that integrates optics, mechanics, electrics, and control. Due to the existence of hysteresis nonlinearity, mechanical resonance, and all kinds of disturbances, precise tracking control of a PFSM is a challenging task. This paper presents a comprehensive study of modeling, controller design, and simulation evaluation for a PFSM system. First a general model of a PFSM system integrating mechanical dynamics, electrical dynamics, and hysteresis nonlinearity is proposed, and then a robust adaptive controller is developed under both unknown hysteresis nonlinearities and parameter uncertainties. The parameters needed directly in the formulation of the controller are adaptively estimated. The proposed control law ensures the uniform boundedness of all signals in the closed-loop system. Furthermore, a stability analysis of the control system is performed to guarantee that the output tracking error converges to zero asymptotically. Finally, simulation tests with different motion trajectories are conducted to verify the effectiveness of the proposed method.
Antonella Agodi
2005-03-01
Full Text Available
Background. Food contamination by trichothecene mycotoxins is considered to be an emerging public health problem. The aim of this study was to validate a rapid sonification protocol, previously set up for cereal Fusarium DNA extraction from fast food samples, produced by a centre for research and development in the food industry in Catania, Sicily, and to validate it for a diagnostic PCR assay targeted at tri5, the key gene of trichothecene biosynthesis.
Methods. DNA from reference Fusarium spp. strains and from fast food samples was prepared, setting up an extraction protocol adapted using some modifications based on a method recently described. Validation experiments were performed: serial dilution of DNA extracted from fungal cultures were added to food samples and then DNA was extracted. Specific primer pairs were used to detect F. graminearum and F. culmorum DNA in species-specific assays as well as trichothecene-producing Fusarium spp. in a groupspecific system.
Results. All genomic DNA extracted from trichothecene-producing Fusarium spp. as well as from DNA-spiked fast food samples and from food still in it’s original condition resulted in the correct amplification. The detection limit was 1 x 10-4 μg of DNA. All fungal and food samples tested gave highly consistent results in repeatability assays, thus demonstrating the within-lab and within/between-day precision of the method.
Conclusions. Information on the epidemiology of trichothecene producing Fusarium through the food chain and the identification of the most frequently contaminated components of fast food are essential in order to develop effective public health strategies for minimising consumer exposure to trichothecenes. Key words: Fusarium, fast food, trichothecenes, PCR
Windowed multipole for cross section Doppler broadening
Josey, C.; Ducru, P.; Forget, B.; Smith, K.
2016-02-01
This paper presents an in-depth analysis on the accuracy and performance of the windowed multipole Doppler broadening method. The basic theory behind cross section data is described, along with the basic multipole formalism followed by the approximations leading to windowed multipole method and the algorithm used to efficiently evaluate Doppler broadened cross sections. The method is tested by simulating the BEAVRS benchmark with a windowed multipole library composed of 70 nuclides. Accuracy of the method is demonstrated on a single assembly case where total neutron production rates and 238U capture rates compare within 0.1% to ACE format files at the same temperature. With regards to performance, clock cycle counts and cache misses were measured for single temperature ACE table lookup and for windowed multipole. The windowed multipole method was found to require 39.6% more clock cycles to evaluate, translating to a 7.9% performance loss overall. However, the algorithm has significantly better last-level cache performance, with 3 fewer misses per evaluation, or a 65% reduction in last-level misses. This is due to the small memory footprint of the windowed multipole method and better memory access pattern of the algorithm.
Adaptive AFM scan speed control for high aspect ratio fast structure tracking
Ahmad, Ahmad; Schuh, Andreas; Rangelow, Ivo W. [Department of Microelectronic and Nanoelectronic Systems, Faculty of Electrical Engineering and Information Technology Ilmenau University of Technology, Gustav-Kirchhoffstr. 1, 98684 Ilmenau (Germany)
2014-10-15
Improved imaging rates in Atomic Force Microscopes (AFM) are of high interest for disciplines such as life sciences and failure analysis of semiconductor wafers, where the sample topology shows high aspect ratios. Also, fast imaging is necessary to cover a large surface under investigation in reasonable times. Since AFMs are composed of mechanical components, they are associated with comparably low resonance frequencies that undermine the effort to increase the acquisition rates. In particular, high and steep structures are difficult to follow, which causes the cantilever to temporarily loose contact to or crash into the sample. Here, we report on a novel approach that does not affect the scanner dynamics, but adapts the lateral scanning speed of the scanner. The controller monitors the control error signal and, only when necessary, decreases the scan speed to allow the z-piezo more time to react to changes in the sample's topography. In this case, the overall imaging rate can be significantly increased, because a general scan speed trade-off decision is not needed and smooth areas are scanned fast. In contrast to methods trying to increase the z-piezo bandwidth, our method is a comparably simple approach that can be easily adapted to standard systems.
Fast and Adaptive Lossless On-Board Hyperspectral Data Compression System for Space Applications
Aranki, Nazeeh; Bakhshi, Alireza; Keymeulen, Didier; Klimesh, Matthew
2009-01-01
Efficient on-board lossless hyperspectral data compression reduces the data volume necessary to meet NASA and DoD limited downlink capabilities. The techniques also improves signature extraction, object recognition and feature classification capabilities by providing exact reconstructed data on constrained downlink resources. At JPL a novel, adaptive and predictive technique for lossless compression of hyperspectral data was recently developed. This technique uses an adaptive filtering method and achieves a combination of low complexity and compression effectiveness that far exceeds state-of-the-art techniques currently in use. The JPL-developed 'Fast Lossless' algorithm requires no training data or other specific information about the nature of the spectral bands for a fixed instrument dynamic range. It is of low computational complexity and thus well-suited for implementation in hardware, which makes it practical for flight implementations of pushbroom instruments. A prototype of the compressor (and decompressor) of the algorithm is available in software, but this implementation may not meet speed and real-time requirements of some space applications. Hardware acceleration provides performance improvements of 10x-100x vs. the software implementation (about 1M samples/sec on a Pentium IV machine). This paper describes a hardware implementation of the JPL-developed 'Fast Lossless' compression algorithm on a Field Programmable Gate Array (FPGA). The FPGA implementation targets the current state of the art FPGAs (Xilinx Virtex IV and V families) and compresses one sample every clock cycle to provide a fast and practical real-time solution for Space applications.
Atomic Forces for Geometry-Dependent Point Multipole and Gaussian Multipole Models
Elking, Dennis M.; Perera, Lalith; Duke, Robert; Darden, Thomas; Pedersen, Lee G.
2010-01-01
In standard treatments of atomic multipole models, interaction energies, total molecular forces, and total molecular torques are given for multipolar interactions between rigid molecules. However, if the molecules are assumed to be flexible, two additional multipolar atomic forces arise due to 1) the transfer of torque between neighboring atoms, and 2) the dependence of multipole moment on internal geometry (bond lengths, bond angles, etc.) for geometry-dependent multipole models. In the curr...
Fast Adapting Ensemble: A New Algorithm for Mining Data Streams with Concept Drift
Agustín Ortíz Díaz
2015-01-01
Full Text Available The treatment of large data streams in the presence of concept drifts is one of the main challenges in the field of data mining, particularly when the algorithms have to deal with concepts that disappear and then reappear. This paper presents a new algorithm, called Fast Adapting Ensemble (FAE, which adapts very quickly to both abrupt and gradual concept drifts, and has been specifically designed to deal with recurring concepts. FAE processes the learning examples in blocks of the same size, but it does not have to wait for the batch to be complete in order to adapt its base classification mechanism. FAE incorporates a drift detector to improve the handling of abrupt concept drifts and stores a set of inactive classifiers that represent old concepts, which are activated very quickly when these concepts reappear. We compare our new algorithm with various well-known learning algorithms, taking into account, common benchmark datasets. The experiments show promising results from the proposed algorithm (regarding accuracy and runtime, handling different types of concept drifts.
Arne Schroeder
2012-08-01
Full Text Available This paper proposes a compression of far field matrices in the fast multipole method and its multilevel extension for electromagnetic problems. The compression is based on a spherical harmonic representation of radiation patterns in conjunction with a radiating mode expression of the surface current. The method is applied to study near field effects and the far field of an antenna placed on a ship surface. Furthermore, the electromagnetic scattering of an electrically large plate is investigated. It is demonstrated, that the proposed technique leads to a significant memory saving, making multipole algorithms even more efficient without compromising the accuracy.
Multipole surface plasmons in metallic nanohole arrays
Nishida, Munehiro; Kadoya, Yutaka
2015-01-01
The quasi-bound electromagnetic modes for the arrays of nanoholes perforated in thin gold film are analyzed both numerically by the rigorous coupled wave analysis (RCWA) method and semi-analytically by the coupled mode method. It is shown that when the size of the nanohole occupies large portion of the unit cell, the surface plasmon polaritons (SPPs) at both sides of the film are combined by the higher order waveguide modes of the holes to produce multipole surface plasmons: coupled surface plasmon modes with multipole texture on the electric field distributions. Further, it is revealed that the multipole texture either enhances or suppresses the couplings between SPPs depending on their diffraction orders and also causes band inversion and reconstruction in the coupled SPP band structure. Due to the multipole nature of the quasi-bound modes, multiple dark modes coexist to produce variety of Fano resonance structures on the transmission and reflection spectra.
Fast binarized time-reversed adapted-perturbation (b-TRAP) optical focusing inside scattering media
Ma, Cheng; Liu, Yan; Wang, Lihong V
2015-01-01
Light scattering inhibits high-resolution optical imaging, manipulation and therapy deep inside biological tissue by preventing focusing. To form deep foci, wavefront-shaping and time-reversal techniques that break the optical diffusion limit have been developed. For in vivo applications, such focusing must provide high gain, high speed, and a large number of spatial modes. However, none of the previous techniques meet these requirements simultaneously. Here, we overcome this challenge by rapidly measuring the perturbed optical field within a single camera exposure followed by adaptively time-reversing the phase-binarized perturbation. Consequently, a phase-conjugated wavefront is synthesized within a millisecond, two orders of magnitude shorter than the digitally achieved record. We demonstrated real-time focusing in dynamic scattering media, and extended laser speckle contrast imaging to new depths. The unprecedented combination of fast response, high gain, and large mode count makes this work a major strid...
Fast-adaptive fiber-optic sensor for ultra-small vibration and deformation measurement
Adaptive fiber-optic interferometer measuring system based on a dynamic hologram recorded in photorefractive CdTe crystal without applying an external electric field is developed. Vectorial mixing of two waves with different polarizations in the anisotropic diffraction geometry allows for the realization of linear regime of phase demodulation at the diffusion hologram. High sensitivity of the interferometer is achieved due to recording of the hologram in reflection geometry at high spatial frequencies in a crystal with sufficient concentration of photorefractive centers. The sensitivity obtained makes possible a broadband detection of ultra-small vibrations with amplitude of less then 0.1 nm. High cut-off frequency of the interferometer achieved using low-power light sources due to fast response of CdTe crystal allows one to eliminate temperature fluctuations and other industrial noises
Fast simulation of transport and adaptive permeability estimation in porous media
Berre, Inga
2005-07-01
The focus of the thesis is twofold: Both fast simulation of transport in porous media and adaptive estimation of permeability are considered. A short introduction that motivates the work on these topics is given in Chapter 1. In Chapter 2, the governing equations for one- and two-phase flow in porous media are presented. Overall numerical solution strategies for the two-phase flow model are also discussed briefly. The concepts of streamlines and time-of-flight are introduced in Chapter 3. Methods for computing streamlines and time-of-flight are also presented in this chapter. Subsequently, in Chapters 4 and 5, the focus is on simulation of transport in a time-of-flight perspective. In Chapter 4, transport of fluids along streamlines is considered. Chapter 5 introduces a different viewpoint based on the evolution of isocontours of the fluid saturation. While the first chapters focus on the forward problem, which consists in solving a mathematical model given the reservoir parameters, Chapters 6, 7 and 8 are devoted to the inverse problem of permeability estimation. An introduction to the problem of identifying spatial variability in reservoir permeability by inversion of dynamic production data is given in Chapter 6. In Chapter 7, adaptive multiscale strategies for permeability estimation are discussed. Subsequently, Chapter 8 presents a level-set approach for improving piecewise constant permeability representations. Finally, Chapter 9 summarizes the results obtained in the thesis; in addition, the chapter gives some recommendations and suggests directions for future work. Part II In Part II, the following papers are included in the order they were completed: Paper A: A Streamline Front Tracking Method for Two- and Three-Phase Flow Including Capillary Forces. I. Berre, H. K. Dahle, K. H. Karlsen, and H. F. Nordhaug. In Fluid flow and transport in porous media: mathematical and numerical treatment (South Hadley, MA, 2001), volume 295 of Contemp. Math., pages 49
Fast Adaptive S-ALOHA Scheme for Event-driven Machine-to-Machine Communications
Wu, Huasen; La, Richard J; Liu, Xin; Zhang, Youguang
2012-01-01
Machine-to-Machine (M2M) communication is now playing a market-changing role in a wide range of business world. However, in event-driven M2M communications, a large number of devices activate within a short period of time, which in turn causes high radio congestions and severe access delay. To solve this problem, we propose a Fast Adaptive S-ALOHA (FASA) scheme for M2M communication systems with bursty traffic. The statistics of consecutive idle and collision slots are used in FASA to accelerate the tracking process of network status, instead of observing in a single slot. Furthermore, the fast convergence property of FASA is guaranteed by using drift analysis. Simulation results demonstrate that the proposed FASA scheme achieves near-optimal performance in reducing access delay, which outperforms that of traditional additive schemes such as PB-ALOHA. Moreover, compared to multiplicative schemes, FASA shows its robust performance even under heavy traffic load in addition to better delay performance.
Park, Eunjeong; Nam, Hyo Suk
2009-11-01
The need for fast treatment of patients in critical conditions motivates the use of mobile devices to provide prompt and consistent communication between hospitals and physicians. We propose a framework that supports ubiquitous access to medical systems using personalized mobile services and integrated medical systems. The proposed service-oriented medical framework provides dynamically composed services that are adapted to contextual variables such as the user's role, the network bandwidth, and resources available at mobile devices while supporting task allocation in distributed servers for massive resource-consuming services. It also manages accurate patient data by integrating local medical systems using medical information standards such as Digital Imaging and Communications in Medicine and Health Level 7. We have demonstrated the effectiveness of our framework by building a prototype of context-based adaptation of computerized tomography image retrieval for acute stroke treatments, which allows images to be viewed on mobile devices with WiMax wireless network. The proposed medical framework reduces hospital delays of patients and facilitates treatments in the absence of medical specialists. PMID:19775976
Fast Model Adaptation for Automated Section Classification in Electronic Medical Records.
Ni, Jian; Delaney, Brian; Florian, Radu
2015-01-01
Medical information extraction is the automatic extraction of structured information from electronic medical records, where such information can be used for improving healthcare processes and medical decision making. In this paper, we study one important medical information extraction task called section classification. The objective of section classification is to automatically identify sections in a medical document and classify them into one of the pre-defined section types. Training section classification models typically requires large amounts of human labeled training data to achieve high accuracy. Annotating institution-specific data, however, can be both expensive and time-consuming; which poses a big hurdle for adapting a section classification model to new medical institutions. In this paper, we apply two advanced machine learning techniques, active learning and distant supervision, to reduce annotation cost and achieve fast model adaptation for automated section classification in electronic medical records. Our experiment results show that active learning reduces the annotation cost and time by more than 50%, and distant supervision can achieve good model accuracy using weakly labeled training data only. PMID:26262005
GPU-based ultra-fast direct aperture optimization for online adaptive radiation therapy
Online adaptive radiation therapy (ART) has great promise to significantly reduce normal tissue toxicity and/or improve tumor control through real-time treatment adaptations based on the current patient anatomy. However, the major technical obstacle for clinical realization of online ART, namely the inability to achieve real-time efficiency in treatment re-planning, has yet to be solved. To overcome this challenge, this paper presents our work on the implementation of an intensity-modulated radiation therapy (IMRT) direct aperture optimization (DAO) algorithm on the graphics processing unit (GPU) based on our previous work on the CPU. We formulate the DAO problem as a large-scale convex programming problem, and use an exact method called the column generation approach to deal with its extremely large dimensionality on the GPU. Five 9-field prostate and five 5-field head-and-neck IMRT clinical cases with 5 x 5 mm2 beamlet size and 2.5 x 2.5 x 2.5 mm3 voxel size were tested to evaluate our algorithm on the GPU. It takes only 0.7-3.8 s for our implementation to generate high-quality treatment plans on an NVIDIA Tesla C1060 GPU card. Our work has therefore solved a major problem in developing ultra-fast (re-)planning technologies for online ART.
Ren, Shenghan; Hu, Haihong; Li, Gen; Cao, Xu; Zhu, Shouping; Chen, Xueli; Liang, Jimin
2016-04-01
Bioluminescence tomography (BLT) has been a valuable optical molecular imaging technique to non-invasively depict the cellular and molecular processes in living animals with high sensitivity and specificity. Due to the inherent ill-posedness of BLT, a priori information of anatomical structure is usually incorporated into the reconstruction. The structural information is usually provided by computed tomography (CT) or magnetic resonance imaging (MRI). In order to obtain better quantitative results, BLT reconstruction with heterogeneous tissues needs to segment the internal organs and discretize them into meshes with the finite element method (FEM). It is time-consuming and difficult to handle the segmentation and discretization problems. In this paper, we present a fast reconstruction method for BLT based on multi-atlas registration and adaptive voxel discretization to relieve the complicated data processing procedure involved in the hybrid BLT/CT system. A multi-atlas registration method is first adopted to estimate the internal organ distribution of the imaged animal. Then, the animal volume is adaptively discretized into hexahedral voxels, which are fed into FEM for the following BLT reconstruction. The proposed method is validated in both numerical simulation and an in vivo study. The results demonstrate that the proposed method can reconstruct the bioluminescence source efficiently with satisfactory accuracy. PMID:27446674
Adaptive control of piezoelectric fast steering mirror for high precision tracking application
A piezoelectric fast steering mirror (PFSM) is a complex, strong coupling nonlinear system that integrates optics, mechanics, electrics, and control. Due to the existence of hysteresis nonlinearity, mechanical resonance, and all kinds of disturbances, precise tracking control of a PFSM is a challenging task. This paper presents a comprehensive study of modeling, controller design, and simulation evaluation for a PFSM system. First a general model of a PFSM system integrating mechanical dynamics, electrical dynamics, and hysteresis nonlinearity is proposed, and then a robust adaptive controller is developed under both unknown hysteresis nonlinearities and parameter uncertainties. The parameters needed directly in the formulation of the controller are adaptively estimated. The proposed control law ensures the uniform boundedness of all signals in the closed-loop system. Furthermore, a stability analysis of the control system is performed to guarantee that the output tracking error converges to zero asymptotically. Finally, simulation tests with different motion trajectories are conducted to verify the effectiveness of the proposed method. (paper)
Ren, Shenghan; Hu, Haihong; Li, Gen; Cao, Xu; Zhu, Shouping; Chen, Xueli; Liang, Jimin
2016-01-01
Bioluminescence tomography (BLT) has been a valuable optical molecular imaging technique to non-invasively depict the cellular and molecular processes in living animals with high sensitivity and specificity. Due to the inherent ill-posedness of BLT, a priori information of anatomical structure is usually incorporated into the reconstruction. The structural information is usually provided by computed tomography (CT) or magnetic resonance imaging (MRI). In order to obtain better quantitative results, BLT reconstruction with heterogeneous tissues needs to segment the internal organs and discretize them into meshes with the finite element method (FEM). It is time-consuming and difficult to handle the segmentation and discretization problems. In this paper, we present a fast reconstruction method for BLT based on multi-atlas registration and adaptive voxel discretization to relieve the complicated data processing procedure involved in the hybrid BLT/CT system. A multi-atlas registration method is first adopted to estimate the internal organ distribution of the imaged animal. Then, the animal volume is adaptively discretized into hexahedral voxels, which are fed into FEM for the following BLT reconstruction. The proposed method is validated in both numerical simulation and an in vivo study. The results demonstrate that the proposed method can reconstruct the bioluminescence source efficiently with satisfactory accuracy.
Suppressing CMB low multipoles with ISW effect
Recent results of Planck data reveal that the power [1,2] in the low multipoles of the CMB angular power spectrum, approximately up to l = 30, is significantly lower than the theoretically predicted in the best fit ΛCDM model. There are different known physical effects that can affect the power at low multipoles, such as features in the primordial power spectrum (PPS) in some models of inflation and ISW effect. In this paper we investigate the possibility of invoking the Integrated Sachs-Wolfe (ISW) effect to explain the power deficit at low multipoles. The ISW effect that originates from the late time expansion history of the universe is rich in possibilities given the limited understanding of the origin of dark energy (DE). It is a common understanding that the ISW effect adds to the power at the low multipoles of the CMB angular power spectrum. In this paper we carry out an analytic study to show that there are some expansion histories in which the ISW effect, instead of adding power, provides negative contribution to the power at low multipoles. Guided by the analytic study, we present examples of the features required in the late time expansion history of the universe that could explain the power deficiency through the ISW effect. We also show that an ISW origin of power deficiency is consistent, at present, with other cosmological observations that probe the expansion history such as distance modulus, matter power spectrum and the evolution of cluster number count. We also show that the ISW effect may be distinguished from power deficit originating from features in the PPS using the measurements of the CMB polarization spectrum at low multipoles expected from Planck. We conclude that the power at low multipoles of the CMB anisotropy could well be closely linked to Dark Energy puzzle in cosmology and this observation could be actually pointing to richer phenomenology of DE beyond the cosmological constant Λ
Eddy-current-induced multipole field calculations
Time-varying magnetic fields of magnets in booster accelerators induce substantial eddy currents in the vacuum chambers. The eddy currents in turn act to produce various multipole fields that act on the beam. These fields must be taken into account when doing a lattice design. In the APS booster, the relatively long dipole magnets (3 meters) are linearly ramped to accelerate the injected 325 MeV beam to 7 GeV. Substantial dipole and sextupole fields are generated in the elliptical vacuum chamber from the induced eddy currents. In this note, formulas for the induced dipole and sextupole fields are derived for elliptical and rectangular vacuum chambers for a time-varying dipole field. A discussion is given on how to generalize this derivation method to include eddy-current-induced multipole fields from higher multipole magnets (quadrupole, sextupole, etc.). Finally, transient effects are considered
Electromagnetic multipole fields of neutron stars
There is now indisputable evidence that some pulsars possess space velocities so high that internal asymmetries in the dynamics of their formation are strongly implied. We develop in this paper a complete formalism for the calculation of the only such mechanism that has yet been subjected to quantitative analysis: electromagnetic recoil radiation. To make the general problem tractable without doing violence to the physics, we have made the following simplifying assumptions: (1) the magnetic induction B in athin shell enclosing the surface can be satisfactorily approximated by a sum of vacuum multipole fields; (2) the star is spherical, and all parts are in good electrical contact; (3) vertical-bar Ω X r vertical-barvery-much-less-thanc everywhere within the star; and (4) the star is surrounded by a vacuum. Our qualitative conclusions hold even if these assumptions are violated, but corrections to our quantitative results required by a relaxation of our assumptions are not easily computed.Given this simple electrodynamic model of a neutron star, we solve the following problems: (1) What electric multipoles are induced by each magnetic multipole. (2) What is the general formula for the recoil produced by the projection on the rotational axis of a net linear momentum flux produced by the rotation of any two magnetic multipoles. (3) What is the set of centered multipoles that represents the field of an arbitrary off-centered multipole. We use these general results go perform a detailed analysis of the linear momentum radiated by an off-centered dipole. We find a force larger by a factor 6 than that obtained for the special case treated in the best previous calculation. In spite of this considerable increase in the computed strengrh of the effect, we still believe it to be too weak to produce the large space velocities observed for pulsars. For the mechanism to be effective, the pulsar must be born rotating near the breakup velocity
A Fast Adaptive Artificial Neural Network Controller for Flexible Link Manipulators
Amin Riad Maouche
2016-01-01
Full Text Available This paper describes a hybrid approach to the problem of controlling flexible link manipulators in the dynamic phase of the trajectory. A flexible beam/arm is an appealing option for civil and military applications, such as space-based robot manipulators. However, flexibility brings with it unwanted oscillations and severe chattering which may even lead to an unstable system. To tackle these challenges, a novel control architecture scheme is presented. First, a neural network controller based on the robot’s dynamic equation of motion is elaborated. Its aim is to produce a fast and stable control of the joint position and velocity and damp the vibration of each arm. Then, an adaptive Cerebellar Model Articulation Controller (CMAC is implemented to balance unmodeled dynamics, enhancing the precision of the control. Efficiency of the new controller obtained is tested on a two-link flexible manipulator. Simulation results on a dynamic trajectory with a sinusoidal form show the effectiveness of the proposed control strategy.
Adaptation of fast responding power supply for radial position control in SST-1
A high current, fast responding power supply was installed in 2005 for vertical stabilization of elongated plasmas in SST-1 tokamak. Presently, during initial experiments of SST-1 tokamak the need for radial control during current build-up was envisaged. For this purpose the existing power supply was suitable and the same was re-commissioned and control adaptations were carried as per experimental requirements. This paper highlights the capabilities of the power supply and details the modifications in the control interfaces and test programs for the radial control purpose. Details of the operation of the power supply along with control interfaces with performance measurements are provided. The re-commissioning provided an opportunity in the trouble shooting methods and sequential operation of the system. With the operational use on the actual coil the mutual effects are understood better and appropriate test programs are prepared. The power supply provided satisfactory performance for the intended use. In additional the system is suitable to simulate a plasma current loop to enable the testing and calibration of Rogowski coil used for plasma current measurement. (author)
Zuidmeer-Jongejan, Laurian; Fernandez-Rivas, Montserrat; Poulsen, Lars K.;
2012-01-01
ABSTRACT: The FAST project (Food Allergy Specific Immunotherapy) aims at the development of safe and effective treatment of food allergies, targeting prevalent, persistent and severe allergy to fish and peach. Classical allergen-specific immunotherapy (SIT), using subcutaneous injections with...... aqueous food extracts may be effective but has proven to be accompanied by too many anaphylactic side-effects. FAST aims to develop a safe alternative by replacing food extracts with hypoallergenic recombinant major allergens as the active ingredients of SIT. Both severe fish and peach allergy are caused...... lower their threshold for fish or peach intake, thereby decreasing their anxiety and dependence on rescue medication....
Dipole leakage and low CMB multipoles
A number of studies of WMAP-7 have highlighted that the power at the low multipoles in CMB power spectrum is lower than their theoretically predicted value. Angular correlation between the orientations of these low multipoles have also been discovered. While these observations may have cosmological ramification, it is important to investigate possible observational artifacts that can mimic them. The CMB dipole, which is much higher than the quadrupole, can get leaked to the higher multipoles due to the non-circular beam of the CMB experiment. In this paper, an analytical method has been developed and simulations are carried out to study the effect of the non-circular beam on power leakage from the dipole. It has been shown that the small, but non-negligible power from the dipole can get transferred to the quadrupole and the higher multipoles due to the non-circular beam. Simulations have also been carried out for Planck scan strategy, and comparative results between WMAP and Planck have been presented
Multipole expansion approach to Ostwald ripening
A multipole expansion theory is used to derive a kinetic equation for a collection of droplets undergoing the Ostwald ripening or coarsening in two dimensions. The theory allows the morphological change of coarsening droplets as well as the migration behavior. Numerical computations are performed to give a good estimate to the solution of Laplace equation even in a few mode truncation
Dipole leakage and low CMB multipoles
Das, Santanu
2012-01-01
A number of studies of WMAP-7 have highlighted that the power at the low multipoles in CMB power spectrum are lower than their theoretically predicted values. Angular correlation between the orientations of these low multipoles have also been discovered. While these observations may have cosmological ramification, it is important to investigate possible observational artifacts that can mimic them. The CMB dipole which is almost 550 times higher than the quadrupole can get leaked to the higher multipoles due to the non-circular beam of the CMB experiment. In this paper an analytical method has been developed and simulations are carried out to study the effect of the non-circular beam on power leakage from the dipole. It has been shown that the small, but non-negligible power from the dipole can get transferred to the quadrupole and the higher multipoles due to the non-circular beam. Simulations have also been carried out for Planck scan strategy and comparative results between WMAP and Planck have been present...
Thermal conductivity and resonant multipole interactions.
Nyeland, C.; Mason, E. A.; Monchick, L.
1972-01-01
Investigation of the influence of resonant multipole interactions on exchanges of rotational energy in molecular collisions, by means of a simple two-state impact-parameter approximation. It is found that dipole-quadrupole and quadrupole-quadrupole interactions can have a significant effect for molecules with low moments of inertia.
Equilibrium and stability of the multipole
The equilibrium and linear ballooning mode stability of a plasma in an axisymmetric closed field line device such as the multipole is investigated. Two models of the plasma are used, a kinetic model and an ideal MHD model, and the results are compared. Numerical calculations are made of the equilibrium and of ballooning mode stability criteria in the Wisconsin Levitated Octupole device
Dailey, Rachael E; Fontaine, Christine M; Avery, Julie P
2016-09-01
Most organisms undergo changes in their environment, both predictably and unpredictably, which require them to alter priorities in nutrient allocation with regards to food availability. Species that more predictably encounter extended periods of limited food resources or intake while mitigating the negative effects of starvation are considered to be fasting adapted. Northern elephant seals (NES) are one such species and routinely undergo extended periods of fasting for breeding, molting, as well as a post-weaning fast at 6-8weeks of age. However, during unusual times of nutritional deprivation, animals may enter stage III fasting. While fasting and foraging in this species has been extensively studied, realimentation following fasting beyond normal life history parameters has not been investigated. In this study, changes in ghrelin, growth hormone (GH), and insulin-like growth factor (IGF)-I were compared across 8weeks of realimentation following emaciation in three age classes: neonates, post-molt pups, and yearlings. Longitudinal changes in hormone profiles indicate that neonate and post-molt pups are slow to recover mass and positive energy balance despite an energy dense diet fed at 10% body mass. In addition, ghrelin and GH concentrations remained elevated in post-molt pups compared to other age classes. Changes in hormone concentrations early in realimentation indicate that yearling animals recover more rapidly from periods of nutritional deprivation than do younger animals. Overall, this suggests that the ability to regulate metabolic homeostasis with regards to nutrient allocation may develop over time, even in a species that is considered to be fasting adapted. PMID:27288636
HU Bing; FANG Zhi-chu
2008-01-01
Due to actuator time delay existing in an adaptive control of the active balancing system for a fastspeed-varying Jeffcott rotor, if an unsynchronized control force (correction imbalance) is applied to the system,it may lead to degradation in control efficiency and instability of the control system. In order to avoid theseshortcomings, a simple adaptive controller was designed for a strictly positive real rotor system with actuatortime delay, then a Lyapunov-Krasovskii functional was constructed after an appropriate transform of this sys-tem model, the stability conditions of this adaptive control system with actuator time delay were derived. Afteradding a filter function, the active balancing system for the fast speed-varying Jeffcott rotor with actuator timedelay can easily be converted to a strictly positive real system, and thus it can use the above adaptive controllersatisfying the stability conditions. Finally, numerical simulations show that the adaptive controller proposedworks very well to perform the active balancing for the fast speed-varying Jeffcott rotor with actuator timedelay.
Multipole Analysis of Circular Cylindircal Magnetic Systems
J Selvaggi
2006-01-09
This thesis deals with an alternate method for computing the external magnetic field from a circular cylindrical magnetic source. The primary objective is to characterize the magnetic source in terms of its equivalent multipole distribution. This multipole distribution must be valid at points close to the cylindrical source and a spherical multipole expansion is ill-equipped to handle this problem; therefore a new method must be introduced. This method, based upon the free-space Green's function in cylindrical coordinates, is developed as an alternative to the more familiar spherical harmonic expansion. A family of special functions, called the toroidal functions or Q-functions, are found to exhibit the necessary properties for analyzing circular cylindrical geometries. In particular, the toroidal function of zeroth order, which comes from the integral formulation of the free-space Green's function in cylindrical coordinates, is employed to handle magnetic sources which exhibit circular cylindrical symmetry. The toroidal functions, also called Q-functions, are the weighting coefficients in a ''Fourier series-like'' expansion which represents the free-space Green's function. It is also called a toroidal expansion. This expansion can be directly employed in electrostatic, magnetostatic, and electrodynamic problems which exhibit cylindrical symmetry. Also, it is shown that they can be used as an alternative to the Elliptic integral formulation. In fact, anywhere that an Elliptic integral appears, one can replace it with its corresponding Q-function representation. A number of problems, using the toroidal expansion formulation, are analyzed and compared to existing known methods in order to validate the results. Also, the equivalent multipole distribution is found for most of the solved problems along with its corresponding physical interpretation. The main application is to characterize the external magnetic field due to a six
Drumea, Georges,; Frezza-Buet, Hervé
2007-01-01
International audience In this paper, an original method named GNG-T, extended from GNG-U algorithm by Fritzke is presented. The method performs continuously vector quantization over a distribution that changes over time. It deals with both sudden changes and continuous ones, and is thus suited for video tracking framework, where continuous tracking is required as well as fast adaptation to incoming and outgoing people. The central mechanism relies on the management of quantization resolut...
Nuclear photonics at ultra-high counting rates and higher multipole excitations
Thirolf, P. G.; Habs, D.; Filipescu, D.; Gernhäuser, R.; Günther, M. M.; Jentschel, M.; Marginean, N.; Pietralla, N.
2012-07-01
Next-generation γ beams from laser Compton-backscattering facilities like ELI-NP (Bucharest)] or MEGa-Ray (Livermore) will drastically exceed the photon flux presently available at existing facilities, reaching or even exceeding 1013 γ/sec. The beam structure as presently foreseen for MEGa-Ray and ELI-NP builds upon a structure of macro-pulses (˜120 Hz) for the electron beam, accelerated with X-band technology at 11.5 GHz, resulting in a micro structure of 87 ps distance between the electron pulses acting as mirrors for a counterpropagating intense laser. In total each 8.3 ms a γ pulse series with a duration of about 100 ns will impinge on the target, resulting in an instantaneous photon flux of about 1018 γ/s, thus introducing major challenges in view of pile-up. Novel γ optics will be applied to monochromatize the γ beam to ultimately ΔE/E˜10-6. Thus level-selective spectroscopy of higher multipole excitations will become accessible with good contrast for the first time. Fast responding γ detectors, e.g. based on advanced scintillator technology (e.g. LaBr3(Ce)) allow for measurements with count rates as high as 106-107 γ/s without significant drop of performance. Data handling adapted to the beam conditions could be performed by fast digitizing electronics, able to sample data traces during the micro-pulse duration, while the subsequent macro-pulse gap of ca. 8 ms leaves ample time for data readout. A ball of LaBr3 detectors with digital readout appears to best suited for this novel type of nuclear photonics at ultra-high counting rates.
Nuclear photonics at ultra-high counting rates and higher multipole excitations
Thirolf, P. G.; Habs, D.; Filipescu, D.; Gernhaeuser, R.; Guenther, M. M.; Jentschel, M.; Marginean, N.; Pietralla, N. [Fakultaet f. Physik, Ludwig-Maximilians-Universitaet Muenchen, Garching (Germany); Fakultaet f. Physik, Ludwig-Maximilians-Universitaet Muenchen, Garching, Germany and Max-Planck-Institute f. Quantum Optics, Garching (Germany); IFIN-HH, Bucharest-Magurele (Romania); Physik Department E12,Technische Universitaet Muenchen, Garching (Germany); Max-Planck-Institute f. Quantum Optics, Garching (Germany); Institut Laue-Langevin, Grenoble (France); Physik Department E12,Technische Universitaet Muenchen, Garching (Germany); Institut f. Kernphysik, Technische Universitaet Darmstadt (Germany)
2012-07-09
Next-generation {gamma} beams from laser Compton-backscattering facilities like ELI-NP (Bucharest)] or MEGa-Ray (Livermore) will drastically exceed the photon flux presently available at existing facilities, reaching or even exceeding 10{sup 13}{gamma}/sec. The beam structure as presently foreseen for MEGa-Ray and ELI-NP builds upon a structure of macro-pulses ({approx}120 Hz) for the electron beam, accelerated with X-band technology at 11.5 GHz, resulting in a micro structure of 87 ps distance between the electron pulses acting as mirrors for a counterpropagating intense laser. In total each 8.3 ms a {gamma} pulse series with a duration of about 100 ns will impinge on the target, resulting in an instantaneous photon flux of about 10{sup 18}{gamma}/s, thus introducing major challenges in view of pile-up. Novel {gamma} optics will be applied to monochromatize the {gamma} beam to ultimately {Delta}E/E{approx}10{sup -6}. Thus level-selective spectroscopy of higher multipole excitations will become accessible with good contrast for the first time. Fast responding {gamma} detectors, e.g. based on advanced scintillator technology (e.g. LaBr{sub 3}(Ce)) allow for measurements with count rates as high as 10{sup 6}-10{sup 7}{gamma}/s without significant drop of performance. Data handling adapted to the beam conditions could be performed by fast digitizing electronics, able to sample data traces during the micro-pulse duration, while the subsequent macro-pulse gap of ca. 8 ms leaves ample time for data readout. A ball of LaBr{sub 3} detectors with digital readout appears to best suited for this novel type of nuclear photonics at ultra-high counting rates.
Fast 3D EM scattering and radiation solvers based on MLFMA
Hu Jun; Nie Zaiping; Lei Lin; Hu Jie; Gong Xiaodong; Zhao Huapeng
2008-01-01
As the fastest integral equation solver to date, the multilevel fast multipole algorithm (MLFMA)has been applied successfully to solve electromagnetic scattering and radiation from 3D electrically large objects.But for very large-scale problems, the storage and CPU time required in MLFMA are still expensive. Fast 3D electromagnetic scattering and radiation solvers are introduced based on MLFMA. A brief review of MLFMA is first given. Then, four fast methods including higher-order MLFMA (HO-MLFMA), fast far field approximation combined with adaptive ray propagation MLFMA (FAFFA-ARP-MLFMA), local MLFMA and parallel MLFMA are introduced. Some typical numerical results demonstrate the efficiency of these fast methods.
A Fast, Locally Adaptive, Interactive Retrieval Algorithm for the Analysis of DIAL Measurements
Samarov, D. V.; Rogers, R.; Hair, J. W.; Douglass, K. O.; Plusquellic, D.
2010-12-01
Differential absorption light detection and ranging (DIAL) is a laser-based tool which is used for remote, range-resolved measurement of particular gases in the atmosphere, such as carbon-dioxide and methane. In many instances it is of interest to study how these gases are distributed over a region such as a landfill, factory, or farm. While a single DIAL measurement only tells us about the distribution of a gas along a single path, a sequence of consecutive measurements provides us with information on how that gas is distributed over a region, making DIAL a natural choice for such studies. DIAL measurements present a number of interesting challenges; first, in order to convert the raw data to concentration it is necessary to estimate the derivative along the path of the measurement. Second, as the distribution of gases across a region can be highly heterogeneous it is important that the spatial nature of the measurements be taken into account. Finally, since it is common for the set of collected measurements to be quite large it is important for the method to be computationally efficient. Existing work based on Local Polynomial Regression (LPR) has been developed which addresses the first two issues, but the issue of computational speed remains an open problem. In addition to the latter, another desirable property is to allow user input into the algorithm. In this talk we present a novel method based on LPR which utilizes a variant of the RODEO algorithm to provide a fast, locally adaptive and interactive approach to the analysis of DIAL measurements. This methodology is motivated by and applied to several simulated examples and a study out of NASA Langley Research Center (LaRC) looking at the estimation of aerosol extinction in the atmosphere. A comparison study of our method against several other algorithms is also presented. References Chaudhuri, P., Marron, J.S., Scale-space view of curve estimation, Annals of Statistics 28 (2000) 408-428. Duong, T., Cowling
Theory and application of plane elliptic multipoles for static magnetic fields
Schnizer, P. [Gesellschaft fuer Schwerionenforschung, Darmstadt (Germany)], E-mail: p.schnizer@gsi.de; Schnizer, B. [Technische Universitaet Graz (Austria)], E-mail: schnizer@itp.tu-graz.ac.at; Akishin, P. [Joint Institute for Nuclear Research, Dubna, Moscow Region (Russian Federation)], E-mail: akishin@jinr.ru; Fischer, E. [Gesellschaft fuer Schwerionenforschung, Darmstadt (Germany)], E-mail: e.fischer@gsi.de
2009-08-21
Standard textbooks on beam dynamics study the impact of the magnetic field quality on the beam using field representations based on circular multipoles. Iron dominated magnets, however, typically provide a good field region with a non-circular aspect ratio (i.e. an ellipse whose axis a is significantly larger than the axis b); a boundary not ideal for circular multipoles. The development of superconductors, originally driven to reach fields above {approx}2T, allows using them today in completely different fields: iron dominated DC magnets, to save the energy for coil powering as well as repeatedly fast ramped magnets. The cold mass of magnets, housed in common cryostats sectors, makes it tedious to implement additional correction magnets at a later stage, as it requires to warm up the sections where the magnets should be installed as well as unwelding the cryostat. Thus the field homogeneity of the magnets and its influence on the beam has to be thoroughly studied during the project planning phase. Elliptic multipoles, a new type of field expansion for static or quasi-static (here magnetic) two-dimensional fields, are proposed and investigated, which are particular solutions of the potential equation in plane elliptic coordinates obtained by the method of separation. The proper subsets of these particular solutions appropriate for representing static real or complex fields regular within an ellipse are identified. Formulas are given for computing expansion coefficients from given fields. The advantage of this new approach is that the expansion is valid, convergent and accurate in a larger domain, namely in an ellipse circumscribed to the reference circle of the common circular multipoles in polar coordinates. Formulas are derived for calculating the circular multipoles from the elliptical ones. The effectiveness of the approach was tested on many different magnet designs and is illustrated here on the dipole design chosen for the core synchrotron (SIS 100) of the
Multipole Traps as Tools in Environmental Studies
Mihalcea, Bogdan M; Giurgiu, Liviu C; Groza, Andreea; Surmeian, Agavni; Ganciu, Mihai; Filinov, Vladimir; Lapitsky, Dmitry; Deputatova, Lidiya; Vasilyak, Leonid; Pecherkin, Vladimir; Vladimirov, Vladimir; Syrovatka, Roman
2015-01-01
Trapping of microparticles, nanoparticles and aerosols is an issue of major interest for physics and chemistry. We present a setup intended for microparticle trapping in multipole linear Paul trap geometries, operating under Standard Ambient Temperature and Pressure (SATP) conditions. A 16-electrode linear trap geometry has been designed and tested, with an aim to confine a larger number of particles with respect to quadrupole traps and thus enhance the signal to noise ratio, as well as to study microparticle dynamical stability in electrodynamic fields. Experimental tests and numerical simulations suggest that multipole traps are very suited for high precision mass spectrometry measurements in case of different microparticle species or to identify the presence of certain aerosols and polluting agents in the atmosphere. Particle traps represent versatile tools for environment monitoring or for the study of many-body Coulomb systems and dusty plasmas.
Poloidal OHMIC heating in a multipole
Holly, D.J.
1982-01-01
The feasibility of using poloidal currents to heat plasmas confined by a multipole field has been examined experimentaly in Tokapole II. The machine is operated as a toroidal octupole, with a time-varying toroidal magnetic field driving poloidal plasma currents I/sub plasma/ - 20 kA to give densities n/sub e/ - 10/sup 13/ cm/sup -3/ and temperatures T/sub e/ - 30 eV.
Poloidal OHMIC heating in a multipole
The feasibility of using poloidal currents to heat plasmas confined by a multipole field has been examined experimentaly in Tokapole II. The machine is operated as a toroidal octupole, with a time-varying toroidal magnetic field driving poloidal plasma currents I/sub plasma/ - 20 kA to give densities n/sub e/ - 1013 cm-3 and temperatures T/sub e/ - 30 eV
Multipole solutions in metric-affine gravity
Socorro, J; Macías, A; Mielke, E W; Socorro, José; Lämmerzahl, Claus; Macías, Alfredo; Mielke, Eckehard W.
1998-01-01
Above Planck energies, the spacetime might become non--Riemannian, as it is known fron string theory and inflation. Then geometries arise in which nonmetricity and torsion appear as field strengths, side by side with curvature. By gauging the affine group, a metric affine gauge theory emerges as dynamical framework. Here, by using the harmonic map ansatz, a new class of multipole like solutions in the metric affine gravity theory (MAG) is obtained.
A new adaptive fast motion estimation algorithm based on local motion similarity degree (LMSD)
LIU Long; HAN Chongzhao; BAI Yan
2005-01-01
In the motion vector field adaptive search technique (MVFAST) and the predictive motion vector field adaptive search technique (PMVFAST), the size of the largest motion vector from the three adjacent blocks (left, top, top-right) is compared with the threshold to select different search scheme. But a suitable search center and search pattern will not be selected in the adaptive search technique when the adjacent motion vectors are not coherent in local region. This paper presents an efficient adaptive search algorithm. The motion vector variation degree (MVVD) is considered a reasonable factor for adaptive search selection. By the relationship between local motion similarity degree (LMSD) and the variation degree of motion vector (MVVD), the motion vectors are classified as three categories according to corresponding LMSD; then different proposed search schemes are adopted for motion estimation. The experimental results show that the proposed algorithm has a significant computational speedup compared with MVFAST and PMVFAST algorithms, and offers a similar, even better performance.
Alsaadi, Fuad E.
2016-03-01
Optical wireless systems are promising candidates for next-generation indoor communication networks. Optical wireless technology offers freedom from spectrum regulations and, compared to current radio-frequency networks, higher data rates and increased security. This paper presents a fast adaptation method for multibeam angle and delay adaptation systems and a new spot-diffusing geometry, and also considers restrictions needed for complying with eye safety regulations. The fast adaptation algorithm reduces the computational load required to reconfigure the transmitter in the case of transmitter and/or receiver mobility. The beam clustering approach enables the transmitter to assign power to spots within the pixel's field of view (FOV) and increases the number of such spots. Thus, if the power per spot is restricted to comply with eye safety standards, the new approach, in which more spots are visible within the FOV of the pixel, leads to enhanced signal-to-noise ratio (SNR). Simulation results demonstrate that the techniques proposed in this paper lead to SNR improvements that enable reliable operation at data rates as high as 15 Gbit/s. These results are based on simulation and not on actual measurements or experiments.
Source integrals of multipole moments for static space-times
Hernandez-Pastora, J L; Ruiz, E
2016-01-01
The definition of Komar for the mass of a relativistic source is used as a starting point to introduce volume integrals for Relativistic Multipole Moments (RMM). A certain generalization of the classical Gauss theorem is used to rewrite these multipole moments as integrals over a surface at the infinity. Therefore it is shown that the above generalization leads to Asymptotic Relativistic Multipole Moments (ARMM), recovering the multipoles of Geroch or Thorne, when the integrals are evaluated in asympotically cartesian harmonic coordinates. Relationships regarding the Thorne definition and the classical theory of moments are shown.
肖立民; 许希斌; 姚彦
2001-01-01
Diversity combining technologies are analyzed for fastfrequency-hopping spread spectrum systems during partial-band noise jamming to develop a novel combining receiver called an Adaptive Threshold Clipper Combining Receiver (ATCCR). The optimal clipping level for an ATCCR is analyzed, computed, and compared with several other diversity combining technologies. Since the ATCCR can estimate the power of the jamming and the number of jammed frequency cells to adaptively adjust the clipper's threshold, the system performance using the adaptive threshold clipper combining technique can be greatly improved.
Fast and Flexible Tracking and Mitigating a Jamming Signal with an Adaptive Notch Filter
BORIO DANIELE; O'DRISCOLL CILLIAN; Fortuny Guasch, Joaquim
2014-01-01
GNSS jammers are small portable devices able to broadcast disruptive interference and overpower the much weaker GNSS signals. The authors consider the use of an adaptive notch filter as an effective solution for mitigating jamming effects.
Fast Spectral Velocity Estimation Using Adaptive Techniques: In-Vivo Results
Gran, Fredrik; Jakobsson, Andreas; Udesen, Jesper; Jensen, Jørgen Arendt
Adaptive spectral estimation techniques are known to provide good spectral resolution and contrast even when the observation window(OW) is very sbort. In this paper two adaptive techniques are tested and compared to the averaged perlodogram (Welch) for blood velocity estimation. The Blood Power...... spectral Capon (BPC) method is based on a standard minimum variance technique adapted to account for both averaging over slowtime and depth. The Blood Amplitude and Phase Estimation technique (BAPES) is based on finding a set of matched filters (one for each velocity component of interest) and filtering...... the blood process over slow-time and averaging over depth to find the power spectral density estimate. In this paper, the two adaptive methods are explained, and performance Is assessed in controlled steady How experiments and in-vivo measurements. The three methods were tested on a circulating How...
Role of PGC-1{alpha} in exercise and fasting induced adaptations in mouse liver
Haase, Tobias Nørresø; Jørgensen, Stine Ringholm; Leick, Lotte;
2011-01-01
The transcriptional coactivator peroxisome proliferator activated receptor (PPAR)-¿ coactivator (PGC)-1a plays a role in regulation of several metabolic pathways. By use of whole body PGC-1a knockout (KO) mice we investigated the role of PGC-1a in fasting, acute exercise and exercise training...... induced regulation of key proteins in gluconeogenesis and metabolism in the liver. In both wild type (WT) and PGC-1a KO mice liver, the mRNA content of the gluconeogenic proteins glucose-6-phosphatase (G6Pase) and phosphoenolpyruvate carboxykinase (PEPCK) was upregulated during fasting. Pyruvate...... carboxylase (PC) remained unchanged after fasting in WT mice, but was upregulated in PGC-1a KO mice. In response to a single exercise bout G6Pase mRNA was upregulated in both genotypes, whereas no significant changes were detected in PEPCK or PC mRNA. While G6Pase and PC protein remained unchanged, liver...
Khosravinia, H
2016-06-01
Four hundred and eighty mixed-sex broiler chicks aged 3 h after hatching were allotted according to a completely random design in a 6 × 2 × 2 factorial schedule into two groups of 12 replications of 20 chicks each. The main experimental factors were fasting for 0, 6, 12, 24, 36 and 48 h after chick placement and calcium gluconate (Ca-glu) injection (0 and 0.6 ml). Live body weight (BW) of chicks decreased linearly (Y = 43.36-0.109BW0 h , r(2) = 0.876) as neonatal fasting extended. Injection of 0.6 ml Ca-glu at 3 h post-hatching did not affect weight loss of chicks. Yolk residuals (YR) utilized linearly (Y = 5.75-0.062YR, r(2) = 0.956) by 0.062 g/h in neonate fasted chicks up to 48 h, showing no effect of Ca-glu injection. Neonatal fasting periods longer than 12 h increased liver weight (p < 0.05). The mean absolute and proportional (% of BW0 h ) breast and leg weight were reduced linearly as neonatal fasting extended (p < 0.05). Serum glucose concentration increased up to 6 h and then reduced linearly to 150 mg/dl after 48-h fasting. The Ca-glu treatment influenced serum glucose level for a short period up to 6 h of fasting. Serum Ca concentration sharply increased up to threefolds in the birds received Ca-glu injection resulting in acute hypercalcemia, then decreased to the initial level after 24-h feed withdrawal (p < 0.05). The mean serum level for creatinine, uric acid, cholesterol, HDL, albumins and total proteins significantly increased during the fasting periods of 6 to 48 h and significantly elevated in the birds receiving 0.6-ml Ca-glu injection compared with the non-treated chicks (p < 0.05). It was concluded that subcutaneous administration of 0.6 ml Ca-glu in the chick's neck did not suitably support the increased metabolic demands for glucose and calcium in feed-deprived neonate chicks. PMID:26344414
Fast Muting Adaptation for LTE-A HetNets with Remote Radio Heads
Soret, Beatriz; Pedersen, Klaus I.; Kolding, Troels Emil; Kroener, Hans; Maniatis, Ioannis
Management in macro-RRH cases is developed. Based on minimal information collection, we propose an algorithm that adjusts the ABS ratio on a fast basis, aiming at balancing the instantaneous load betweeen the macro and the RRH layer. Performance results with bursty traffic and low and high load conditions...
Multipole electrostatics in hydration free energy calculations.
Shi, Yue; Wu, Chuanjie; Ponder, Jay W; Ren, Pengyu
2011-04-15
Hydration free energy (HFE) is generally used for evaluating molecular solubility, which is an important property for pharmaceutical and chemical engineering processes. Accurately predicting HFE is also recognized as one fundamental capability of molecular mechanics force field. Here, we present a systematic investigation on HFE calculations with AMOEBA polarizable force field at various parameterization and simulation conditions. The HFEs of seven small organic molecules have been obtained alchemically using the Bennett Acceptance Ratio method. We have compared two approaches to derive the atomic multipoles from quantum mechanical calculations: one directly from the new distributed multipole analysis and the other involving fitting to the electrostatic potential around the molecules. Wave functions solved at the MP2 level with four basis sets (6-311G*, 6-311++G(2d,2p), cc-pVTZ, and aug-cc-pVTZ) are used to derive the atomic multipoles. HFEs from all four basis sets show a reasonable agreement with experimental data (root mean square error 0.63 kcal/mol for aug-cc-pVTZ). We conclude that aug-cc-pVTZ gives the best performance when used with AMOEBA, and 6-311++G(2d,2p) is comparable but more efficient for larger systems. The results suggest that the inclusion of diffuse basis functions is important for capturing intermolecular interactions. The effect of long-range correction to van der Waals interaction on the hydration free energies is about 0.1 kcal/mol when the cutoff is 12Å, and increases linearly with the number of atoms in the solute/ligand. In addition, we also discussed the results from a hybrid approach that combines polarizable solute with fixed-charge water in the HFE calculation. PMID:20925089
Abdellah, Skoudarli; Mokhtar, Nibouche; Amina, Serir
2015-11-01
The H.264/AVC video coding standard is used in a wide range of applications from video conferencing to high-definition television according to its high compression efficiency. This efficiency is mainly acquired from the newly allowed prediction schemes including variable block modes. However, these schemes require a high complexity to select the optimal mode. Consequently, complexity reduction in the H.264/AVC encoder has recently become a very challenging task in the video compression domain, especially when implementing the encoder in real-time applications. Fast mode decision algorithms play an important role in reducing the overall complexity of the encoder. In this paper, we propose an adaptive fast intermode algorithm based on motion activity, temporal stationarity, and spatial homogeneity. This algorithm predicts the motion activity of the current macroblock from its neighboring blocks and identifies temporal stationary regions and spatially homogeneous regions using adaptive threshold values based on content video features. Extensive experimental work has been done in high profile, and results show that the proposed source-coding algorithm effectively reduces the computational complexity by 53.18% on average compared with the reference software encoder, while maintaining the high-coding efficiency of H.264/AVC by incurring only 0.097 dB in total peak signal-to-noise ratio and 0.228% increment on the total bit rate.
Cordero, Ninoska; Maza, Felipe; Navea-Perez, Helen; Aravena, Andrés; Marquez-Fontt, Bárbara; Navarrete, Paola; Figueroa, Guillermo; González, Mauricio; Latorre, Mauricio; Reyes-Jara, Angélica
2016-01-01
Listeria monocytogenes has become one of the principal foodborne pathogens worldwide. The capacity of this bacterium to grow at low temperatures has opened an interesting field of study in terms of the identification and classification of new strains of L. monocytogenes with different growth capacities at low temperatures. We determined the growth rate at 8°C of 110 strains of L. monocytogenes isolated from different food matrices. We identified a group of slow and fast strains according to their growth rate at 8°C and performed a global transcriptomic assay in strains previously adapted to low temperature. We then identified shared and specific transcriptional mechanisms, metabolic and cellular processes of both groups; bacterial motility was the principal process capable of differentiating the adaptation capacity of L. monocytogenes strains with different ranges of tolerance to low temperatures. Strains belonging to the fast group were less motile, which may allow these strains to achieve a greater rate of proliferation at low temperature. PMID:26973610
IBA calculations of multipole mixing ratios
Multipole mixing ratios of gamma rays serve as a sensitive probe of nuclear structure. Directional measurements, mostly by γγ angular correlation and nuclear orientation, yield the mixing ratios delta with sign. Unfortunately the results are often double valued. Furthermore, comparison between experimental and theoretical conversion coefficients can give |delta|. The E2/M1 mixing ratio delta(E2/M1) is of special importance in studying collective phenomena and models for their description. Present study is carried out mostly in terms of IBA-1, which makes no distinction between proton and neutron bosons
Roth, John R.
2010-01-01
Experiments of Pränting and Andersson demonstrate how bacteria adapt to the growth limitation caused by antibiotic resistance mutations. The process of adaptation relies on gene copy number changes that arise at high rates, including duplications (10−4/cell/generation), amplifications (10−2/cell/generation) and mutant copy loss (10−2/cell/division). Reversible increases in copy number improve growth by small steps and provide more targets for rare sequence alterations (10−9/cell/division) tha...
Multipole stack for the 4 rings of the PS Booster
1976-01-01
The PS Booster (originally 800 MeV, now 1.4 GeV) saw first beam in 1972, routine operation began in 1973. The strive for ever higher intensities required the addition of multipoles. Manufacture of 8 stacks of multipoles was launched in 1974, for installation in 1976. For details, see 7511120X.
Fast heating and cooling in nanoimprint using a spring-loaded adapter in a preheated press
Schift, Helmut; Bellini, Sandro; Gobrecht, Jens;
2007-01-01
By using a spring-loaded adapter, instant heating and cooling of wafer-type substrates was implemented in standard hot embossing equipment. This was possible by using the well-known concept of a clamped stack of stamp and substrate, pre-assembled in an alignment fixture. A number of thermoplastic...
Fast automated protein NMR data collection and assignment by ADAPT-NMR on Bruker spectrometers
Lee, Woonghee; Hu, Kaifeng; Tonelli, Marco; Bahrami, Arash; Neuhardt, Elizabeth; Glass, Karen C.; Markley, John L.
2013-11-01
ADAPT-NMR (Assignment-directed Data collection Algorithm utilizing a Probabilistic Toolkit in NMR) supports automated NMR data collection and backbone and side chain assignment for [U-13C, U-15N]-labeled proteins. Given the sequence of the protein and data for the orthogonal 2D 1H-15N and 1H-13C planes, the algorithm automatically directs the collection of tilted plane data from a variety of triple-resonance experiments so as to follow an efficient pathway toward the probabilistic assignment of 1H, 13C, and 15N signals to specific atoms in the covalent structure of the protein. Data collection and assignment calculations continue until the addition of new data no longer improves the assignment score. ADAPT-NMR was first implemented on Varian (Agilent) spectrometers [A. Bahrami, M. Tonelli, S.C. Sahu, K.K. Singarapu, H.R. Eghbalnia, J.L. Markley, PLoS One 7 (2012) e33173]. Because of broader interest in the approach, we present here a version of ADAPT-NMR for Bruker spectrometers. We have developed two AU console programs (ADAPT_ORTHO_run and ADAPT_NMR_run) that run under TOPSPIN Versions 3.0 and higher. To illustrate the performance of the algorithm on a Bruker spectrometer, we tested one protein, chlorella ubiquitin (76 amino acid residues), that had been used with the Varian version: the Bruker and Varian versions achieved the same level of assignment completeness (98% in 20 h). As a more rigorous evaluation of the Bruker version, we tested a larger protein, BRPF1 bromodomain (114 amino acid residues), which yielded an automated assignment completeness of 86% in 55 h. Both experiments were carried out on a 500 MHz Bruker AVANCE III spectrometer equipped with a z-gradient 5 mm TCI probe. ADAPT-NMR is available at http://pine.nmrfam.wisc.edu/ADAPT-NMR in the form of pulse programs, the two AU programs, and instructions for installation and use.
Antonella Agodi; Martina Barchitta; Margherita Ferrante; Salvatore Sciacca; Ludwig Niessen
2005-01-01
Background. Food contamination by trichothecene mycotoxins is considered to be an emerging public health problem. The aim of this study was to validate a rapid sonification protocol, previously set up for cereal Fusarium DNA extraction from fast food samples, produced by a centre for research and development in the food industry in Catania, Sicily, and to validate it for a diagnostic PCR assay targeted at tri5, the key gene of trichothecene biosynthesis.&...
Generation of fast neturon spectra using an adaptive Gauss-Kronrod Quadrature algorithm
Triplett, Brian Scott
A lattice physics calculation is often the first step in analyzing a nuclear reactor. This calculation condenses regions of the reactor into average parameters (i.e., group constants) that can be used in coarser full-core, time-dependent calculations. This work presents a high-fidelity deterministic method for calculating the neutron energy spectrum in an infinite medium. The spectrum resulting from this calculation can be used to generate accurate group constants. This method includes a numerical algorithm based on Gauss-Kronrod Quadrature to determine the neutron transfer source to a given energy while controlling numerical error. This algorithm was implemented in a pointwise transport solver program called Pointwise Fast Spectrum Generator (PWFSG). PWFSG was benchmarked against the Monte Carlo program MCNP and another pointwise spectrum generation program, CENTRM, for a set of fast reactor infinite medium example cases. PWFSG showed good agreement with MCNP, yielding coefficients of determination above 98% for all example cases. In addition, PWFSG had 6 to 8 times lower flux estimation error than CENTRM in the cases examined. With run-times comparable to CENTRM, PWFSG represents a robust set of methods for generation of fast neutron spectra with increased accuracy without increased computational cost.
Microscopic Theory of Multipole Ordering in f-Electron Systems
Takashi Hotta
2012-01-01
Full Text Available A microscopic framework to determine multipole ordering in f-electron systems is provided on the basis of the standard quantum field theory. For the construction of the framework, a seven-orbital Hubbard Hamiltonian with strong spin-orbit coupling is adopted as a prototype model. A type of multipole and ordering vector is determined from the divergence of multipole susceptibility, which is evaluated in a random phase approximation. As an example of the application of the present framework, a multipole phase diagram on a three-dimensional simple cubic lattice is discussed for the case of n=2, where n denotes the average f-electron number per site. Finally, future problems concerning multipole ordering and fluctuations are briefly discussed.
Multipole Stack for the 800 MeV PS Booster
1975-01-01
The 800 MeV PS Booster had seen first beam in its 4 superposed rings in 1972, routine operation began in 1973. In the strive for ever higher beam intensities, the need for additional multipole lenses became evident. After detailed studies, the manufacture of 8 stacks of multipoles was launched in 1974. Each stack consists of 4 superposed multipoles and each multipole has 4 concentric shells. From the innermost to the outermost shell, Type A contains octupole, skew-octupole, sextupole, skew-sextupole. Type B contains skew-octupole, skew-sextupole, vertical dipole, horizontal dipole. Completion of installation in 1976 opened the way to higher beam intensities. M. Battiaz is seen here with a multipole stack and its many electrical connections.
Revising the Multipole Moments of Numerical Spacetimes and its Consequences
Pappas, George; Apostolatos, Theocharis A.
2012-06-01
Identifying the relativistic multipole moments of a spacetime of an astrophysical object that has been constructed numerically is of major interest, both because the multipole moments are intimately related to the internal structure of the object, and because the construction of a suitable analytic metric that mimics a numerical metric should be based on the multipole moments of the latter one in order to yield a reliable representation. In this Letter, we show that there has been a widespread delusion in the way the multipole moments of a numerical metric are read from the asymptotic expansion of the metric functions. We show how one should read correctly the first few multipole moments (starting from the quadrupole mass moment) and how these corrected moments improve the efficiency of describing the metric functions with analytic metrics that have already been used in the literature, as well as other consequences of using the correct moments.
Revising the multipole moments of numerical spacetimes, and its consequences
Pappas, George
2012-01-01
Identifying the relativistic multipole moments of a spacetime of an astrophysical object that has been constructed numerically is of major interest, both because the multipole moments are intimately related to the internal structure of the object, and because the construction of a suitable analytic metric that mimics a numerical metric should be based on the multipole moments of the latter one, in order to yield a reliable representation. In this note we show that there has been a widespread delusion in the way the multipole moments of a numerical metric are read from the asymptotic expansion of the metric functions. We show how one should read correctly the first few multipole moments (starting from the quadrupole mass-moment), and how these corrected moments improve the efficiency of describing the metric functions with analytic metrics that have already been used in the literature, as well as other consequences of using the correct moments.
Computation of electrostatic lenses and multipoles by the first order finite element method
The use of the first order finite element method for the computation of electrostatic rotationally symmetric electron lenses and electrostatic multipole elements is described. Attention is paid to the correct evaluation of the coefficients of FEM equations. Specific features in the program implementation include the use of thin electrodes, the use of dense fine mesh with smoothly varying step size, and the solution of the resulting set of linear equations with a fast and efficient ICCG method. The program is equipped with a suitable user interface allowing data input and modification graphically or by editing in clearly arranged menus, automeshing procedure for the generation of the fine mesh, and allowing display and plotter output of computed results. In a related program for electrostatic multipoles higher order harmonic components can be calculated. The use and the accuracy of the programs is illustrated by an example of the computation of a two-tube electrostatic lens, and by the evaluation of the axial mult ipole function for higher multipoles. (orig.)
Research on fast real-time adaptive audio mixing in multimedia conference
FAN Xing; GU Wei-kang; YE Xiu-qing
2005-01-01
In multimedia conference, the capability of audio processing is basic and requires more for real-time criteria. In this article, we categorize and analyze the schemes, and provide several multipoint speech audio mixing schemes using weighted algorithm, which meet the demand of practical needs for real-time multipoint speech mixing, for which the ASW and AEW schemes are especially recommended. Applying the adaptive algorithms, the high-performance schemes we provide do not use the saturation operation widely used in multimedia processing. Therefore, no additional noise will be added to the output. The above adaptive algorithms have relatively low computational complexity and good hearing perceptibility. The schemes are designed for parallel processing, and can be easily implemented with hardware, such as DSPs, and widely applied in multimedia conference systems.
Adaptive prediction in digitally controlled buck converter with fast load transient response
Lee, TLA; Chan, PCH
2012-01-01
An adaptive prediction scheme based on linear extrapolation for digitally controlled voltage-mode buck-type switching converter is presented. A major drawback of conventional digitally controlled switching converters is bandwidth limitation due to the additional phase lag in the digital feedback control loop. By predicting the future error voltage, the ADC sampling time delay is compensated in order to achieve a higher bandwidth even with a modest sampling rate. Both simulation and measuremen...
Jawerth, Bjoern; Sweldens, Wim
1993-01-01
We present ideas on how to use wavelets in the solution of boundary value ordinary differential equations. Rather than using classical wavelets, we adapt their construction so that they become (bi)orthogonal with respect to the inner product defined by the operator. The stiffness matrix in a Galerkin method then becomes diagonal and can thus be trivially inverted. We show how one can construct an O(N) algorithm for various constant and variable coefficient operators.
Amplifying effect of a release mechanism for fast adaptation in the hair bundle.
Sul, Bora; Iwasa, Kuni H
2009-07-01
A "release" mechanism, which has been experimentally observed as the fast component in the hair bundle's response to mechanical stimulation, appears similar to common mechanical relaxation with a damping effect. This observation is puzzling because such a response is expected to have an amplifying role in the mechanoelectrical transduction process in hair cells. Here it is shown that a release mechanism can indeed have a role in amplification, if it is associated with negative stiffness due to the gating of the mechonoelectric transducer channel. PMID:19603855
Investigation of the self field resonance instability of Ioffe-multipole stabilized electron rings
Experiments carried out on the RECE-Berta field reversing electron-ring device have identified a new type of anomalous particle loss caused by externally applied quadrupole-Ioffe magnetic field perturbations. Magnetic probe measurements of the ring magnetic fields and field reversal factor delta (= delta B/sub ring on axis/B/sub 0/) reveal three characteristic anomalous features in the ring decay when quadrupole fields are applied: (1) critical field reversal levels of delta/sub A/ < delta/sub B/ separating; (2) regimes of steady state anomalous decay between and above these levels, and (3) sudden dumps of fast ring electrons immediately before delta/sub A/ and delta/sub B/. Normal classical ring decay caused by collisions of fast ring electrons with neutral background gas is observed when delta < delta/sub A/. As part of an extensive series of experiments, the scaling of the three anomalous features with gas pressure and quadrupole field strength is studied. Additional measurements of x-ray emissions from a radial tungsten limiter target show that the dump of delta/sub A/ is coincident with a highly localized dump of fast electrons onto the limiter. Rings subjected to hexapole field perturbations show qualitatively similar decay anomalies. To obtain an explanation of these losses, the resonant coupling of particle orbital degrees of freedom by the externally applied symmetry breaking Ioffe-multipole field perturbations is investigated theoretically. The conditions for orbital instability of particles moving in the magnetic field of the electron ring and external mirror-multipole coils are calculated analytically, generalizing the results of Cohen et al. to all multipole orders. Unstable orbits occur when particle orbital frequencies w/sub z/ w/sub r/, and w/sub theta/, obey a resonance condition Nw/sub z/ + Lw/sub r/ + Mw/sub theata/ = 0 for theoretically calculated harmonic numbers (N,L,M)
Life in the Fast Lane: The Evolution of an Adaptive Vehicle Control System
Jochem, Todd; Pomerleau, Dean
1996-01-01
Giving robots the ability to operate in the real world has been, and continues to be, one of the most difficult tasks in AI research. Since 1987, researchers at Carnegie Mellon University have been investigating one such task. Their research has been focused on using adaptive, vision-based systems to increase the driving performance of the Navlab line of on-road mobile robots. This research has led to the development of a neural network system that can learn to drive on many road types simply...
Yang Zhiwei
2010-01-01
Full Text Available We propose a subspace-tracking-based space-time adaptive processing technique for airborne radar applications. By applying a modified approximated power iteration subspace tracing algorithm, the principal subspace in which the clutter-plus-interference reside is estimated. Therefore, the moving targets are detected by projecting the data on the minor subspace which is orthogonal to the principal subspace. The proposed approach overcomes the shortcomings of the existing methods and has satisfactory performance. Simulation results confirm that the performance improvement is achieved at very small secondary sample support, a feature that is particularly attractive for applications in heterogeneous environments.
Fast Fourier and Wavelet Transforms for Wavefront Reconstruction in Adaptive Optics
Dowla, F U; Brase, J M; Olivier, S S
2000-07-28
Wavefront reconstruction techniques using the least-squares estimators are computationally quite expensive. We compare wavelet and Fourier transforms techniques in addressing the computation issues of wavefront reconstruction in adaptive optics. It is shown that because the Fourier approach is not simply a numerical approximation technique unlike the wavelet method, the Fourier approach might have advantages in terms of numerical accuracy. However, strictly from a numerical computations viewpoint, the wavelet approximation method might have advantage in terms of speed. To optimize the wavelet method, a statistical study might be necessary to use the best basis functions or ''approximation tree.''
Collaborative Software Development in Support of Fast Adaptive AeroSpace Tools (FAAST)
Kleb, William L.; Nielsen, Eric J.; Gnoffo, Peter A.; Park, Michael A.; Wood, William A.
2003-01-01
A collaborative software development approach is described. The software product is an adaptation of proven computational capabilities combined with new capabilities to form the Agency's next generation aerothermodynamic and aerodynamic analysis and design tools. To efficiently produce a cohesive, robust, and extensible software suite, the approach uses agile software development techniques; specifically, project retrospectives, the Scrum status meeting format, and a subset of Extreme Programming's coding practices are employed. Examples are provided which demonstrate the substantial benefits derived from employing these practices. Also included is a discussion of issues encountered when porting legacy Fortran 77 code to Fortran 95 and a Fortran 95 coding standard.
D. Madurasinghe
2009-01-01
Full Text Available The proposed technique allows the radar receiver to accurately estimate the range of a large number of targets using a transmitter of opportunity as long as the location of the transmitter is known. The technique does not depend on the use of communication satellites or GPS systems, instead it relies on the availability of the direct transmit copy of the signal from the transmitter and the reflected paths off the various targets. An array-based space-fast time adaptive processor is implemented in order to estimate the path difference between the direct signal and the delayed signal, which bounces off the target. This procedure allows us to estimate the target distance as well as bearing.
GPU-based ultra-fast direct aperture optimization for online adaptive radiation therapy
Men, Chunhua; Jiang, Steve B
2010-01-01
Online adaptive radiation therapy (ART) has great promise to significantly reduce normal tissue toxicity and/or improve tumor control through real-time treatment adaptations based on the current patient anatomy. However, the major technical obstacle for clinical realization of online ART, namely the inability to achieve real-time efficiency in treatment re-planning, has yet to be solved. To overcome this challenge, this paper presents our work on the implementation of an intensity modulated radiation therapy (IMRT) direct aperture optimization (DAO) algorithm on graphics processing unit (GPU) based on our previous work on CPU. We formulate the DAO problem as a large-scale convex programming problem, and use an exact method called column generation approach to deal with its extremely large dimensionality on GPU. Five 9-field prostate and five 5-field head-and-neck IMRT clinical cases with 5\\times5 mm2 beamlet size and 2.5\\times2.5\\times2.5 mm3 voxel size were used to evaluate our algorithm on GPU. It takes onl...
SIRT1 Gain of Function Does Not Mimic or Enhance the Adaptations to Intermittent Fasting
Marie Boutant
2016-03-01
Full Text Available Caloric restriction (CR has been shown to prevent the onset of insulin resistance and to delay age-related physiological decline in mammalian organisms. SIRT1, a NAD+-dependent deacetylase enzyme, has been suggested to mediate the adaptive responses to CR, leading to the speculation that SIRT1 activation could be therapeutically used as a CR-mimetic strategy. Here, we used a mouse model of moderate SIRT1 overexpression to test whether SIRT1 gain of function could mimic or boost the metabolic benefits induced by every-other-day feeding (EODF. Our results indicate that SIRT1 transgenesis does not affect the ability of EODF to decrease adiposity and improve insulin sensitivity. Transcriptomic analyses revealed that SIRT1 transgenesis and EODF promote very distinct adaptations in individual tissues, some of which can be even be metabolically opposite, as in brown adipose tissue. Therefore, whereas SIRT1 overexpression and CR both improve glucose metabolism and insulin sensitivity, the etiologies of these benefits are largely different.
SIRT1 Gain of Function Does Not Mimic or Enhance the Adaptations to Intermittent Fasting.
Boutant, Marie; Kulkarni, Sameer S; Joffraud, Magali; Raymond, Frédéric; Métairon, Sylviane; Descombes, Patrick; Cantó, Carles
2016-03-01
Caloric restriction (CR) has been shown to prevent the onset of insulin resistance and to delay age-related physiological decline in mammalian organisms. SIRT1, a NAD(+)-dependent deacetylase enzyme, has been suggested to mediate the adaptive responses to CR, leading to the speculation that SIRT1 activation could be therapeutically used as a CR-mimetic strategy. Here, we used a mouse model of moderate SIRT1 overexpression to test whether SIRT1 gain of function could mimic or boost the metabolic benefits induced by every-other-day feeding (EODF). Our results indicate that SIRT1 transgenesis does not affect the ability of EODF to decrease adiposity and improve insulin sensitivity. Transcriptomic analyses revealed that SIRT1 transgenesis and EODF promote very distinct adaptations in individual tissues, some of which can be even be metabolically opposite, as in brown adipose tissue. Therefore, whereas SIRT1 overexpression and CR both improve glucose metabolism and insulin sensitivity, the etiologies of these benefits are largely different. PMID:26923584
Structural phase transitions in multipole traps
Marciante, Mathieu; Calisti, Annette; Knoop, Martina
2012-01-01
A small number of laser-cooled ions trapped in a linear radiofrequency multipole trap forms a hollow tube structure. We have studied, by means of molecular dynamics simulations, the structural transition from a double ring to a single ring of ions. We show that the single-ring configuration has the advantage to inhibit the thermal transfer from the rf-excited radial components of the motion to the axial component, allowing to reach the Doppler limit temperature along the direction of the trap axis. Once cooled in this particular configuration, the ions experience an angular dependency of the confinement if the local adiabaticity parameter exceeds the empirical limit. Bunching of the ion structures can then be observed and an analytic expression is proposed to take into account for this behaviour.
Electric multipole plasmons in deformed sodium clusters
Kleinig, W; Reinhard, P G
2001-01-01
The random-phase-approximation (RPA) method with separable residual forces (SRPA) is proposed for the description of multipole electric oscillations of valence electrons in deformed alkali metal clusters. Both the deformed mean field and residual interaction are derived self-consistently from the Kohn-Sham functional. SRPA drastically simplifies the computational effort which is urgent if not decisive for deformed systems. The method is applied to the description of dipole, quadrupole and octupole plasmons in deformed sodium clusters of a moderate size. We demonstrate that, in clusters with the size N>50, Landau damping successfully competes with deformation splitting and even becomes decisive in forming the width and gross structure of the dipole plasmon. Besides, the plasmon is generated by excitations from both ground state and shape isomers. In such clusters familiar experimental estimates for deformation splitting of dipole plasmon are useless.
Characterization of Multipole Microwave Cavity Filters
Reinhold, Philip; Schuster, David; Bishop, Lev S.
2013-03-01
An essential requirement for a quantum information processor is the ability to controllably couple and decouple individual qubits with each other. With superconducting circuit QED, this can be implemented by coupling multiple qubits to a transmission line cavity bus, and can be controlled by moving the qubit frequencies in and out of resonance with the bus. As coherence times increase, and the number of qubits attached to a bus grows larger, the problem of spurious coupling while detuned will become more important. We propose using principles from microwave filter design to create new couplers with higher contrast ratios in the effective qubit-qubit coupling. We present progress towards circuit implementations of multipole qubit coupling architectures.
Multipole strength distribution in 112Cd
Low-lying excited states of 112Cd have been investigated by means of the inelastic scattering of 30.7 and 60.0 MeV protons with an overall energy resolution of approximately 3x10-4. Most of the known levels, and several others observed for the first time, appear well resolved in final spectra up to about 4 MeV excitation energy. Spin assignments are proposed by comparing measured angular distributions with coupled-channels calculations. A large fragmentation of the hexadecapole strength is found. This is in contrast with the behaviour of the quadrupole and octupole strengths which are mostly concentrated in the first 2+ and 3- levels, respectively. The deduced multipole transition strengths are compared with the predictions of the interacting boson model in the sd, sdf, and sdg boson schemes. (orig.)
Stability of multipole-mode solitons in thermal nonlinear media
We study the stability of multipole-mode solitons in one-dimensional thermal nonlinear media. We show how the sample geometry impacts the stability of multipole-mode solitons and reveals that the tripole and quadrupole can be made stable in their whole domain of existence, provided that the sample width exceeds a critical value. In spite of such geometry-dependent soliton stability, we find that the maximal number of peaks in stable multipole-mode solitons in thermal media is the same as that in nonlinear materials with finite-range nonlocality.
Fast and intuitive programming of adaptive laser cutting of lace enabled by machine vision
Vaamonde, Iago; Souto-López, Álvaro; García-Díaz, Antón
2015-07-01
A machine vision system has been developed, validated, and integrated in a commercial laser robot cell. It permits an offline graphical programming of laser cutting of lace. The user interface allows loading CAD designs and aligning them with images of lace pieces. Different thread widths are discriminated to generate proper cutting program templates. During online operation, the system aligns CAD models of pieces and lace images, pre-checks quality of lace cuts and adapts laser parameters to thread widths. For pieces detected with the required quality, the program template is adjusted by transforming the coordinates of every trajectory point. A low-cost lace feeding system was also developed for demonstration of full process automation.
Xue, Xiaoming; Zhou, Jianzhong; Xu, Yanhe; Zhu, Wenlong; Li, Chaoshun
2015-10-01
Ensemble empirical mode decomposition (EEMD) represents a significant improvement over the original empirical mode decomposition (EMD) method for eliminating the mode mixing problem. However, the added white noises generate some tough problems including the high computational cost, the determination of the two critical parameters (the amplitude of the added white noise and the number of ensemble trials), and the contamination of the residue noise in the signal reconstruction. To solve these problems, an adaptively fast EEMD (AFEEMD) method combined with complementary EEMD (CEEMD) is proposed in this paper. In the proposed method, the two critical parameters are respectively fixed as 0.01 times standard deviation of the original signal and two ensemble trials. Instead, the upper frequency limit of the added white noise is the key parameter which needs to be prescribed beforehand. Unlike the original EEMD method, only two high-frequency white noises are added to the signal to be investigated with anti-phase in AFEEMD. Furthermore, an index termed relative root-mean-square error is employed for the adaptive selection of the proper upper frequency limit of the added white noises. Simulation test and vibration signals based fault diagnosis of rolling element bearing under different fault types are utilized to demonstrate the feasibility and effectiveness of the proposed method. The analysis results indicate that the AFEEMD method represents a sound improvement over the original EEMD method, and has strong practicability.
Stability of point-vortex multipoles revisited
Kizner, Ziv
2011-06-01
The point-vortex tripoles and pentapoles with zero total circulation are considered in the rigid-lid barotropic, equivalent-barotropic, and quasigeostrophic two-layer models. A tripole is assembled by three symmetrically arranged collinear vortices, while a pentapole by five vortices, of which one is located at the center of a square and four in the vertices of the square. The vortices on the sides, termed satellite vortices, are equal in strength and opposite in sign to the central vortex. To fulfill the zero-total-circulation condition, the central vortex is taken to be twice as strong as each of the satellite vortices in a tripole and four times as strong in a pentapole. In the two-layer model, two cases are distinguished, namely, the flat multipoles whose vortices are all located in the same layer and the carousel multipoles whose central vortex and satellite vortices reside in different layers. In all the models, the tripoles are shown to be nonlinearly stable and pentapoles, generally, unstable. Carousel pentapoles comparable in their size with the Rossby radius, and smaller, are exceptional in that they are stable to centrally symmetric perturbations (and, presumably, to arbitrary perturbations). The simple proof of the tripole stability is based on the fact that among the possible three-vortex configurations with zero total circulation characterized by the same (fixed) value of the Hamiltonian, there exists only one tripole, and, within the iso-Hamiltonian sheet, the squared linear momentum vanishes at this unique tripole only. This approach, being in essence universal for all models, works only with tripoles. For instance, a quadrupole cannot be treated in such a way, because there is a continuum of configurations of four vortices with zero total circulation on which the squared impulse vanishes. Dealing with pentapoles, we consider the perturbations that do not violate the central symmetry of the vortex configuration, fix the angular momentum, and examine
A fast, accurate and stable optimization algorithm is very important for inverse planning of intensity-modulated radiation therapy (IMRT), and for implementing dose-adaptive radiotherapy in the future. Conventional numerical search algorithms with positive beam weight constraints generally require numerous iterations and may produce suboptimal dose results due to trapping in local minima regions of the objective function landscape. A direct solution of the inverse problem using conventional quadratic objective functions without positive beam constraints is more efficient but it will result in unrealistic negative beam weights. We review here a direct solution of the inverse problem that is efficient and does not yield unphysical negative beam weights. In fast inverse dose optimization (FIDO) method the objective function for the optimization of a large number of beamlets is reformulated such that the optimization problem is reducible to a linear set of equations. The optimal set of intensities is then found through a matrix inversion, and negative beamlet intensities are avoided without the need for externally imposed ad hoc conditions. In its original version [S. P. Goldman, J. Z. Chen, and J. J. Battista, in Proceedings of the XIVth International Conference on the Use of Computers in Radiation Therapy, 2004, pp. 112-115; S. P. Goldman, J. Z. Chen, and J. J. Battista, Med. Phys. 32, 3007 (2005)], FIDO was tested on single two-dimensional computed tomography (CT) slices with sharp KERMA beams without scatter, in order to establish a proof of concept which demonstrated that FIDO could be a viable method for the optimization of cancer treatment plans. In this paper we introduce the latest advancements in FIDO that now include not only its application to three-dimensional volumes irradiated by beams with full scatter but include as well a complete implementation of clinical dose-volume constraints including maximum and minimum dose as well as equivalent uniform dose
Polarizable Atomic Multipole-based Molecular Mechanics for Organic Molecules
Ren, Pengyu; Wu, Chuanjie; Ponder, Jay W.
2011-01-01
An empirical potential based on permanent atomic multipoles and atomic induced dipoles is reported for alkanes, alcohols, amines, sulfides, aldehydes, carboxylic acids, amides, aromatics and other small organic molecules. Permanent atomic multipole moments through quadrupole moments have been derived from gas phase ab initio molecular orbital calculations. The van der Waals parameters are obtained by fitting to gas phase homodimer QM energies and structures, as well as experimental densities ...
Spherical spaces for cosmic topology and multipole selection rules
Kramer, Peter
2010-01-01
Spherical manifolds yield cosmic spaces with positive curvature. They result by closing pieces from the sphere used by Einstein for his initial cosmology. Harmonic analysis on the manifolds aims at explaining the observed low amplitudes at small multipole orders of the cosmic microwave background. We analyze assumptions of point symmetry and randomness for spherical spaces. There emerge four spaces named orbifolds, with low volume fraction from the sphere and sharp multipole selection rules i...
Searching the Force Field Electrostatic Multipole Parameter Space.
Jakobsen, Sofie; Jensen, Frank
2016-04-12
We show by tensor decomposition analyses that the molecular electrostatic potential for amino acid peptide models has an effective rank less than twice the number of atoms. This rank indicates the number of parameters that can be derived from the electrostatic potential in a statistically significant way. Using this as a guideline, we investigate different strategies for deriving a reduced set of atomic charges, dipoles, and quadrupoles capable of reproducing the reference electrostatic potential with a low error. A full combinatorial search of selected parameter subspaces for N-methylacetamide and a cysteine peptide model indicates that there are many different parameter sets capable of providing errors close to that of the global minimum. Among the different reduced multipole parameter sets that have low errors, there is consensus that atoms involved in π-bonding require higher order multipole moments. The possible correlation between multipole parameters is investigated by exhaustive searches of combinations of up to four parameters distributed in all possible ways on all possible atomic sites. These analyses show that there is no advantage in considering combinations of multipoles compared to a simple approach where the importance of each multipole moment is evaluated sequentially. When combined with possible weighting factors related to the computational efficiency of each type of multipole moment, this may provide a systematic strategy for determining a computational efficient representation of the electrostatic component in force field calculations. PMID:26925529
Park, Susanna B; Davare, Marco; Falla, Marika; Kennedy, William R; Selim, Mona M; Wendelschafer-Crabb, Gwen; Koltzenburg, Martin
2016-06-01
Sensory feedback from cutaneous mechanoreceptors in the fingertips is important in effective object manipulation, allowing appropriate scaling of grip and load forces during precision grip. However, the role of mechanoreceptor subtypes in these tasks remains incompletely understood. To address this issue, psychophysical tasks that may specifically assess function of type I fast-adapting (FAI) and slowly adapting (SAI) mechanoreceptors were used with object manipulation experiments to examine the regulation of grip force control in an experimental model of graded reduction in tactile sensitivity (healthy volunteers wearing 2 layers of latex gloves). With gloves, tactile sensitivity decreased significantly from 1.9 ± 0.4 to 12.3 ± 2.2 μm in the Bumps task assessing function of FAI afferents but not in a grating orientation task assessing SAI afferents (1.6 ± 0.1 to 1.8 ± 0.2 mm). Six axis force/torque sensors measured peak grip (PGF) and load (PLF) forces generated by the fingertips during a grip-lift task. With gloves there was a significant increase of PGF (14 ± 6%), PLF (17 ± 5%), and grip and load force rates (26 ± 8%, 20 ± 8%). A variable-weight series task was used to examine sensorimotor memory. There was a 20% increase in PGF when the lift of a light object was preceded by a heavy relative to a light object. This relationship was not significantly altered when lifting with gloves, suggesting that the addition of gloves did not change sensorimotor memory effects. We conclude that FAI fibers may be important for the online force scaling but not for the buildup of a sensorimotor memory. PMID:27052582
Bargatze, L. F.
2015-12-01
Active Data Archive Product Tracking (ADAPT) is a collection of software routines that permits one to generate XML metadata files to describe and register data products in support of the NASA Heliophysics Virtual Observatory VxO effort. ADAPT is also a philosophy. The ADAPT concept is to use any and all available metadata associated with scientific data to produce XML metadata descriptions in a consistent, uniform, and organized fashion to provide blanket access to the full complement of data stored on a targeted data server. In this poster, we present an application of ADAPT to describe all of the data products that are stored by using the Common Data File (CDF) format served out by the CDAWEB and SPDF data servers hosted at the NASA Goddard Space Flight Center. These data servers are the primary repositories for NASA Heliophysics data. For this purpose, the ADAPT routines have been used to generate data resource descriptions by using an XML schema named Space Physics Archive, Search, and Extract (SPASE). SPASE is the designated standard for documenting Heliophysics data products, as adopted by the Heliophysics Data and Model Consortium. The set of SPASE XML resource descriptions produced by ADAPT includes high-level descriptions of numerical data products, display data products, or catalogs and also includes low-level "Granule" descriptions. A SPASE Granule is effectively a universal access metadata resource; a Granule associates an individual data file (e.g. a CDF file) with a "parent" high-level data resource description, assigns a resource identifier to the file, and lists the corresponding assess URL(s). The CDAWEB and SPDF file systems were queried to provide the input required by the ADAPT software to create an initial set of SPASE metadata resource descriptions. Then, the CDAWEB and SPDF data repositories were queried subsequently on a nightly basis and the CDF file lists were checked for any changes such as the occurrence of new, modified, or deleted
Multipole matrix elements for 208Pb
New measurements of inelastic proton scattering to low-lying states of 208Pb at 200 and 400 MeV are reported. Deformation lengths extracted from angular distributions for the 3- (2.614 MeV), 51- (3.198 MeV), 52- (3.209 MeV), 2+ (4.086 MeV) and 4+ (4.324 MeV) states are in good accord with values extracted at other incident proton energies. The fact that the deformation lengths are independent of incident proton energy within experimental uncertainty provides support for the validity of the collective DWBA for medium energy proton scattering to strongly excited states. Advantage is taken of this to extract statistically more precise values of the ratio of neutron to proton multipole matrix elements (M/sub n//M/sub p/). Different methods of determining the appropriate average value of M/sub n//M/sub p/ are discussed. 3 refs., 1 tab
Exchange splitting of the interaction energy and the multipole expansion of the wave function
The exchange splitting J of the interaction energy of the hydrogen atom with a proton is calculated using the conventional surface-integral formula Jsurf[Φ], the volume-integral formula of the symmetry-adapted perturbation theory JSAPT[Φ], and a variational volume-integral formula Jvar[Φ]. The calculations are based on the multipole expansion of the wave function Φ, which is divergent for any internuclear distance R. Nevertheless, the resulting approximations to the leading coefficient j0 in the large-R asymptotic series J(R) = 2e−R−1R(j0 + j1R−1 + j2R−2 + ⋯) converge with the rate corresponding to the convergence radii equal to 4, 2, and 1 when the Jvar[Φ], Jsurf[Φ], and JSAPT[Φ] formulas are used, respectively. Additionally, we observe that also the higher jk coefficients are predicted correctly when the multipole expansion is used in the Jvar[Φ] and Jsurf[Φ] formulas. The symmetry adapted perturbation theory formula JSAPT[Φ] predicts correctly only the first two coefficients, j0 and j1, gives a wrong value of j2, and diverges for higher jn. Since the variational volume-integral formula can be easily generalized to many-electron systems and evaluated with standard basis-set techniques of quantum chemistry, it provides an alternative for the determination of the exchange splitting and the exchange contribution of the interaction potential in general
Cardiac magnetic source imaging based on current multipole model
Tang Fa-Kuan; Wang Qian; Hua Ning; Lu Hong; Tang Xue-Zheng; Ma Ping
2011-01-01
It is widely accepted that the heart current source can be reduced into a current multipole. By adopting three linear inverse methods, the cardiac magnetic imaging is achieved in this article based on the current multipole model expanded to the first order terms. This magnetic imaging is realized in a reconstruction plane in the centre of human heart, where the current dipole array is employed to represent realistic cardiac current distribution. The current multipole as testing source generates magnetic fields in the measuring plane, serving as inputs of cardiac magnetic inverse problem. In the heart-torso model constructed by boundary element method, the current multipole magnetic field distribution is compared with that in the homogeneous infinite space, and also with the single current dipole magnetic field distribution.Then the minimum-norm least-squares (MNLS) method, the optimal weighted pseuDOInverse method (OWPIM), and the optimal constrained linear inverse method (OCLIM) are selected as the algorithms for inverse computation based on current multipole model innovatively, and the imaging effects of these three inverse methods are compared. Besides,two reconstructing parameters, residual and mean residual, are also discussed, and their trends under MNLS, OWPIM and OCLIM each as a function of SNR are obtained and compared.
HIFI-C: a robust and fast method for determining NMR couplings from adaptive 3D to 2D projections
We describe a novel method for the robust, rapid, and reliable determination of J couplings in multi-dimensional NMR coupling data, including small couplings from larger proteins. The method, 'High-resolution Iterative Frequency Identification of Couplings' (HIFI-C) is an extension of the adaptive and intelligent data collection approach introduced earlier in HIFI-NMR. HIFI-C collects one or more optimally tilted two-dimensional (2D) planes of a 3D experiment, identifies peaks, and determines couplings with high resolution and precision. The HIFI-C approach, demonstrated here for the 3D quantitative J method, offers vital features that advance the goal of rapid and robust collection of NMR coupling data. (1) Tilted plane residual dipolar couplings (RDC) data are collected adaptively in order to offer an intelligent trade off between data collection time and accuracy. (2) Data from independent planes can provide a statistical measure of reliability for each measured coupling. (3) Fast data collection enables measurements in cases where sample stability is a limiting factor (for example in the presence of an orienting medium required for residual dipolar coupling measurements). (4) For samples that are stable, or in experiments involving relatively stronger couplings, robust data collection enables more reliable determinations of couplings in shorter time, particularly for larger biomolecules. As a proof of principle, we have applied the HIFI-C approach to the 3D quantitative J experiment to determine N-C' RDC values for three proteins ranging from 56 to 159 residues (including a homodimer with 111 residues in each subunit). A number of factors influence the robustness and speed of data collection. These factors include the size of the protein, the experimental set up, and the coupling being measured, among others. To exhibit a lower bound on robustness and the potential for time saving, the measurement of dipolar couplings for the N-C' vector represents a realistic
Gupta, V; Wang, Y; Romero, A; Heijmen, B; Hoogeman, M [Erasmus MC Cancer Institute, Rotterdam (Netherlands); Myronenko, A; Jordan, P [Accuray Incorporated, Sunnyvale, United States. (United States)
2014-06-01
Purpose: Various studies have demonstrated that online adaptive radiotherapy by real-time re-optimization of the treatment plan can improve organs-at-risk (OARs) sparing in the abdominal region. Its clinical implementation, however, requires fast and accurate auto-segmentation of OARs in CT scans acquired just before each treatment fraction. Autosegmentation is particularly challenging in the abdominal region due to the frequently observed large deformations. We present a clinical validation of a new auto-segmentation method that uses fully automated non-rigid registration for propagating abdominal OAR contours from planning to daily treatment CT scans. Methods: OARs were manually contoured by an expert panel to obtain ground truth contours for repeat CT scans (3 per patient) of 10 patients. For the non-rigid alignment, we used a new non-rigid registration method that estimates the deformation field by optimizing local normalized correlation coefficient with smoothness regularization. This field was used to propagate planning contours to repeat CTs. To quantify the performance of the auto-segmentation, we compared the propagated and ground truth contours using two widely used metrics- Dice coefficient (Dc) and Hausdorff distance (Hd). The proposed method was benchmarked against translation and rigid alignment based auto-segmentation. Results: For all organs, the auto-segmentation performed better than the baseline (translation) with an average processing time of 15 s per fraction CT. The overall improvements ranged from 2% (heart) to 32% (pancreas) in Dc, and 27% (heart) to 62% (spinal cord) in Hd. For liver, kidneys, gall bladder, stomach, spinal cord and heart, Dc above 0.85 was achieved. Duodenum and pancreas were the most challenging organs with both showing relatively larger spreads and medians of 0.79 and 2.1 mm for Dc and Hd, respectively. Conclusion: Based on the achieved accuracy and computational time we conclude that the investigated auto
HIFI-C: a robust and fast method for determining NMR couplings from adaptive 3D to 2D projections
Cornilescu, Gabriel, E-mail: gabrielc@nmrfam.wisc.edu; Bahrami, Arash; Tonelli, Marco; Markley, John L.; Eghbalnia, Hamid R. [University of Wisconsin-Madison, National Magnetic Resonance Facility at Madison (United States)], E-mail: eghbalni@nmrfam.wisc.edu
2007-08-15
We describe a novel method for the robust, rapid, and reliable determination of J couplings in multi-dimensional NMR coupling data, including small couplings from larger proteins. The method, 'High-resolution Iterative Frequency Identification of Couplings' (HIFI-C) is an extension of the adaptive and intelligent data collection approach introduced earlier in HIFI-NMR. HIFI-C collects one or more optimally tilted two-dimensional (2D) planes of a 3D experiment, identifies peaks, and determines couplings with high resolution and precision. The HIFI-C approach, demonstrated here for the 3D quantitative J method, offers vital features that advance the goal of rapid and robust collection of NMR coupling data. (1) Tilted plane residual dipolar couplings (RDC) data are collected adaptively in order to offer an intelligent trade off between data collection time and accuracy. (2) Data from independent planes can provide a statistical measure of reliability for each measured coupling. (3) Fast data collection enables measurements in cases where sample stability is a limiting factor (for example in the presence of an orienting medium required for residual dipolar coupling measurements). (4) For samples that are stable, or in experiments involving relatively stronger couplings, robust data collection enables more reliable determinations of couplings in shorter time, particularly for larger biomolecules. As a proof of principle, we have applied the HIFI-C approach to the 3D quantitative J experiment to determine N-C' RDC values for three proteins ranging from 56 to 159 residues (including a homodimer with 111 residues in each subunit). A number of factors influence the robustness and speed of data collection. These factors include the size of the protein, the experimental set up, and the coupling being measured, among others. To exhibit a lower bound on robustness and the potential for time saving, the measurement of dipolar couplings for the N-C' vector
HIFI-C: a robust and fast method for determining NMR couplings from adaptive 3D to 2D projections.
Cornilescu, Gabriel; Bahrami, Arash; Tonelli, Marco; Markley, John L; Eghbalnia, Hamid R
2007-08-01
We describe a novel method for the robust, rapid, and reliable determination of J couplings in multi-dimensional NMR coupling data, including small couplings from larger proteins. The method, "High-resolution Iterative Frequency Identification of Couplings" (HIFI-C) is an extension of the adaptive and intelligent data collection approach introduced earlier in HIFI-NMR. HIFI-C collects one or more optimally tilted two-dimensional (2D) planes of a 3D experiment, identifies peaks, and determines couplings with high resolution and precision. The HIFI-C approach, demonstrated here for the 3D quantitative J method, offers vital features that advance the goal of rapid and robust collection of NMR coupling data. (1) Tilted plane residual dipolar couplings (RDC) data are collected adaptively in order to offer an intelligent trade off between data collection time and accuracy. (2) Data from independent planes can provide a statistical measure of reliability for each measured coupling. (3) Fast data collection enables measurements in cases where sample stability is a limiting factor (for example in the presence of an orienting medium required for residual dipolar coupling measurements). (4) For samples that are stable, or in experiments involving relatively stronger couplings, robust data collection enables more reliable determinations of couplings in shorter time, particularly for larger biomolecules. As a proof of principle, we have applied the HIFI-C approach to the 3D quantitative J experiment to determine N-C' RDC values for three proteins ranging from 56 to 159 residues (including a homodimer with 111 residues in each subunit). A number of factors influence the robustness and speed of data collection. These factors include the size of the protein, the experimental set up, and the coupling being measured, among others. To exhibit a lower bound on robustness and the potential for time saving, the measurement of dipolar couplings for the N-C' vector represents a realistic
Purpose: Various studies have demonstrated that online adaptive radiotherapy by real-time re-optimization of the treatment plan can improve organs-at-risk (OARs) sparing in the abdominal region. Its clinical implementation, however, requires fast and accurate auto-segmentation of OARs in CT scans acquired just before each treatment fraction. Autosegmentation is particularly challenging in the abdominal region due to the frequently observed large deformations. We present a clinical validation of a new auto-segmentation method that uses fully automated non-rigid registration for propagating abdominal OAR contours from planning to daily treatment CT scans. Methods: OARs were manually contoured by an expert panel to obtain ground truth contours for repeat CT scans (3 per patient) of 10 patients. For the non-rigid alignment, we used a new non-rigid registration method that estimates the deformation field by optimizing local normalized correlation coefficient with smoothness regularization. This field was used to propagate planning contours to repeat CTs. To quantify the performance of the auto-segmentation, we compared the propagated and ground truth contours using two widely used metrics- Dice coefficient (Dc) and Hausdorff distance (Hd). The proposed method was benchmarked against translation and rigid alignment based auto-segmentation. Results: For all organs, the auto-segmentation performed better than the baseline (translation) with an average processing time of 15 s per fraction CT. The overall improvements ranged from 2% (heart) to 32% (pancreas) in Dc, and 27% (heart) to 62% (spinal cord) in Hd. For liver, kidneys, gall bladder, stomach, spinal cord and heart, Dc above 0.85 was achieved. Duodenum and pancreas were the most challenging organs with both showing relatively larger spreads and medians of 0.79 and 2.1 mm for Dc and Hd, respectively. Conclusion: Based on the achieved accuracy and computational time we conclude that the investigated auto
Objective: There is substantial debate regarding the appropriate protocol for ventilatory management in fast-track cardiac anesthesia (FTCA). This study was carried out to assess and compare the risks and benefits of respiratory weaning based on adaptive support ventilation (ASV) and synchronized intermittent mandatory ventilation (SIMV) after uncomplicated cardiac surgery. Methodology: In a randomized clinical trial, after receiving approval of the Department Research Committee and informed consent from study subjects, 100 patients undergoing elective coronary artery bypass graft (CABG) surgery with cardiopulmonary bypass (CPB) were enrolled during a 4-month period at a university-based hospital. After surgery and admission to the intensive care unit (ICU), patients were randomized to ASV and SIMV groups. Arterial blood gas (ABG) and hemodynamic variables, respiratory and ventilator characteristics including lung compliance, rapid shallow breathing index (RSBI), tidal volume (TV), respiratory rate (RR), peak inspiratory pressure (P peak), mean airway pressure (p mean), Pao2/FIo2, duration of mechanical ventilation and tracheal intubation, and length of ICU stay were recorded and compared between the two groups. The data were analyzed in 82 patients after considering the exclusion criteria. Results: There were no differences between ASV and SIMV groups in demographics and preoperative characteristics. The duration of tracheal intubation and the length of ICU stay were similar in both groups. There were no statistically and clinically relevant differences between the two groups in ABG, hemodynamic changes, and respiratory and ventilator characteristics during ICU stay. Conclusion: Although ASV may facilitate postoperative respiratory management in FTCA, both ASV and SIMV provide similarly safe and practicable respiratory weaning in the cardiac ICU. The evaluation of potential advantages in patient outcomes and resource utilization of respiratory weaning based on ASV
Beam injection with pulsed multipole magnet at UVSOR-III
Yamamoto, N., E-mail: naoto@nagoya-u.jp [Synchrotron Radiation Research Center, Nagoya University, Furo-cho, Chikusa-ku, Nagoya, Aich 464-8603 (Japan); Zen, H. [Institute of Advanced Energy, Kyoto University, Gokasho, Uji, Kyoto 611-0011 (Japan); Hosaka, M. [Synchrotron Radiation Research Center, Nagoya University, Furo-cho, Chikusa-ku, Nagoya, Aich 464-8603 (Japan); Konomi, T. [UVSOR, Institute for Molecular Science, 38 Nishigo-Naka, Myodaiji, Okazaki 444-8585 (Japan); Adachi, M. [High Energy Accelerator Research Organization, KEK 1-1 Oho, Tsukuba, Ibaraki 305-0801 (Japan); Hayashi, K.; Yamazaki, J. [UVSOR, Institute for Molecular Science, 38 Nishigo-Naka, Myodaiji, Okazaki 444-8585 (Japan); Takashima, Y. [Synchrotron Radiation Research Center, Nagoya University, Furo-cho, Chikusa-ku, Nagoya, Aich 464-8603 (Japan); Katoh, M. [UVSOR, Institute for Molecular Science, 38 Nishigo-Naka, Myodaiji, Okazaki 444-8585 (Japan)
2014-12-11
In this study, we designed and manufactured a pulsed multipole magnet for beam injection into the UVSOR-III ring. A sextupole-like magnetic field could be excited when using the multipole magnet. To compensate for the residual field at the center of the magnet caused by manufacturing imprecisions, thin ferrite sheets were used. The injection experiments at UVSOR-III demonstrated multi-turn injections with the pulsed multipole magnet. The injection efficiency was 23% and the electron beam was stored up to the normal operation current of 300 mA. Moreover, we confirmed that oscillations of stored beams caused by beam injection were drastically suppressed compared with conventional pulsed dipole injection.
Multipole expansions in four-dimensional hyperspherical harmonics
The technique of vector differentiation is applied to the problem of the derivation of multipole expansions in four-dimensional space. Explicit expressions for the multipole expansion of the function rnCj(r-circumflex) with r = r1 + r2 are given in terms of tensor products of two hyperspherical harmonics depending on the unit vectors r-circumflex1 and r-circumflex2. The multipole decomposition of the function (r1 . r2)n is also derived. The proposed method can be easily generalized to the case of the space with dimensionality larger than four. Several explicit expressions for the four-dimensional Clebsch-Gordan coefficients with particular values of parameters are presented in the closed form
Reply to comment on 'On the origin dependence of multipole moments in electromagnetism'
De Visschere, Patrick
2010-01-01
In this reply, we emphasize that the goal of our original paper was to show that a straightforward application of standard multipole theory does not lead to unphysical results as claimed by Raab and de Lange. We stress once more that an origin shift for calculating the multipoles must be accompanied by a relocation of these multipoles, which compensates the effects of the origin dependence of the multipoles. We point out that the position of the macroscopic boundary is a relevant parameter. W...
Numerical simulation of Electron Energy Loss Spectroscopy using a Generalized Multipole Technique
Kiewidt, Lars [Foundation Institute of Materials Science (IWT), Department of Production Engineering, University of Bremen, Badgasteiner Str. 3, 28359 Bremen (Germany); Karamehmedović, Mirza, E-mail: mirza@iwt.uni-bremen.de [Foundation Institute of Materials Science (IWT), Department of Production Engineering, University of Bremen, Badgasteiner Str. 3, 28359 Bremen (Germany); Matyssek, Christian [Humboldt Universität zu Berlin, AG Theoretische Optik and Photonik, Newtonstr. 15, 12489 Berlin (Germany); Hergert, Wolfram [Martin Luther University Halle, Institute of Physics, Theory Group, von-Seckendorff-Platz 1, 06120 Halle (Germany); Mädler, Lutz [Foundation Institute of Materials Science (IWT), Department of Production Engineering, University of Bremen, Badgasteiner Str. 3, 28359 Bremen (Germany); Institut für Werkstofftechnik, Badgasteiner Str. 3, 28359 Bremen (Germany); Wriedt, Thomas [Institut für Werkstofftechnik, Badgasteiner Str. 3, 28359 Bremen (Germany)
2013-10-15
We numerically simulate low-loss Electron Energy Loss Spectroscopy (EELS) of isolated spheroidal nanoparticles, using an electromagnetic model based on a Generalized Multipole Technique (GMT). The GMT is fast and accurate, and, in principle, flexible regarding nanoparticle shape and the incident electron beam. The implemented method is validated against reference analytical and numerical methods for plane-wave scattering by spherical and spheroidal nanoparticles. Also, simulated electron energy loss (EEL) spectra of spherical and spheroidal nanoparticles are compared to available analytical and numerical solutions. An EEL spectrum is predicted numerically for a prolate spheroidal aluminum nanoparticle. The presented method is the basis for a powerful tool for the computation, analysis and interpretation of EEL spectra of general geometric configurations. - Highlights: • We simulate Electron Energy Loss Spectroscopy using a Generalized Multipole Technique. • We achieve good correspondence with reference methods for spherical nanoparticles. • The presented method is a basis for the computation and interpretation of Electron Energy Loss spectra of general geometric configurations.
Numerical simulation of Electron Energy Loss Spectroscopy using a Generalized Multipole Technique
We numerically simulate low-loss Electron Energy Loss Spectroscopy (EELS) of isolated spheroidal nanoparticles, using an electromagnetic model based on a Generalized Multipole Technique (GMT). The GMT is fast and accurate, and, in principle, flexible regarding nanoparticle shape and the incident electron beam. The implemented method is validated against reference analytical and numerical methods for plane-wave scattering by spherical and spheroidal nanoparticles. Also, simulated electron energy loss (EEL) spectra of spherical and spheroidal nanoparticles are compared to available analytical and numerical solutions. An EEL spectrum is predicted numerically for a prolate spheroidal aluminum nanoparticle. The presented method is the basis for a powerful tool for the computation, analysis and interpretation of EEL spectra of general geometric configurations. - Highlights: • We simulate Electron Energy Loss Spectroscopy using a Generalized Multipole Technique. • We achieve good correspondence with reference methods for spherical nanoparticles. • The presented method is a basis for the computation and interpretation of Electron Energy Loss spectra of general geometric configurations
A genetic algorithm for optimizing multi-pole Debye models of tissue dielectric properties
Models of tissue dielectric properties (permittivity and conductivity) enable the interactions of tissues and electromagnetic fields to be simulated, which has many useful applications in microwave imaging, radio propagation, and non-ionizing radiation dosimetry. Parametric formulae are available, based on a multi-pole model of tissue dispersions, but although they give the dielectric properties over a wide frequency range, they do not convert easily to the time domain. An alternative is the multi-pole Debye model which works well in both time and frequency domains. Genetic algorithms are an evolutionary approach to optimization, and we found that this technique was effective at finding the best values of the multi-Debye parameters. Our genetic algorithm optimized these parameters to fit to either a Cole–Cole model or to measured data, and worked well over wide or narrow frequency ranges. Over 10 Hz–10 GHz the best fits for muscle, fat or bone were each found for ten dispersions or poles in the multi-Debye model. The genetic algorithm is a fast and effective method of developing tissue models that compares favourably with alternatives such as the rational polynomial fit. (paper)
Stability of Multipole-mode Solitons in Thermal Nonlinear Media
Dong, Liangwei; Ye, Fangwei
2010-01-01
We study the stability of multipole-mode solitons in one-dimensional thermal nonlinear media. We show how the sample geometry impacts the stability of mutlipole-mode solitons and reveal that the tripole and quadrupole can be made stable in their whole domain of existence, provided that the sample width exceeds a critical value. In spite of such geometry-dependent soliton stability, we find that the maximal number of peaks in stable multipole-mode solitons in thermal media is the same as that ...
Multipole effects to the opacity of hot dense gold plasma
ZENG Lei; JIN Feng-tao; YUAN Jian-min
2006-01-01
The contributions of the multipole transitions to the opacity of hot dense gold plasma are taken into account by using an average-atom model.The influences of the E2,E3 and E4 transitions on the Rosseland opacity are studied,respectively.Comparisons with Miao's calculation have been made.It shows that using the Taylor series to account for the multipole transitions is no longer valid since ik.r is not much smaller than the unit when the photon energy goes very high.
Multipole matrix elements for Coulomb excitation of rotation aligned brand in 164Er
The possibility of exciting high spin states of the rotation-aligned band is studied. A microscopic model has been employed to estimate large l multipole matrix elements for transitions to high spin states. Possible multipole Coulomb excitation paths are given. Large l multipole transitions to high spin states by inelastic proton scattering are studied. (author). 19 refs., 4 figs
Exchange splitting of the interaction energy and the multipole expansion of the wave function
Gniewek, Piotr
2015-01-01
The exchange splitting $J$ of the interaction energy of the hydrogen atom with a proton is calculated using the conventional surface-integral formula $J_{\\textrm{surf}}[\\varphi]$, the volume-integral formula of the symmetry-adapted perturbation theory $J_{\\textrm{SAPT}}[\\varphi]$, and a variational volume-integral formula $J_{\\textrm{var}}[\\varphi]$. The calculations are based on the multipole expansion of the wave function $\\varphi$, which is divergent for any internuclear distance $R$. Nevertheless, the resulting approximations to the leading coefficient $j_0$ in the large-$R$ asymptotic series $J(R) = 2 e^{-R-1} R ( j_0 + j_1 R^{-1} + j_2 R^{-2} +\\cdots ) $ converge, with the rate corresponding to the convergence radii equal to 4, 2, and 1 when the $J_{\\textrm{var}}[\\varphi]$, $J_{\\textrm{surf}}[\\varphi]$, and $J_{\\textrm{SAPT}}[\\varphi]$ formulas are used, respectively. Additionally, we observe that also the higher $j_k$ coefficients are predicted correctly when the multipole expansion is used in the $J_{...
Mohr, Stephan; Genovese, Luigi; Ratcliff, Laura; Masella, Michel
The quantum mechanics/molecular mechanis (QM/MM) method is a popular approach that allows to perform atomistic simulations using different levels of accuracy. Since only the essential part of the simulation domain is treated using a highly precise (but also expensive) QM method, whereas the remaining parts are handled using a less accurate level of theory, this approach allows to considerably extend the total system size that can be simulated without a notable loss of accuracy. In order to couple the QM and MM regions we use an approximation of the electrostatic potential based on a multipole expansion. The multipoles of the QM region are determined based on the results of a linear scaling Density Functional Theory (DFT) calculation using a set of adaptive, localized basis functions, as implemented within the BigDFT software package. As this determination comes at virtually no extra cost compared to the QM calculation, the coupling between QM and MM region can be done very efficiently. In this presentation I will demonstrate the accuracy of both the linear scaling DFT approach itself as well as of the approximation of the electrostatic potential based on the multipole expansion, and show some first QM/MM applications using the aforementioned approach.
Kinnunen, Sanni; Mänttäri, Satu; Herzig, Karl-Heinz; Nieminen, Petteri; Mustonen, Anne-Mari; Saarela, Seppo
2016-02-01
The raccoon dog (Nyctereutes procyonoides) is a canid with autumnal fattening and passive wintering strategy. We examined the effects of wintertime fasting and seasonality on AMP-activated protein kinase (AMPK), a regulator of metabolism, and its target, acetyl-CoA carboxylase (ACC) on the species. Twelve farmed raccoon dogs (eleven females/one male) were divided into two groups: half were fasted for ten weeks in December-March (winter fasted) and the others were fed ad libitum (winter fed). A third group (autumn fed, eight females) was fed ad libitum and sampled in December. Total AMPK, ACC and their phosphorylated forms (pAMPK, pACC) were measured from hypothalamus, liver, intra-abdominal (iWAT) and subcutaneous white adipose tissues (sWAT). The fasted animals lost 32% and the fed 20% of their body mass. Hypothalamic AMPK expression was lower and pACC levels higher in the winter groups compared to the autumn fed group. Liver pAMPK was lower in the winter fasted group, with consistently decreased ACC and pACC. AMPK and pAMPK were down-regulated in sWAT and iWAT of both winter groups, with a parallel decline in pACC in sWAT. The responses of AMPK and ACC to fasting were dissimilar to the effects observed previously in non-seasonal mammals and hibernators. Differences between the winter fed and autumn fed groups indicate that the functions of AMPK and ACC could be regulated in a season-dependent manner. Furthermore, the distinctive effects of prolonged fasting and seasonal adaptation on AMPK-ACC pathway could contribute to the wintering strategy of the raccoon dog. PMID:26603554
Multipole moments for embedding potentials: Exploring different atomic allocation algorithms.
S Nørby, Morten; Magnus Haugaard Olsen, Jógvan; Kongsted, Jacob; Aagard Jensen, Hans Jørgen
2016-07-01
Polarizable quantum mechanical (QM)/molecular mechanics (MM)-embedding methods are currently among the most promising methods for computationally feasible, yet reliable, production calculations of localized excitations and molecular response properties of large molecular complexes, such as proteins and RNA/DNA, and of molecules in solution. Our aim is to develop a computational methodology for distributed multipole moments and their associated multipole polarizabilities which is accurate, computationally efficient, and with smooth convergence with respect to multipole order. As the first step toward this goal, we herein investigate different ways of obtaining distributed atom-centered multipole moments that are used in the construction of the electrostatic part of the embedding potential. Our objective is methods that not only are accurate and computationally efficient, but which can be consistently extended with site polarizabilities including internal charge transfer terms. We present a new way of dealing with well-known problems in relation to the use of basis sets with diffuse functions in conventional atomic allocation algorithms, avoiding numerical integration schemes. Using this approach, we show that the classical embedding potential can be systematically improved, also when using basis sets with diffuse functions, and that very accurate embedding potentials suitable for QM/MM embedding calculations can be acquired. © 2016 Wiley Periodicals, Inc. PMID:27187063
Nonrelativistic atom-photon interaction beyond the multipole approximation
Boers, D.J.; Wijers, C.M.J.
2003-01-01
We investigate the interaction between the hydrogen atomic orbitals and the quantized modes of the electromagnetic field within the domain of nonrelativistic quantum electrodynamics in the Coulomb gauge. Contrary to the conventional dipole approximation and higher-order multipole approximations, whi
Concept of multipole magnetic ﬁeld rotation in ECRIS
M H Rashid; R K Bhandari
2002-11-01
The conventional type of magnetic well is formed by superposition of two types of magnetic ﬁeld, axial bumpy ﬁeld and radial multipole ﬁeld. It is used to contain plasma that consists of neutrals, ions and electrons. These particles are in constant motion in the well and energetic electrons create plasma by violent collisions with neutrals and ions. The conﬁned electrons are constantly heated by ECR technique in the presence of magnetic ﬁeld. In this paper it has been shown theoretically that how the electron motion is inﬂuenced in terms of heating, containment and azimuthal uniformity of plasma, by the axial rotation of the multipole magnetic ﬁeld [1,2]. Afterwards, the feasibility of achieving a rotating magnetic multipole ﬁeld is discussed to some extent. And it is seen that it is not beyond the capability of the scientiﬁc community in the present scenario of the advanced technology. Presently, it can be achieved for lesser ﬁeld and slightly larger size of the multipole electromagnet and can be used for improvement of the ECR ion source (ECRIS).
Multipole analyses and photo-decay couplings at intermediate energies
The authors describe the results of several multipole analyses of pion-photoproduction data to 2 GeV in the lab photon energy. Comparisons are made with previous analyses. The photo-decay couplings for the delta are examined in detail. Problems in the representation of photoproduction data are discussed, with an emphasis on the recent LEGS data. 16 refs., 4 tabs
Two-dimensional multipole solitons in nonlocal nonlinear media.
Rotschild, Carmel; Segev, Mordechai; Xu, Zhiyong; Kartashov, Yaroslav V; Torner, Lluis; Cohen, Oren
2006-11-15
We present the experimental observation of scalar multipole solitons in highly nonlocal nonlinear media, including dipole, tripole, quadrupole, and necklace-type solitons, organized as arrays of out-of-phase bright spots. These complex solitons are metastable, but with a large parameters range where the instability is weak, permitting their experimental observation. PMID:17072407
Tune spread due to magnetic multipoles in RHIC
Analytical expressions have been obtained of the amplitude and momentum dependence of the transverse tunes due to magnetic multipoles and orbit misalignment. Based on these expressions, compensation methods are developed to minimize the tune spread in RHIC with the β* = lm design
Spiralling solitons and multipole localized modes in nonlocal nonlinear media
Buccoliero, Daniel; Lopez-Aguayo, Servando; Skupin, Stefan;
2007-01-01
We analyze the propagation of rotating multi-soliton localized structures in optical media with spatially nonlocal nonlinearity. We demonstrate that nonlocality stabilizes the azimuthal breakup of rotating dipole as well as multipole localized soliton modes. We compare the results for two differe...
Linear-scaling multipole-accelerated Gaussian and finite-element Coulomb method
Watson, Mark A.; Kurashige, Yuki; Nakajima, Takahito; Hirao, Kimihiko
2008-02-01
A linear-scaling implementation of the Gaussian and finite-element Coulomb (GFC) method is presented for the rapid computation of the electronic Coulomb potential. The current work utilizes the fast multipole method (FMM) for the evaluation of the Poisson equation boundary condition. The FMM affords significant savings for small- and medium-sized systems and overcomes the bottleneck in the GFC method for very large systems. Compared to an exact analytical treatment of the boundary, more than 100-fold speedups are observed for systems with more than 1000 basis functions without any significant loss of accuracy. We present CPU times to demonstrate the effectiveness of the linear-scaling GFC method for both one-dimensional polyalanine chains and the challenging case of three-dimensional diamond fragments.
Chow, Brian Y.; Chuong, Amy S.; Klapoetke, Nathan C; Boyden, Edward S.
2011-01-01
The life and operation of cells involve many physiological processes that take place over fast timescales of milliseconds to minutes. Genetically-encoded technologies for driving or suppressing specific fast physiological processes in intact cells, perhaps embedded within intact tissues in living organisms, are critical for the ability to understand how these physiological processes contribute to emergent cellular and organismal functions and behaviors. Such “synthetic physiology” tools are o...
Exchange splitting of the interaction energy and the multipole expansion of the wave function
Gniewek, Piotr, E-mail: pgniewek@tiger.chem.uw.edu.pl; Jeziorski, Bogumił, E-mail: jeziorsk@chem.uw.edu.pl [Faculty of Chemistry, University of Warsaw, Pasteura 1, 02-093 Warsaw (Poland)
2015-10-21
The exchange splitting J of the interaction energy of the hydrogen atom with a proton is calculated using the conventional surface-integral formula J{sub surf}[Φ], the volume-integral formula of the symmetry-adapted perturbation theory J{sub SAPT}[Φ], and a variational volume-integral formula J{sub var}[Φ]. The calculations are based on the multipole expansion of the wave function Φ, which is divergent for any internuclear distance R. Nevertheless, the resulting approximations to the leading coefficient j{sub 0} in the large-R asymptotic series J(R) = 2e{sup −R−1}R(j{sub 0} + j{sub 1}R{sup −1} + j{sub 2}R{sup −2} + ⋯) converge with the rate corresponding to the convergence radii equal to 4, 2, and 1 when the J{sub var}[Φ], J{sub surf}[Φ], and J{sub SAPT}[Φ] formulas are used, respectively. Additionally, we observe that also the higher j{sub k} coefficients are predicted correctly when the multipole expansion is used in the J{sub var}[Φ] and J{sub surf}[Φ] formulas. The symmetry adapted perturbation theory formula J{sub SAPT}[Φ] predicts correctly only the first two coefficients, j{sub 0} and j{sub 1}, gives a wrong value of j{sub 2}, and diverges for higher j{sub n}. Since the variational volume-integral formula can be easily generalized to many-electron systems and evaluated with standard basis-set techniques of quantum chemistry, it provides an alternative for the determination of the exchange splitting and the exchange contribution of the interaction potential in general.
Harmonic Decomposition of Orbit Data for Multipole Analysis
Yang Ming Jen
2005-01-01
The unprecedented position resolution provided by the newly commissioned Recycer BPM system is opening up a new chapter of beam based multipole analysis at Fermilab. The closed orbit data, taken with circulating beam and averaged over many consecutive turns, has been shown to have the resolution of a few microns. The result of polynomial fit to BPM position data, as a function of dipole kick sizes, is used to separate orbit data into first, second, and third order. Combining both the in-plane and cross-plane orbit data it is possible to determine the multipole content within each half cell. This paper presents the algorithm behind the procedure, the data collected from the Fermilab Recycler Ring, and the final analysis result.
Multipole moments of the generalized Quevedo-Mashhoon metric
Frutos-Alfaro, Francisco
2016-01-01
Among the known exact solutions of Einsteins vacuum field equations the generalized Quevedo-Mashhoon (QM) metric with parameters M, a and qn (n being an even integer) might be the most suitable one for the description of the exterior gravitational field of some real non-collapsed body. Originally this metric was derived by some Hoenselaers-Kinnersley-Xanthopoulos transformation with the Erez-Rosen metric as seed metric. We verified directly the validity of the corresponding Ernst equations by means of a REDUCE program. By means of the Fodor-Hoenselaers-Perjes formalism we computed the first ten mass and spin multipole moments of the QM metric that obviously have not been published before. Corresponding moments were derived for the static vacuum solutions of Gutsunaev-Manko and Hernandez-Martin. A direct comparison between the multipole moments of these non-isometric spacetimes is given.
IBIS: A hollow-cathode multipole boundary ion source
The plasma production and containment system for a high-power continuously operating magnetic multipole ion source has been designed and constructed. Preliminary tests on this system prior to high voltage extraction of large beams indicate advantageous performance for neutral-beam injection applications. The source has produced 80 A to the extractor region at 0.33 A/cm2 with a discharge of 330 A at 80 V. Density uniformity is better than 1% over a 16-cm diameter, dropping to -4% at 18 cm, with plasma noise of less than 3%. Gas utilizaion efficiency and atomic (H+) species output are anticipated to be high due to a source length of 40 cm. This quiet efficient performance is attributed to the use of a hollow-tube LaB6 cathode and an improved magnetic multipole confinement system