A Multi-signal Variant for the GPU-based Parallelization of Growing Self-Organizing Networks

Among the many possible approaches for the parallelization of self-organizing networks, and in particular of growing self-organizing networks, perhaps the most common one is producing an optimized, parallel implementation of the standard sequential algorithms reported in the literature. In this paper we explore an alternative approach, based on a new algorithm variant specifically designed to match the features of the large-scale, fine-grained parallelism of GPUs, in which multiple input signals are processed at once. Comparative tests have been performed, using both parallel and sequential implementations of the new algorithm variant, in particular for a growing self-organizing network that reconstructs surfaces from point clouds. The experimental results show that this approach allows harnessing in a more effective way the intrinsic parallelism that the self-organizing networks algorithms seem intuitively to suggest, obtaining better performances even with networks of smaller size.Comment: 17 page

Rapid planetesimal formation in turbulent circumstellar discs

The initial stages of planet formation in circumstellar gas discs proceed via dust grains that collide and build up larger and larger bodies (Safronov 1969). How this process continues from metre-sized boulders to kilometre-scale planetesimals is a major unsolved problem (Dominik et al. 2007): boulders stick together poorly (Benz 2000), and spiral into the protostar in a few hundred orbits due to a head wind from the slower rotating gas (Weidenschilling 1977). Gravitational collapse of the solid component has been suggested to overcome this barrier (Safronov 1969, Goldreich & Ward 1973, Youdin & Shu 2002). Even low levels of turbulence, however, inhibit sedimentation of solids to a sufficiently dense midplane layer (Weidenschilling & Cuzzi 1993, Dominik et al. 2007), but turbulence must be present to explain observed gas accretion in protostellar discs (Hartmann 1998). Here we report the discovery of efficient gravitational collapse of boulders in locally overdense regions in the midplane. The boulders concentrate initially in transient high pressures in the turbulent gas (Johansen, Klahr, & Henning 2006), and these concentrations are augmented a further order of magnitude by a streaming instability (Youdin & Goodman 2005, Johansen, Henning, & Klahr 2006, Johansen & Youdin 2007) driven by the relative flow of gas and solids. We find that gravitationally bound clusters form with masses comparable to dwarf planets and containing a distribution of boulder sizes. Gravitational collapse happens much faster than radial drift, offering a possible path to planetesimal formation in accreting circumstellar discs.Comment: To appear in Nature (30 August 2007 issue). 18 pages (in referee mode), 3 figures. Supplementary Information can be found at 0708.389

Investigation of initiation of gigantic jets connecting thunderclouds to the ionosphere

