BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

Gerald Schubert

Gerald Schubert contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
cond-mat.str-elastro-ph.EPgr-qccond-mat.dis-nncond-mat.mes-hallPerformanceDistributed, Parallel, and Cluster Computing

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
17works
0followers
7topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

17 published item(s)

preprint2016arXiv

Deep space experiment to measure $G$

Responding to calls from the National Science Foundation (NSF) for new proposals to measure the gravitational constant $G$, we offer an interesting experiment in deep space employing the classic gravity train mechanism. Our setup requires three bodies: a larger layered solid sphere with a cylindrical hole through its center, a much smaller retroreflector which will undergo harmonic motion within the hole and a host spacecraft with laser ranging capabilities to measure round trip light-times to the retroreflector but ultimately separated a significant distance away from the sphere-retroreflector apparatus. Measurements of the period of oscillation of the retroreflector in terms of host spacecraft clock time using existing technology could give determinations of $G$ nearly three orders of magnitude more accurate than current measurements here on Earth. However, significant engineering advances in the release mechanism of the apparatus from the host spacecraft will likely be necessary. Issues with regard to the stability of the system are briefly addressed.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Measurements of Newton's gravitational constant and the length of day

About a dozen measurements of Newton's gravitational constant, G, since 1962 have yielded values that differ by far more than their reported random plus systematic errors. We find that these values for G are oscillatory in nature, with a period of P = 5.899 +/- 0.062 yr, an amplitude of (1.619 +/- 0.103) x 10^{-14} m^3 kg^{-1} s^{-2}, and mean-value crossings in 1994 and 1997. However, we do not suggest that G is actually varying by this much, this quickly, but instead that something in the measurement process varies. Of other recently reported results, to the best of our knowledge, the only measurement with the same period and phase is the Length of Day (LOD - defined as a frequency measurement such that a positive increase in LOD values means slower Earth rotation rates and therefore longer days). The aforementioned period is also about half of a solar activity cycle, but the correlation is far less convincing. The 5.9 year periodic signal in LOD has previously been interpreted as due to fluid core motions and inner-core coupling. We report the G/LOD correlation, whose statistical significance is 0.99764 assuming no difference in phase, without claiming to have any satisfactory explanation for it. Least unlikely, perhaps, are currents in the Earth's fluid core that change both its moment of inertia (affecting LOD) and the circumstances in which the Earth-based experiments measure G. In this case, there might be correlations with terrestrial magnetic field measurements.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Charge carrier localisation in disordered graphene nanoribbons

We study the electronic properties of actual-size graphene nanoribbons subjected to substitutional disorder particularly with regard to the experimentally observed metal-insulator transition. Calculating the local, mean and typical density of states, as well as the time-evolution of the particle density we comment on a possible disorder-induced localisation of charge carriers at and close to the Dirac point within a percolation transition scenario.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

The fate of topological-insulator surface states under strong disorder

Three-dimensional topological insulators feature Dirac-like surface states which are topologically protected against the influence of weak quenched disorder. Here we investigate the effect of surface disorder beyond the weak-disorder limit using large-scale numerical simulations. We find two qualitatively distinct regimes: Moderate disorder destroys the Dirac cone and induces diffusive metallic behavior at the surface. Even more remarkably, for strong surface disorder a Dirac cone reappears, as new weakly disordered "surface" states emerge in the sample beneath the disordered surface layer, which can be understood in terms of an interface between a topological and an Anderson insulator. Together, this demonstrates the drastic effect of disorder on topological surface states, which cannot be captured within effective two-dimensional models for the surface states alone.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Hybrid-parallel sparse matrix-vector multiplication with explicit communication overlap on current multicore-based systems

We evaluate optimized parallel sparse matrix-vector operations for several representative application areas on widespread multicore-based cluster configurations. First the single-socket baseline performance is analyzed and modeled with respect to basic architectural properties of standard multicore chips. Beyond the single node, the performance of parallel sparse matrix-vector operations is often limited by communication overhead. Starting from the observation that nonblocking MPI is not able to hide communication cost using standard MPI implementations, we demonstrate that explicit overlap of communication and computation can be achieved by using a dedicated communication thread, which may run on a virtual core. Moreover we identify performance benefits of hybrid MPI/OpenMP programming due to improved load balancing even without explicit communication overlap. We compare performance results for pure MPI, the widely used "vector-like" hybrid programming strategies, and explicit overlap on a modern multicore-based cluster and a Cray XE6 system.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Jupiter's Moment of Inertia: A Possible Determination by JUNO