The initiation of giant electrical discharges called as "gigantic jets" connecting thunderclouds to the ionosphere is investigated by numerical simulation method in this paper. Using similarity relations, the triggering conditions of streamer formation in laboratory situations are extended to form a criterion of initiation of gigantic jets. The energy source causing a gigantic jet is considered due to the quasi-electrostatic field generated by thunderclouds. The electron dynamics from ionization threshold to streamer initiation are simulated by the Monte Carlo technique. It is found that gigantic jets are initiated at a height of ~18-24 km. This is in agreement with the observations. The method presented in this paper could be also applied to the analysis of the initiation of other discharges such as blue jets and red sprites.Comment: 12th International Congress on Plasma Physics, 25-29 October 2004, Nice (France

Desynchronization and Wave Pattern Formation in MPI-Parallel and Hybrid Memory-Bound Programs

Analytic, first-principles performance modeling of distributed-memory parallel codes is notoriously imprecise. Even for applications with extremely regular and homogeneous compute-communicate phases, simply adding communication time to computation time does often not yield a satisfactory prediction of parallel runtime due to deviations from the expected simple lockstep pattern caused by system noise, variations in communication time, and inherent load imbalance. In this paper, we highlight the specific cases of provoked and spontaneous desynchronization of memory-bound, bulk-synchronous pure MPI and hybrid MPI+OpenMP programs. Using simple microbenchmarks we observe that although desynchronization can introduce increased waiting time per process, it does not necessarily cause lower resource utilization but can lead to an increase in available bandwidth per core. In case of significant communication overhead, even natural noise can shove the system into a state of automatic overlap of communication and computation, improving the overall time to solution. The saturation point, i.e., the number of processes per memory domain required to achieve full memory bandwidth, is pivotal in the dynamics of this process and the emerging stable wave pattern. We also demonstrate how hybrid MPI-OpenMP programming can prevent desirable desynchronization by eliminating the bandwidth bottleneck among processes. A Chebyshev filter diagonalization application is used to demonstrate some of the observed effects in a realistic setting.Comment: 18 pages, 8 figure

Calculations of binding affinity between C8-substituted GTP analogs and the bacterial cell-division protein FtsZ

The FtsZ protein is a self-polymerizing GTPase that plays a central role in bacterial cell division. Several C8-substituted GTP analogs are known to inhibit the polymerization of FtsZ by competing for the same binding site as its endogenous activating ligand GTP. Free energy calculations of the relative binding affinities to FtsZ for a set of five C8-substituted GTP analogs were performed. The calculated values agree well with the available experimental data, and the main contribution to the free energy differences is determined to be the conformational restriction of the ligands. The dihedral angle distributions around the glycosidic bond of these compounds in water are known to vary considerably depending on the physicochemical properties of the substituent at C8. However, within the FtsZ protein, this substitution has a negligible influence on the dihedral angle distributions, which fall within the narrow range of −140° to −90° for all investigated compounds. The corresponding ensemble average of the coupling constants 3J(C4,H1′) is calculated to be 2.95 ± 0.1 Hz. The contribution of the conformational selection of the GTP analogs upon binding was quantified from the corresponding populations. The obtained restraining free energy values follow the same trend as the relative binding affinities to FtsZ, indicating their dominant contribution

Single-Molecule Force Spectroscopy: Experiments, Analysis, and Simulations

International audienceThe mechanical properties of cells and of subcellular components are important to obtain a mechanistic molecular understanding of biological processes. The quantification of mechanical resistance of cells and biomolecules using biophysical methods matured thanks to the development of nanotechnologies such as optical and magnetic tweezers, the biomembrane force probe and atomic force microscopy (AFM). The quantitative nature of force spectroscopy measurements has converted AFM into a valuable tool in biophysics. Force spectroscopy allows the determination of the forces required to unfold protein domains and to disrupt individual receptor/ligand bonds. Molecular simulation as a computational microscope allows investigation of similar biological processes with an atomistic detail. In this chapter, we first provide a step-by-step protocol of force spectroscopy including sample preparation, measurement and analysis of force spectroscopy using AFM and its interpretation in terms of available theories. Next, we present the background for molecular dynamics (MD) simulations focusing on steered molecular dynamics (SMD) and the importance of bridging of computational tools with experimental technique

Unified Homogenization Theory for Magnetoinductive and Electromagnetic Waves in Split Ring Metamaterials

A unified homogenization procedure for split ring metamaterials taking into account time and spatial dispersion is introduced. The procedure is based on two coupled systems of equations. The first one comes from an approximation of the metamaterial as a cubic arrangement of coupled LC circuits, giving the relation between currents and local magnetic field. The second equation comes from macroscopic Maxwell equations, and gives the relation between the macroscopic magnetic field and the average magnetization of the metamaterial. It is shown that electromagnetic and magnetoinductive waves propagating in the metamaterial are obtained from this analysis. Therefore, the proposed time and spatially dispersive permeability accounts for the characterization of the complete spectrum of waves of the metamaterial. Finally, it is shown that the proposed theory is in good quantitative and qualitative agreement with full wave simulations.Comment: 4 pages, 3 figure

All-atom/coarse-grained hybrid predictions of distribution coefficients in SAMPL5

We present blind predictions submitted to the SAMPL5 challenge on calculating distribution coefficients. The predictions were based on estimating the solvation free energies in water and cyclohexane of the 53 compounds in the challenge. These free energies were computed using alchemical free energy simulations based on a hybrid all-atom/coarse-grained model. The compounds were treated with the general Amber force field, whereas the solvent molecules were treated with the Elba coarse-grained model. Considering the simplicity of the solvent model and that we approximate the distribution coefficient with the partition coefficient of the neutral species, the predictions are of good accuracy. The correlation coefficient, R is 0.64, 82 % of the predictions have the correct sign and the mean absolute deviation is 1.8 log units. This is on a par with or better than the other simulation-based predictions in the challenge. We present an analysis of the deviations to experiments and compare the predictions to another submission that used all-atom solvent. ELECTRONIC SUPPLEMENTARY MATERIAL: The online version of this article (doi:10.1007/s10822-016-9926-z) contains supplementary material, which is available to authorized users

Dynamics of Disks and Warps

This chapter reviews theoretical work on the stellar dynamics of galaxy disks. All the known collective global instabilities are identified, and their mechanisms described in terms of local wave mechanics. A detailed discussion of warps and other bending waves is also given. The structure of bars in galaxies, and their effect on galaxy evolution, is now reasonably well understood, but there is still no convincing explanation for their origin and frequency. Spiral patterns have long presented a special challenge, and ideas and recent developments are reviewed. Other topics include scattering of disk stars and the survival of thin disks.Comment: Chapter accepted to appear in Planets, Stars and Stellar Systems, vol 5, ed G. Gilmore. 32 pages, 17 figures. Includes minor corrections made in proofs. Uses emulateapj.st