The moment of inertia of a giant planet reveals important information about the planet's internal density structure and this information is not identical to that contained in the gravitational moments. The forthcoming Juno mission to Jupiter might determine Jupiter's normalized moment of inertia NMoI=C/MR^2 by measuring Jupiter's pole precession and the Lense-Thirring acceleration of the spacecraft (C is the axial moment of inertia, and M and R are Jupiter's mass and mean radius, respectively). We investigate the possible range of NMoI values for Jupiter based on its measured gravitational field using a simple core/envelope model of the planet assuming that J_2 and J_4 are perfectly known and are equal to their measured values. The model suggests that for fixed values of J_2 and J_4 a range of NMOI values between 0.2629 and 0.2645 can be found. The Radau-Darwin relation gives a NMoI value that is larger than the model values by less than 1%. A low NMoI of ~ 0.236, inferred from a dynamical model (Ward & Canup, 2006, ApJ, 640, L91) is inconsistent with this range, but the range is model dependent. Although we conclude that the NMoI is tightly constrained by the gravity coefficients, a measurement of Jupiter's NMoI to a few tenths of percent by Juno could provide an important constraint on Jupiter's internal structure. We carry out a simplified assessment of the error involved in Juno's possible determination of Jupiter's NMoI.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Metal-to-insulator transition and electron-hole puddle formation in disordered graphene nanoribbons

The experimentally observed metal-to-insulator transition in hydrogenated graphene is numerically confirmed for actual sized graphene samples and realistic impurity concentrations. The eigenstates of our tight-binding model with substitutional disorder corroborate the formation of electron-hole-puddles with characteristic length scales comparable to the ones found in experiments. The puddles cause charge inhomogeneities and tend to suppress Anderson localization. Even though, monitoring the charge carrier quantum dynamics and performing a finite-size scaling of the local density of states distribution, we find strong evidence for the existence of localized states in graphene nanoribbons with short-range but also correlated long-range disorder.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Shapes and gravitational fields of rotating two-layer Maclaurin ellipsoids: Application to planets and satellites

The exact solution for the shape and gravitational field of a rotating two-layer Maclaurin ellipsoid of revolution is compared with predictions of the theory of figures up to third order in the small rotational parameter of the theory of figures. An explicit formula is derived for the external gravitational coefficient $J_2$ of the exact solution. A new approach to the evaluation of the theory of figures based on numerical integration of ordinary differential equations is presented. The classical Radau-Darwin formula is found not to be valid for the rotational parameter ε_2 = Ω^2/(2πGρ_2) >= 0.17 since the formula then predicts a surface eccentricity that is smaller than the eccentricity of the core-envelope boundary. Interface eccentricity must be smaller than surface eccentricity. In the formula for $ε_2$, $Ω$ is the angular velocity of the two-layer body, $ρ_2$ is the density of the outer layer, and G is the gravitational constant. For an envelope density of 3000 kg m^(-3) the failure of the Radau-Darwin formula corresponds to a rotation period of about 3 hr. Application of the exact solution and the theory of figures is made to models of Earth, Mars, Uranus, and Neptune. The two-layer model with constant densities in the layers can provide realistic approximations to terrestrial planets and icy outer planet satellites. The two-layer model needs to be generalized to allow for a continuous envelope (outer layer) radial density profile in order to realistically model a gas or ice giant planet.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Distribution of the local density of states as a criterion for Anderson localization: Numerically exact results for various lattices in two and three dimensions

Numerical approaches to Anderson localization face the problem of having to treat large localization lengths while being restricted to finite system sizes. We show that by finite-size scaling of the probability distribution of the local density of states (LDOS) this long-standing problem can be overcome. To this end we reexamine the approach, propose numerical refinements, and apply it to study the dependence of the distribution of the LDOS on the dimensionality and coordination number of the lattice. Particular attention is given to the graphene lattice. We show that the system-size dependence of the LDOS distribution is indeed an unambiguous sign of Anderson localization, irrespective of the dimension and lattice structure. The numerically exact LDOS data obtained by us agree with a log-normal distribution over up to ten orders of magnitude and thereby fulfill a nontrivial symmetry relation previously derived for the non-linear $σ$-model.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Interior Models of Uranus and Neptune

'Empirical' models (pressure vs. density) of Uranus and Neptune interiors constrained by the gravitational coefficients J_2, J_4, the planetary radii and masses, and Voyager solid-body rotation periods are presented. The empirical pressure-density profiles are then interpreted in terms of physical equations of state of hydrogen, helium, ice (H_2O), and rock (SiO_2) to test the physical plausibility of the models. The compositions of Uranus and Neptune are found to be similar with somewhat different distributions of the high-Z material. The big difference between the two planets is that Neptune requires a non-solar envelope while Uranus is best matched with a solar composition envelope. Our analysis suggests that the heavier elements in both Uranus' and Neptune's interior might increase gradually towards the planetary centers. Indeed it is possible to fit the gravitational moments without sharp compositional transitions.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Parallel sparse matrix-vector multiplication as a test case for hybrid MPI+OpenMP programming

We evaluate optimized parallel sparse matrix-vector operations for two representative application areas on widespread multicore-based cluster configurations. First the single-socket baseline performance is analyzed and modeled with respect to basic architectural properties of standard multicore chips. Going beyond the single node, parallel sparse matrix-vector operations often suffer from an unfavorable communication to computation ratio. Starting from the observation that nonblocking MPI is not able to hide communication cost using standard MPI implementations, we demonstrate that explicit overlap of communication and computation can be achieved by using a dedicated communication thread, which may run on a virtual core. We compare our approach to pure MPI and the widely used "vector-like" hybrid programming strategy.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Uranus and Neptune: Shape and Rotation

Both Uranus and Neptune are thought to have strong zonal winds with velocities of several hundred meters per second. These wind velocities, however, assume solid-body rotation periods based on Voyager 2 measurements of periodic variations in the planets' radio signals and of fits to the planets' magnetic fields; 17.24h and 16.11h for Uranus and Neptune, respectively. The realization that the radio period of Saturn does not represent the planet's deep interior rotation and the complexity of the magnetic fields of Uranus and Neptune raise the possibility that the Voyager 2 radio and magnetic periods might not represent the deep interior rotation periods of the ice giants. Moreover, if there is deep differential rotation within Uranus and Neptune no single solid-body rotation period could characterize the bulk rotation of the planets. We use wind and shape data to investigate the rotation of Uranus and Neptune. The shapes (flattening) of the ice giants are not measured, but only inferred from atmospheric wind speeds and radio occultation measurements at a single latitude. The inferred oblateness values of Uranus and Neptune do not correspond to bodies rotating with the Voyager rotation periods. Minimization of wind velocities or dynamic heights of the 1 bar isosurfaces, constrained by the single occultation radii and gravitational coefficients of the planets, leads to solid-body rotation periods of ~16.58h for Uranus and ~17.46h for Neptune. Uranus might be rotating faster and Neptune slower than Voyager rotation speeds. We derive shapes for the planets based on these rotation rates. Wind velocities with respect to these rotation periods are essentially identical on Uranus and Neptune and wind speeds are slower than previously thought. Alternatively, if we interpret wind measurements in terms of differential rotation on cylinders there are essentially no residual atmospheric winds.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Anderson disorder in graphene nanoribbons: A local distribution approach

Disorder effects strongly influence the transport properties of graphene based nanodevices even to the point of Anderson localization. Focusing on the local density of states and its distribution function, we analyze the localization properties of actual size graphene nanoribbons. In particular we determine the time evolution and localization length of the single particle wave function in dependence on the ribbon extension and edge geometry, as well as on the disorder type and strength.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Numerical approaches to time evolution of complex quantum systems

We examine several numerical techniques for the calculation of the dynamics of quantum systems. In particular, we single out an iterative method which is based on expanding the time evolution operator into a finite series of Chebyshev polynomials. The Chebyshev approach benefits from two advantages over the standard time-integration Crank-Nicholson scheme: speedup and efficiency. Potential competitors are semiclassical methods such as the Wigner-Moyal or quantum tomographic approaches. We outline the basic concepts of these techniques and benchmark their performance against the Chebyshev approach by monitoring the time evolution of a Gaussian wave packet in restricted one-dimensional (1D) geometries. Thereby the focus is on tunnelling processes and the motion in anharmonic potentials. Finally we apply the prominent Chebyshev technique to two highly non-trivial problems of current interest: (i) the injection of a particle in a disordered 2D graphene nanoribbon and (ii) the spatiotemporal evolution of polaron states in finite quantum systems. Here, depending on the disorder/electron-phonon coupling strength and the device dimensions, we observe transmission or localisation of the matter wave.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Performance limitations for sparse matrix-vector multiplications on current multicore environments

The increasing importance of multicore processors calls for a reevaluation of established numerical algorithms in view of their ability to profit from this new hardware concept. In order to optimize the existent algorithms, a detailed knowledge of the different performance-limiting factors is mandatory. In this contribution we investigate sparse matrix-vector multiplication, which is the dominant operation in many sparse eigenvalue solvers. Two conceptually different storage schemes and computational kernels have been conceived in the past to target cache-based and vector architectures, respectively. Starting from a series of microbenchmarks we apply the gained insight on optimized sparse MVM implementations, whose serial and OpenMP-parallel performance we review on state-of-the-art multicore systems.

0 citations0 reviewsNo code linkedOpen access
preprint2005arXiv

Comparative numerical study of Anderson localization in disordered electron systems

Taking into account that a proper description of disordered systems should focus on distribution functions, the authors develop a powerful numerical scheme for the determination of the probability distribution of the local density of states (LDOS), which is based on a Chebyshev expansion with kernel polynomial refinement and allows the study of large finite clusters (up to $100^3$). For the three-dimensional Anderson model it is demonstrated that the distribution of the LDOS shows a significant change at the disorder induced delocalisation-localisation transition. Consequently, the so-called typical density of states, defined as the geometric mean of the LDOS, emerges as a natural order parameter. The calculation of the phase diagram of the Anderson model proves the efficiency and reliability of the proposed approach in comparison to other localisation criteria, which rely, e.g., on the decay of the wavefunction or the inverse participation number.

0 citations0 reviewsNo code linkedOpen access