-
Third-order intrinsic alignment of SDSS BOSS LOWZ galaxies
Authors:
Laila Linke,
Susan Pyne,
Benjamin Joachimi,
Christos Georgiou,
Kai Hoffmann,
Rachel Mandelbaum,
Sukhdeep Singh
Abstract:
Cosmic shear is a powerful probe of cosmology, but it is affected by the intrinsic alignment (IA) of galaxy shapes with the large-scale structure. Upcoming surveys like Euclid and Vera C. Rubin Observatory's Legacy Survey of Space and Time (LSST) require an accurate understanding of IA, particularly for higher-order cosmic shear statistics that are vital for extracting the most cosmological inform…
▽ More
Cosmic shear is a powerful probe of cosmology, but it is affected by the intrinsic alignment (IA) of galaxy shapes with the large-scale structure. Upcoming surveys like Euclid and Vera C. Rubin Observatory's Legacy Survey of Space and Time (LSST) require an accurate understanding of IA, particularly for higher-order cosmic shear statistics that are vital for extracting the most cosmological information. In this paper, we report the first detection of third-order IA correlations using the LOWZ galaxy sample from the Sloan Digital Sky Survey (SDSS) Baryon Oscillation Spectroscopic Survey (BOSS). We compare our measurements with predictions from the MICE cosmological simulation and an analytical NLA-inspired model informed by second-order correlations. We also explore the dependence of the third-order correlation on the galaxies' luminosity. We find that the amplitude $A_\mathrm{IA}$ of the IA signal is non-zero at the $4.7σ$ ($7.6σ$) level for scales between $6 h^{-1} \mathrm{Mpc}$ ($1 h^{-1} \mathrm{Mpc}$) and $20 h^{-1} \mathrm{Mpc}$. For scales above $6 h^{-1}\mathrm{Mpc}$ the inferred AIA agrees both with the prediction from the simulation and estimates from second-order statistics within $1σ$ but deviations arise at smaller scales. Our results demonstrate the feasibility of measuring third-order IA correlations and using them for constraining IA models. The agreement between second- and third-order IA constraints also opens the opportunity for a consistent joint analysis and IA self-calibration, promising tighter parameter constraints for upcoming cosmological surveys.
△ Less
Submitted 7 June, 2024;
originally announced June 2024.
-
An iterative closest point algorithm for marker-free 3D shape registration of continuum robots
Authors:
Matthias K. Hoffmann,
Julian Mühlenhoff,
Zhaoheng Ding,
Thomas Sattel,
Kathrin Flaßkamp
Abstract:
Continuum robots have emerged as a promising technology in the medical field due to their potential of accessing deep sited locations of the human body with low surgical trauma. When deriving physics-based models for these robots, evaluating the models poses a significant challenge due to the difficulty in accurately measuring their intricate shapes. In this work, we present an optimization based…
▽ More
Continuum robots have emerged as a promising technology in the medical field due to their potential of accessing deep sited locations of the human body with low surgical trauma. When deriving physics-based models for these robots, evaluating the models poses a significant challenge due to the difficulty in accurately measuring their intricate shapes. In this work, we present an optimization based 3D shape registration algorithm for estimation of the backbone shape of slender continuum robots as part of a pho togrammetric measurement. Our approach to estimating the backbones optimally matches a parametric three-dimensional curve to images of the robot. Since we incorporate an iterative closest point algorithm into our method, we do not need prior knowledge of the robots position within the respective images. In our experiments with artificial and real images of a concentric tube continuum robot, we found an average maximum deviation of the reconstruction from simulation data of 0.665 mm and 0.939 mm from manual measurements. These results show that our algorithm is well capable of producing high accuracy positional data from images of continuum robots.
△ Less
Submitted 24 May, 2024;
originally announced May 2024.
-
Euclid. V. The Flagship galaxy mock catalogue: a comprehensive simulation for the Euclid mission
Authors:
Euclid Collaboration,
F. J. Castander,
P. Fosalba,
J. Stadel,
D. Potter,
J. Carretero,
P. Tallada-Crespí,
L. Pozzetti,
M. Bolzonella,
G. A. Mamon,
L. Blot,
K. Hoffmann,
M. Huertas-Company,
P. Monaco,
E. J. Gonzalez,
G. De Lucia,
C. Scarlata,
M. -A. Breton,
L. Linke,
C. Viglione,
S. -S. Li,
Z. Zhai,
Z. Baghkhani,
K. Pardede,
C. Neissner
, et al. (344 additional authors not shown)
Abstract:
We present the Flagship galaxy mock, a simulated catalogue of billions of galaxies designed to support the scientific exploitation of the Euclid mission. Euclid is a medium-class mission of the European Space Agency optimised to determine the properties of dark matter and dark energy on the largest scales of the Universe. It probes structure formation over more than 10 billion years primarily from…
▽ More
We present the Flagship galaxy mock, a simulated catalogue of billions of galaxies designed to support the scientific exploitation of the Euclid mission. Euclid is a medium-class mission of the European Space Agency optimised to determine the properties of dark matter and dark energy on the largest scales of the Universe. It probes structure formation over more than 10 billion years primarily from the combination of weak gravitational lensing and galaxy clustering data. The breath of Euclid's data will also foster a wide variety of scientific analyses. The Flagship simulation was developed to provide a realistic approximation to the galaxies that will be observed by Euclid and used in its scientific analyses. We ran a state-of-the-art N-body simulation with four trillion particles, producing a lightcone on the fly. From the dark matter particles, we produced a catalogue of 16 billion haloes in one octant of the sky in the lightcone up to redshift z=3. We then populated these haloes with mock galaxies using a halo occupation distribution and abundance matching approach, calibrating the free parameters of the galaxy mock against observed correlations and other basic galaxy properties. Modelled galaxy properties include luminosity and flux in several bands, redshifts, positions and velocities, spectral energy distributions, shapes and sizes, stellar masses, star formation rates, metallicities, emission line fluxes, and lensing properties. We selected a final sample of 3.4 billion galaxies with a magnitude cut of H_E<26, where we are complete. We have performed a comprehensive set of validation tests to check the similarity to observational data and theoretical models. In particular, our catalogue is able to closely reproduce the main characteristics of the weak lensing and galaxy clustering samples to be used in the mission's main cosmological analysis. (abridged)
△ Less
Submitted 22 May, 2024;
originally announced May 2024.
-
Robust Tumor Segmentation with Hyperspectral Imaging and Graph Neural Networks
Authors:
Mayar Lotfy,
Anna Alperovich,
Tommaso Giannantonio,
Bjorn Barz,
Xiaohan Zhang,
Felix Holm,
Nassir Navab,
Felix Boehm,
Carolin Schwamborn,
Thomas K. Hoffmann,
Patrick J. Schuler
Abstract:
Segmenting the boundary between tumor and healthy tissue during surgical cancer resection poses a significant challenge. In recent years, Hyperspectral Imaging (HSI) combined with Machine Learning (ML) has emerged as a promising solution. However, due to the extensive information contained within the spectral domain, most ML approaches primarily classify individual HSI (super-)pixels, or tiles, wi…
▽ More
Segmenting the boundary between tumor and healthy tissue during surgical cancer resection poses a significant challenge. In recent years, Hyperspectral Imaging (HSI) combined with Machine Learning (ML) has emerged as a promising solution. However, due to the extensive information contained within the spectral domain, most ML approaches primarily classify individual HSI (super-)pixels, or tiles, without taking into account their spatial context. In this paper, we propose an improved methodology that leverages the spatial context of tiles for more robust and smoother segmentation. To address the irregular shapes of tiles, we utilize Graph Neural Networks (GNNs) to propagate context information across neighboring regions. The features for each tile within the graph are extracted using a Convolutional Neural Network (CNN), which is trained simultaneously with the subsequent GNN. Moreover, we incorporate local image quality metrics into the loss function to enhance the training procedure's robustness against low-quality regions in the training images. We demonstrate the superiority of our proposed method using a clinical ex vivo dataset consisting of 51 HSI images from 30 patients. Despite the limited dataset, the GNN-based model significantly outperforms context-agnostic approaches, accurately distinguishing between healthy and tumor tissues, even in images from previously unseen patients. Furthermore, we show that our carefully designed loss function, accounting for local image quality, results in additional improvements. Our findings demonstrate that context-aware GNN algorithms can robustly find tumor demarcations on HSI images, ultimately contributing to better surgery success and patient outcome.
△ Less
Submitted 20 November, 2023;
originally announced November 2023.
-
High precision calibration setup for loss measurements in electrical steel sheets
Authors:
K. Pfnuer,
J. Luedke,
K. Hoffmann,
F. Weickert
Abstract:
We present details on the current measurement setup at PTB used for high precision loss calibrations in the frequency range 50 Hz to 1 kHz. A combination of analog and digital feedback control is utilized in accordance with the standard. A detailed measurement uncertainty (MU) analysis based on a systematic model equation is presented and inter-dependencies of model parameters are discussed. Exper…
▽ More
We present details on the current measurement setup at PTB used for high precision loss calibrations in the frequency range 50 Hz to 1 kHz. A combination of analog and digital feedback control is utilized in accordance with the standard. A detailed measurement uncertainty (MU) analysis based on a systematic model equation is presented and inter-dependencies of model parameters are discussed. Experimental results obtained at 50 Hz on NO and GO Epstein samples show excellent agreement between statistical and systematic MU estimation and confirm the MU model analysis. Furthermore, we investigate the influence of external parameters on the loss measurements, like sample loading scheme and the value of maximum demagnetization polarization.
△ Less
Submitted 26 October, 2023;
originally announced November 2023.
-
Multi-Objective Model-Predictive Control for Dielectric Elastomer Wave Harvesters
Authors:
Matthias K. Hoffmann,
Lennart Heib,
Gianluca Rizzello,
Giacomo Moretti,
Kathrin Flaßkamp
Abstract:
This contribution deals with multi-objective model-predictive control (MPC) of a wave energy converter (WEC) device concept, which can harvest energy from sea waves using a dielectric elastomer generator (DEG) power take-off system. We aim to maximise the extracted energy through control while minimising the accumulated damage to the DEG. With reference to system operation in stochastic waves, we…
▽ More
This contribution deals with multi-objective model-predictive control (MPC) of a wave energy converter (WEC) device concept, which can harvest energy from sea waves using a dielectric elastomer generator (DEG) power take-off system. We aim to maximise the extracted energy through control while minimising the accumulated damage to the DEG. With reference to system operation in stochastic waves, we first generate ground truth solutions by solving an optimal control problem, and we analyse the MPC performance to determine a prediction horizon that trades off accuracy and efficiency for computation. Fixed weights in the MPC scheme can produce unpredictable costs for variable sea condition, meaning the average rate of cost accumulation can vary vastly. To steer this cost growth, we propose a heuristic to adapt the algorithm by changing the weighting of the cost functions using for fulfilling the long-time goal of accumulating a small enough damage in a fixed time. A simulated case-study is presented in order to evaluate the performance of the proposed MPC framework and the weight-adaptation algorithm. The proposed heuristic proves to be able to limit the amount of accumulated damage while remaining close to (or even improving) the energy yield obtained with a comparable fixed-weight MPC.
△ Less
Submitted 7 December, 2022;
originally announced December 2022.
-
Path Planning for Concentric Tube Robots: a Toolchain with Application to Stereotactic Neurosurgery
Authors:
Matthias K. Hoffmann,
Willem Esterhuizen,
Karl Worthmann,
Kathrin Flaßkamp
Abstract:
We present a toolchain for solving path planning problems for concentric tube robots through obstacle fields. First, ellipsoidal sets representing the target area and obstacles are constructed from labelled point clouds. Then, the nonlinear and highly nonconvex optimal control problem is solved by introducing a homotopy on the obstacle positions where at one extreme of the parameter the obstacles…
▽ More
We present a toolchain for solving path planning problems for concentric tube robots through obstacle fields. First, ellipsoidal sets representing the target area and obstacles are constructed from labelled point clouds. Then, the nonlinear and highly nonconvex optimal control problem is solved by introducing a homotopy on the obstacle positions where at one extreme of the parameter the obstacles are removed from the operating space, and at the other extreme they are located at their intended positions. We present a detailed example (with more than a thousand obstacles) from stereotactic neurosurgery with real-world data obtained from labelled MPRI scans.
△ Less
Submitted 28 November, 2022;
originally announced November 2022.
-
Accurate dark matter halo elongation from weak-lensing stacking analysis
Authors:
Elizabeth J. Gonzalez,
Kai Hoffmann,
Enrique Gaztañaga,
Diego R. García Lambas,
Pablo Fosalba,
Martin Crocce,
Francisco J. Castander,
Martín Makler
Abstract:
Shape estimates that quantify the halo anisotropic mass distribution are valuable parameters that provide information on their assembly process and evolution. Measurements of the mean shapes for a sample of cluster-sized halos can be used to test halo formation scenarios, as well as improving the modelling of potential biases in constraining cosmological parameters using these systems. In this wor…
▽ More
Shape estimates that quantify the halo anisotropic mass distribution are valuable parameters that provide information on their assembly process and evolution. Measurements of the mean shapes for a sample of cluster-sized halos can be used to test halo formation scenarios, as well as improving the modelling of potential biases in constraining cosmological parameters using these systems. In this work, we test the recovery of halo cluster shapes and masses applying weak-lensing stacking techniques. To this end, we use lensing $shear$ and a new dark matter halo catalogue, derived from the light-cone output of the cosmological simulation MICE-GC. We perform this study by combining the lensing signals obtained for several samples of halos, selected according to their mass and redshift, taking into account the main directions of the dark-matter distributions. In the analysis we test the impact of several potential systematics, such as the adopted modelling, the contribution of the neighbouring mass distributions, miscentering and misalignment effects. Our results show that, when some considerations regarding the halo relaxation state are taken into account, the lensing semi-axis ratio estimates are in agreement within $5\%$ with the mean shapes of the projected dark-matter particle distribution of the stacked halos. The methodology presented provides a useful tool to derive reliable shapes of galaxy clusters and to contrast them with those expected from numerical simulations. Furthermore, our proposed modelling, which takes into account the contribution of neighbouring halos, allows one to constrain the elongation of the surrounding mass distribution.
△ Less
Submitted 28 October, 2022; v1 submitted 19 July, 2022;
originally announced July 2022.
-
Modeling Intrinsic Galaxy Alignment in the MICE Simulation
Authors:
Kai Hoffmann,
Lucas F. Secco,
Jonathan Blazek,
Martin Crocce,
Pau Tallada-Crespí,
Simon Samuroff,
Judit Prat,
Jorge Carretero,
Pablo Fosalba,
Enrique Gaztanaga,
Francisco J. Castander
Abstract:
The intrinsic alignment (IA) of galaxies is potentially a major limitation in deriving cosmological constraints from weak lensing surveys. In order to investigate this effect we assign intrinsic shapes and orientations to galaxies in the light-cone output of the MICE simulation, spanning $\sim5000\,{\rm deg}^2$ and reaching redshift $z=1.4$. This assignment is based on a 'semi-analytic' IA model t…
▽ More
The intrinsic alignment (IA) of galaxies is potentially a major limitation in deriving cosmological constraints from weak lensing surveys. In order to investigate this effect we assign intrinsic shapes and orientations to galaxies in the light-cone output of the MICE simulation, spanning $\sim5000\,{\rm deg}^2$ and reaching redshift $z=1.4$. This assignment is based on a 'semi-analytic' IA model that uses photometric properties of galaxies as well as the spin and shape of their host halos. Advancing on previous work, we include more realistic distributions of galaxy shapes and a luminosity dependent galaxy-halo alignment. The IA model parameters are calibrated against COSMOS and BOSS LOWZ observations. The null detection of IA in observations of blue galaxies is accounted for by setting random orientations for these objects. We compare the two-point alignment statistics measured in the simulation against predictions from the analytical IA models NLA and TATT over a wide range of scales, redshifts and luminosities for red and blue galaxies separately. We find that both models fit the measurements well at scales above $8\,h^{-1}{\rm Mpc}$, while TATT outperforms NLA at smaller scales. The IA parameters derived from our fits are in broad agreement with various observational constraints from red galaxies. Lastly, we build a realistic source sample, mimicking DES Year 3 observations and use it to predict the IA contamination to the observed shear statistics. We find this prediction to be within the measurement uncertainty, which might be a consequence of the random alignment of blue galaxies in the simulation.
△ Less
Submitted 14 December, 2022; v1 submitted 28 June, 2022;
originally announced June 2022.
-
Constraining the Baryonic Feedback with Cosmic Shear Using the DES Year-3 Small-Scale Measurements
Authors:
A. Chen,
G. Aricò,
D. Huterer,
R. Angulo,
N. Weaverdyck,
O. Friedrich,
L. F. Secco,
C. Hernández-Monteagudo,
A. Alarcon,
O. Alves,
A. Amon,
F. Andrade-Oliveira,
E. Baxter,
K. Bechtol,
M. R. Becker,
G. M. Bernstein,
J. Blazek,
A. Brandao-Souza,
S. L. Bridle,
H. Camacho,
A. Campos,
A. Carnero Rosell,
M. Carrasco Kind,
R. Cawthon,
C. Chang
, et al. (117 additional authors not shown)
Abstract:
We use the small scales of the Dark Energy Survey (DES) Year-3 cosmic shear measurements, which are excluded from the DES Year-3 cosmological analysis, to constrain the baryonic feedback. To model the baryonic feedback, we adopt a baryonic correction model and use the numerical package \texttt{Baccoemu} to accelerate the evaluation of the baryonic nonlinear matter power spectrum. We design our ana…
▽ More
We use the small scales of the Dark Energy Survey (DES) Year-3 cosmic shear measurements, which are excluded from the DES Year-3 cosmological analysis, to constrain the baryonic feedback. To model the baryonic feedback, we adopt a baryonic correction model and use the numerical package \texttt{Baccoemu} to accelerate the evaluation of the baryonic nonlinear matter power spectrum. We design our analysis pipeline to focus on the constraints of the baryonic suppression effects, utilizing the implication given by a principal component analysis on the Fisher forecasts. Our constraint on the baryonic effects can then be used to better model and ameliorate the effects of baryons in producing cosmological constraints from the next generation large-scale structure surveys. We detect the baryonic suppression on the cosmic shear measurements with a $\sim 2 σ$ significance. The characteristic halo mass for which half of the gas is ejected by baryonic feedback is constrained to be $M_c > 10^{13.2} h^{-1} M_{\odot}$ (95\% C.L.). The best-fit baryonic suppression is $\sim 5\%$ at $k=1.0 {\rm Mpc}\ h^{-1}$ and $\sim 15\%$ at $k=5.0 {\rm Mpc} \ h^{-1}$. Our findings are robust with respect to the assumptions about the cosmological parameters, specifics of the baryonic model, and intrinsic alignments.
△ Less
Submitted 17 June, 2022;
originally announced June 2022.
-
New catalogue of dark-matter halo properties identified in MICE-GC -- I. Analysis of density profile distributions
Authors:
Elizabeth J. Gonzalez,
Kai Hoffmann,
Diego R. García Lambas,
Enrique Gaztañaga,
Dario Graña,
Pau Tallada-Crespí,
Jorge Carretero,
M. Victoria Santucho,
Pablo Fosalba,
Martin Crocce,
Francisco J. Castander,
Facundo Rodriguez,
Martín Makler
Abstract:
Constraints on dark matter halo masses from weak gravitational lensing can be improved significantly by using additional information about the morphology of their density distribution, leading to tighter cosmological constraints derived from the halo mass function. This work is the first of two in which we investigate the accuracy of halo morphology and mass measurements in 2D and 3D. To this end,…
▽ More
Constraints on dark matter halo masses from weak gravitational lensing can be improved significantly by using additional information about the morphology of their density distribution, leading to tighter cosmological constraints derived from the halo mass function. This work is the first of two in which we investigate the accuracy of halo morphology and mass measurements in 2D and 3D. To this end, we determine several halo physical properties in the MICE-Grand Challenge dark matter only simulation. We present a public catalogue of these properties that includes density profiles and shape parameters measured in 2D and 3D, the halo centre at the peak of the 3D density distribution as well as the gravitational and kinetic energies and angular momentum vectors. The density profiles are computed using spherical and ellipsoidal radial bins, taking into account the halo shapes. We also provide halo concentrations and masses derived from fits to 2D and 3D density profiles using NFW and Einasto models for halos with more than $1000$ particles ($\gtrsim 3 \times 10^{13} h^{-1} M_{\odot}$). We find that the Einasto model provides better fits compared to NFW, regardless of the halo relaxation state and shape. The mass and concentration parameters of the 3D density profiles derived from fits to the 2D profiles are in general biased. Similar biases are obtained when constraining mass and concentrations using a weak-lensing stacking analysis. We show that these biases depend on the radial range and density profile model adopted in the fitting procedure, but not on the halo shape.
△ Less
Submitted 13 May, 2022;
originally announced May 2022.
-
Predicting Sex and Stroke Success -- Computer-aided Player Grunt Analysis in Tennis Matches
Authors:
Lukas Stappen,
Manuel Milling,
Valentin Munst,
Korakot Hoffmann,
Bjorn W. Schuller
Abstract:
Professional athletes increasingly use automated analysis of meta- and signal data to improve their training and game performance. As in other related human-to-human research fields, signal data, in particular, contain important performance- and mood-specific indicators for automated analysis. In this paper, we introduce the novel data set SCORE! to investigate the performance of several features…
▽ More
Professional athletes increasingly use automated analysis of meta- and signal data to improve their training and game performance. As in other related human-to-human research fields, signal data, in particular, contain important performance- and mood-specific indicators for automated analysis. In this paper, we introduce the novel data set SCORE! to investigate the performance of several features and machine learning paradigms in the prediction of the sex and immediate stroke success in tennis matches, based only on vocal expression through players' grunts. The data was gathered from YouTube, labelled under the exact same definition, and the audio processed for modelling. We extract several widely used basic, expert-knowledge, and deep acoustic features of the audio samples and evaluate their effectiveness in combination with various machine learning approaches. In a binary setting, the best system, using spectrograms and a Convolutional Recurrent Neural Network, achieves an unweighted average recall (UAR) of 84.0 % for the player sex prediction task, and 60.3 % predicting stroke success, based only on acoustic cues in players' grunts of both sexes. Further, we achieve a UAR of 58.3 %, and 61.3 %, when the models are exclusively trained on female or male grunts, respectively.
△ Less
Submitted 18 February, 2022;
originally announced February 2022.
-
Dark Energy Survey Year 3 Results: Cosmology from Cosmic Shear and Robustness to Modeling Uncertainty
Authors:
L. F. Secco,
S. Samuroff,
E. Krause,
B. Jain,
J. Blazek,
M. Raveri,
A. Campos,
A. Amon,
A. Chen,
C. Doux,
A. Choi,
D. Gruen,
G. M. Bernstein,
C. Chang,
J. DeRose,
J. Myles,
A. Ferté,
P. Lemos,
D. Huterer,
J. Prat,
M. A. Troxel,
N. MacCrann,
A. R. Liddle,
T. Kacprzak,
X. Fang
, et al. (129 additional authors not shown)
Abstract:
This work and its companion paper, Amon et al. (2021), present cosmic shear measurements and cosmological constraints from over 100 million source galaxies in the Dark Energy Survey (DES) Year 3 data. We constrain the lensing amplitude parameter $S_8\equivσ_8\sqrt{Ω_\textrm{m}/0.3}$ at the 3% level in $Λ$CDM: $S_8=0.759^{+0.025}_{-0.023}$ (68% CL). Our constraint is at the 2% level when using angu…
▽ More
This work and its companion paper, Amon et al. (2021), present cosmic shear measurements and cosmological constraints from over 100 million source galaxies in the Dark Energy Survey (DES) Year 3 data. We constrain the lensing amplitude parameter $S_8\equivσ_8\sqrt{Ω_\textrm{m}/0.3}$ at the 3% level in $Λ$CDM: $S_8=0.759^{+0.025}_{-0.023}$ (68% CL). Our constraint is at the 2% level when using angular scale cuts that are optimized for the $Λ$CDM analysis: $S_8=0.772^{+0.018}_{-0.017}$ (68% CL). With cosmic shear alone, we find no statistically significant constraint on the dark energy equation-of-state parameter at our present statistical power. We carry out our analysis blind, and compare our measurement with constraints from two other contemporary weak-lensing experiments: the Kilo-Degree Survey (KiDS) and Hyper-Suprime Camera Subaru Strategic Program (HSC). We additionally quantify the agreement between our data and external constraints from the Cosmic Microwave Background (CMB). Our DES Y3 result under the assumption of $Λ$CDM is found to be in statistical agreement with Planck 2018, although favors a lower $S_8$ than the CMB-inferred value by $2.3σ$ (a $p$-value of 0.02). This paper explores the robustness of these cosmic shear results to modeling of intrinsic alignments, the matter power spectrum and baryonic physics. We additionally explore the statistical preference of our data for intrinsic alignment models of different complexity. The fiducial cosmic shear model is tested using synthetic data, and we report no biases greater than 0.3$σ$ in the plane of $S_8\timesΩ_\textrm{m}$ caused by uncertainties in the theoretical models.
△ Less
Submitted 13 January, 2022; v1 submitted 27 May, 2021;
originally announced May 2021.
-
Redshift and stellar mass dependence of intrinsic shapes of disc-dominated galaxies from COSMOS observations below $z = 1.0$
Authors:
Kai Hoffmann,
Clotilde Laigle,
Nora Elisa Chisari,
Pau Tallada,
Romain Teyssier,
Yohan Dubois,
Julien Devriendt
Abstract:
The high abundance of disc galaxies without a large central bulge challenges predictions of current hydrodynamic simulations of galaxy formation. We aim to shed light on the formation of these objects by studying the redshift and mass dependence of their intrinsic 3D shape distributions in the COSMOS galaxy survey below redshift $z=1.0$. This distribution is inferred from the observed distribution…
▽ More
The high abundance of disc galaxies without a large central bulge challenges predictions of current hydrodynamic simulations of galaxy formation. We aim to shed light on the formation of these objects by studying the redshift and mass dependence of their intrinsic 3D shape distributions in the COSMOS galaxy survey below redshift $z=1.0$. This distribution is inferred from the observed distribution of 2D shapes, using a reconstruction method which we test using hydrodynamic simulations. Our tests reveal a moderate bias for the inferred average disc circularity and relative thickness, but a large bias on the dispersion of these quantities. Applying the reconstruction method on COSMOS data, we find variations of the average disc circularity and relative thickness with redshift of around $\sim1\%$ and $\sim10\%$ respectively, which is comparable to the error estimates on these quantities. The average relative disc thickness shows a significant mass dependence which can be accounted for by the scaling of disc radius with galaxy mass. We conclude that our data provides no evidence for a strong dependence of the average circularity and absolute thickness of disc-dominated galaxies on redshift and mass that is significant with respect to the statistical uncertainties in our analysis. These findings are expected in the absence of disruptive merging or feedback events that would affect galaxy shapes. They hence support a scenario where present-day discs form early ($z>1.0$) and subsequently undergo a tranquil evolution in isolation. However, more data and a better understanding of systematics are needed to reaffirm our results.
△ Less
Submitted 14 September, 2022; v1 submitted 26 October, 2020;
originally announced October 2020.
-
A robustness measure for singular point and index estimation in discretized orientation and vector fields
Authors:
Karl B. Hoffmann,
Ivo F. Sbalzarini
Abstract:
The identification of singular points or topological defects in discretized vector fields occurs in diverse areas ranging from the polarization of the cosmic microwave background to liquid crystals to fingerprint recognition and bio-medical imaging. Due to their discrete nature, defects and their topological charge cannot depend continuously on each single vector, but they discontinuously change a…
▽ More
The identification of singular points or topological defects in discretized vector fields occurs in diverse areas ranging from the polarization of the cosmic microwave background to liquid crystals to fingerprint recognition and bio-medical imaging. Due to their discrete nature, defects and their topological charge cannot depend continuously on each single vector, but they discontinuously change as soon as a vector changes by more than a threshold. Considering this threshold of admissible change at the level of vectors, we develop a robustness measure for discrete defect estimators. Here, we compare different template paths for defect estimation in discretized vector or orientation fields. Sampling prototypical vector field patterns around defects shows that the robustness increases with the length of template path, but less so in the presence of noise on the vectors. We therefore find an optimal trade-off between resolution and robustness against noise for relatively small templates, except for the "single pixel" defect analysis, which cannot exclude zero robustness. The presented robustness measure paves the way for uncertainty quantification of defects in discretized vector fields.
△ Less
Submitted 30 September, 2020;
originally announced September 2020.
-
The Three Hundred project: shapes and radial alignment of satellite, infalling, and backsplash galaxies
Authors:
Alexander Knebe,
Matias Gamez-Marin,
Frazer R. Pearce,
Weiguang Cui,
Kai Hoffmann,
Marco De Petris,
Chris Power,
Roan Haggar,
Robert Mostoghiu
Abstract:
Using 324 numerically modelled galaxy clusters we investigate the radial and galaxy-halo alignment of dark matter subhaloes and satellite galaxies orbiting within and around them. We find that radial alignment depends on distance to the centre of the galaxy cluster but appears independent of the dynamical state of the central host cluster. Furthermore, we cannot find a relation between radial alig…
▽ More
Using 324 numerically modelled galaxy clusters we investigate the radial and galaxy-halo alignment of dark matter subhaloes and satellite galaxies orbiting within and around them. We find that radial alignment depends on distance to the centre of the galaxy cluster but appears independent of the dynamical state of the central host cluster. Furthermore, we cannot find a relation between radial alignment of the halo or galaxy shape with its own mass. We report that backsplash galaxies, i.e. objects that have already passed through the cluster radius but are now located in the outskirts, show a stronger radial alignment than infalling objects. We further find that there exists a population of well radially aligned objects passing very close to the central cluster's centre which were found to be on highly radial orbit.
△ Less
Submitted 20 May, 2020;
originally announced May 2020.
-
A Multiple Power Law Distribution for Initial Mass Functions
Authors:
Christopher Essex,
Shantanu Basu,
Janett Prehl,
Karl Heinz Hoffmann
Abstract:
We introduce a new multi-power-law distribution for the Initial Mass Function (IMF) to explore its potential properties. It follows on prior work that introduced mechanisms accounting for mass accretion in star formation, developed within the framework of general evolution equations for the mass distribution of accreting and non-accreting (proto)stars. This paper uses the same fundamental framewor…
▽ More
We introduce a new multi-power-law distribution for the Initial Mass Function (IMF) to explore its potential properties. It follows on prior work that introduced mechanisms accounting for mass accretion in star formation, developed within the framework of general evolution equations for the mass distribution of accreting and non-accreting (proto)stars. This paper uses the same fundamental framework to demonstrate that the interplay between a mass-dependent and a time-dependent step-like dropout rate from accretion leads to IMFs that exhibit multiple power laws for an exponential mass growth. While the mass-dependent accretion and its dropout is intrinsic to each star, the time-dependent dropout might be tied to a specific history such as the rapid consumption of nebular material by nearby stars or the sweeping away of some material by shock waves. The time-dependent dropout folded into the mass-dependent process of star formation is shown to have a significant influence on the IMFs.
△ Less
Submitted 23 March, 2020;
originally announced March 2020.
-
Second sound in the BEC-BCS crossover
Authors:
Daniel Kai Hoffmann,
Vijay Pal Singh,
Thomas Paintner,
Manuel Jäger,
Wolfgang Limmer,
Ludwig Mathey,
Johannes Hecker Denschlag
Abstract:
Second sound is an entropy wave which propagates in the superfluid component of a quantum liquid. Because it is an entropy wave, it probes the thermodynamic properties of the quantum liquid which are determined, e.g., by the interaction strength between the particles of the quantum liquid and their temperature. Here, we study second sound propagation for a large range of interaction strengths with…
▽ More
Second sound is an entropy wave which propagates in the superfluid component of a quantum liquid. Because it is an entropy wave, it probes the thermodynamic properties of the quantum liquid which are determined, e.g., by the interaction strength between the particles of the quantum liquid and their temperature. Here, we study second sound propagation for a large range of interaction strengths within the crossover between a Bose-Einstein condensate (BEC) and the Bardeen-Cooper-Schrieffer (BCS) superfluid. In particular, we investigate the strongly-interacting regime where currently theoretical predictions only exist in terms of an interpolation between the BEC, BCS and unitary regimes. Working with a quantum gas of ultracold fermionic $^6$Li atoms with tunable interactions, we show that the second sound speed varies only slightly in the crossover regime. We gain deeper insights into sound propagation and excitation of second sound by varying the excitation procedure which ranges from a sudden force pulse to a gentle heating pulse at the cloud center. These measurements are accompanied by classical-field simulations which help with the interpretation of the experimental data. Furthermore, we determine the spatial extension of the superfluid phase and estimate the superfluid density. In the future, this may be used to construct the so far unknown equation of state throughout the crossover.
△ Less
Submitted 15 March, 2020;
originally announced March 2020.
-
A hybrid echocardiography-computational fluid dynamics framework for ventricular flow simulations
Authors:
Mohammadali Hedayat,
Tatsat R. Patel,
Marek Belohlavek,
Kenneth R. Hoffmann,
Iman Borazjani
Abstract:
Image-based computational fluid dynamics (CFD) has emerged as a powerful tool to study cardiovascular flows while 2D echocardiography (echo) is the most widely used non-invasive imaging modality for diagnosis of heart disease. Here, echo is combined with CFD, i.e., an echo-CFD framework, to study ventricular flows. To achieve this, our previous 3D reconstruction from multiple 2D-echo at standard c…
▽ More
Image-based computational fluid dynamics (CFD) has emerged as a powerful tool to study cardiovascular flows while 2D echocardiography (echo) is the most widely used non-invasive imaging modality for diagnosis of heart disease. Here, echo is combined with CFD, i.e., an echo-CFD framework, to study ventricular flows. To achieve this, our previous 3D reconstruction from multiple 2D-echo at standard cross-sections is extended by 1) reconstructing valves (aortic and mitral valves) from 2D-echo images and the superior wall; 2) optimizing the location and orientation of the surfaces to incorporate the physiological assumption of fixed apex as a reference (fixed) point for reconstruction; and 3) incorporating several smoothing algorithms to remove the nonphysical oscillations (ringing) near the basal section observed in our previous study. The main parameters of the reconstructed left ventricle (LV) are all within the physiologic range. Our results show that the mitral valve can significantly change the flow pattern inside the LV during the diastole phase in terms of ring vortex formation/propagation as well as the flow circulation. In addition, the abnormal shape of unhealthy LV can drastically change the flow field during the diastole phase. Furthermore, the hemodynamic energy loss, as an indicator of the LV pumping performance, for different test cases is calculated which shows a larger energy loss for unhealthy LV compared to the healthy one.
△ Less
Submitted 19 September, 2019;
originally announced September 2019.
-
Multivariate Newton Interpolation
Authors:
Michael Hecht,
Karl B. Hoffmann,
Bevan L. Cheeseman,
Ivo F. Sbalzarini
Abstract:
For $m,n \in \mathbb{N}$, $m\geq 1$ and a given function $f : \mathbb{R}^m\longrightarrow \mathbb{R}$, the polynomial interpolation problem (PIP) is to determine a unisolvent node set $P_{m,n} \subseteq \mathbb{R}^m$ of $N(m,n):=|P_{m,n}|=\binom{m+n}{n}$ points and the uniquely defined polynomial $Q_{m,n,f}\in Π_{m,n}$ in $m$ variables of degree $\mathrm{deg}(Q_{m,n,f})\leq n \in \mathbb{N}$ that…
▽ More
For $m,n \in \mathbb{N}$, $m\geq 1$ and a given function $f : \mathbb{R}^m\longrightarrow \mathbb{R}$, the polynomial interpolation problem (PIP) is to determine a unisolvent node set $P_{m,n} \subseteq \mathbb{R}^m$ of $N(m,n):=|P_{m,n}|=\binom{m+n}{n}$ points and the uniquely defined polynomial $Q_{m,n,f}\in Π_{m,n}$ in $m$ variables of degree $\mathrm{deg}(Q_{m,n,f})\leq n \in \mathbb{N}$ that fits $f$ on $P_{m,n}$, i.e., $Q_{m,n,f}(p) = f(p)$, $\forall\, p \in P_{m,n}$. For $m=1$ the solution to the PIP is well known. In higher dimensions, however, no closed framework was available. We here present a generalization of the classic Newton interpolation from one-dimensional to arbitrary-dimensional spaces. Further we formulate an algorithm, termed PIP-SOLVER, based on a multivariate divided difference scheme that computes the solution $Q_{m,n,f}$ in $\mathcal{O}\big(N(m,n)^2\big)$ time using $\mathcal{O}\big(mN(m,n)\big)$ memory. Further, we introduce unisolvent Newton-Chebyshev nodes and show that these nodes avoid Runge's phenomenon in the sense that arbitrary periodic Sobolev functions $f \in H^k(Ω,\mathbb{R}) \subsetneq C^0(Ω,\mathbb{R})$, $Ω=[-1,1]^m$ of regularity $k >m/2$ can be uniformly approximated, i.e., $ \lim_{n\rightarrow \infty}||\,f -Q_{m,n,f} \,||_{C^0(Ω)}= 0$. Numerical experiments demonstrate the computational performance and approximation accuracy of the PIP-SOLVER in practice. We expect the presented results to be relevant for many applications, including numerical solvers, quadrature, non-linear optimization, polynomial regression, adaptive sampling, Bayesian inference, and spectral analysis.
△ Less
Submitted 19 March, 2020; v1 submitted 11 December, 2018;
originally announced December 2018.
-
Reaction kinetics of ultracold molecule-molecule collisions
Authors:
Daniel K. Hoffmann,
Thomas Paintner,
Wolfgang Limmer,
Dmitry S. Petrov,
Johannes Hecker Denschlag
Abstract:
Studying chemical reactions on a state-to-state level tests and improves our fundamental understanding of chemical processes. For such investigations it is convenient to make use of ultracold atomic and molecular reactants as they can be prepared in well defined internal and external quantum states$^{1-4}$. In general, even cold reactions have many possible final product states$^{5-15}$ and reacti…
▽ More
Studying chemical reactions on a state-to-state level tests and improves our fundamental understanding of chemical processes. For such investigations it is convenient to make use of ultracold atomic and molecular reactants as they can be prepared in well defined internal and external quantum states$^{1-4}$. In general, even cold reactions have many possible final product states$^{5-15}$ and reaction channels are therefore hard to track individually$^{16}$. In special cases, however, only a single reaction channel is essentially participating, as observed e.g. in the recombination of two atoms forming a Feshbach molecule$^{17-19}$ or in atom-Feshbach molecule exchange reactions$^{20,21}$. Here, we investigate a single-channel reaction of two Li$_2$-Feshbach molecules where one of the molecules dissociates into two atoms $2\mathrm{AB}\Rightarrow \mathrm{AB}+\mathrm{A}+\mathrm{B}$. The process is a prototype for a class of four-body collisions where two reactants produce three product particles. We measure the collisional dissociation rate constant of this process as a function of collision energy/ temperature and scattering length. We confirm an Arrhenius-law dependence on the collision energy, an $a^4$ power-law dependence on the scattering length $a$ and determine a universal four body reaction constant.
△ Less
Submitted 21 May, 2018;
originally announced May 2018.
-
A Dual Power Law Distribution for the Stellar Initial Mass Function
Authors:
Karl Heinz Hoffmann,
Christopher Essex,
Shantanu Basu,
Janett Prehl
Abstract:
We introduce a new dual power law (DPL) probability distribution function for the mass distribution of stellar and substellar objects at birth, otherwise known as the initial mass function (IMF). The model contains both deterministic and stochastic elements, and provides a unified framework within which to view the formation of brown dwarfs and stars resulting from an accretion process that starts…
▽ More
We introduce a new dual power law (DPL) probability distribution function for the mass distribution of stellar and substellar objects at birth, otherwise known as the initial mass function (IMF). The model contains both deterministic and stochastic elements, and provides a unified framework within which to view the formation of brown dwarfs and stars resulting from an accretion process that starts from extremely low mass seeds. It does not depend upon a top down scenario of collapsing (Jeans) masses or an initial lognormal or otherwise IMF-like distribution of seed masses. Like the modified lognormal power law (MLP) distribution, the DPL distribution has a power law at the high mass end, as a result of exponential growth of mass coupled with equally likely stopping of accretion at any time interval. Unlike the MLP, a power law decay also appears at the low mass end of the IMF. This feature is closely connected to the accretion stopping probability rising from an initially low value up to a high value. This might be associated with physical effects of ejections sometimes (i.e., rarely) stopping accretion at early times followed by outflow driven accretion stopping at later times, with the transition happening at a critical time (therefore mass). Comparing the DPL to empirical data, the critical mass is close to the substellar mass limit, suggesting that the onset of nuclear fusion plays an important role in the subsequent accretion history of a young stellar object.
△ Less
Submitted 14 May, 2018;
originally announced May 2018.
-
Pair fraction in a finite temperature Fermi gas on the BEC side of the BCS-BEC crossover
Authors:
Thomas Paintner,
Daniel K. Hoffmann,
Manuel Jäger,
Wolfgang Limmer,
Wladimir Schoch,
Benjamin Deissler,
Michele Pini,
Pierbiagio Pieri,
Giancarlo Calvanese Strinati,
Cheng Chin,
Johannes Hecker Denschlag
Abstract:
We investigate pairing in a strongly interacting two-component Fermi gas with positive scattering length. In this regime, pairing occurs at temperatures above the superfluid critical temperature; unbound fermions and pairs coexist in thermal equilibrium. Measuring the total number of these fermion pairs in the gas we systematically investigate the phases in the sectors of pseudogap and preformed-p…
▽ More
We investigate pairing in a strongly interacting two-component Fermi gas with positive scattering length. In this regime, pairing occurs at temperatures above the superfluid critical temperature; unbound fermions and pairs coexist in thermal equilibrium. Measuring the total number of these fermion pairs in the gas we systematically investigate the phases in the sectors of pseudogap and preformed-pair. Our measurements quantitatively test predictions from two theoretical models. Interestingly, we find that already a model based on classical atom-molecule equilibrium describes our data quite well.
△ Less
Submitted 30 April, 2019; v1 submitted 28 March, 2018;
originally announced March 2018.
-
Signatures of Cosmic Reionization on the 21cm 2- and 3-point Correlation Function I: Quadratic Bias Modeling
Authors:
Kai Hoffmann,
Yi Mao,
Jiachuan Xu,
Houjun Mo,
Benjamin D. Wandelt
Abstract:
The three-point correlation function (3PCF) of the 21cm brightness temperature from the Epoch of Reionization (EoR) probes complementary information to the commonly studied two-point correlation function (2PCF) about the morphology of ionized regions. We investigate the 21cm 2PCF and 3PCF in configuration space using semi-numerical simulations and test whether they can be described by the local qu…
▽ More
The three-point correlation function (3PCF) of the 21cm brightness temperature from the Epoch of Reionization (EoR) probes complementary information to the commonly studied two-point correlation function (2PCF) about the morphology of ionized regions. We investigate the 21cm 2PCF and 3PCF in configuration space using semi-numerical simulations and test whether they can be described by the local quadratic bias model. We find that fits of bias model predictions for the 2PCF and 3PCF deviate from our measurements by $\sim 20\%$ at scales above the typical size of ionized regions ($\simeq 30$ Mpc) and at early times with global neutral fractions of $\langle x_{\rm HI} \rangle \gtrsim 0.7$. At later times and smaller scales these deviations increase strongly, indicating a break down of the bias model. The 2PCF and 3PCF fits of the linear bias parameter agree at the $10\%$ level for different EoR model configurations. This agreement holds, when adding redshift space distortions to the simulations. The relation between spatial fluctuations in the matter density and the 21cm signal, as predicted by the bias model, is consistent with direct measurements of this relation in simulations for large smoothing scales ($\gtrsim 30$ Mpc). From this latter test we conclude that negative amplitudes of the 21cm 3PCF result from negative bias parameters, which describe the anti-correlation between the matter over-densities and the 21cm signal during the EoR. However, a more detailed interpretation of the bias parameters may require a description of non-local contributions to the bias model.
△ Less
Submitted 31 May, 2019; v1 submitted 7 February, 2018;
originally announced February 2018.
-
A Quadratic-Time Algorithm for General Multivariate Polynomial Interpolation
Authors:
M. Hecht,
B. L. Cheeseman,
K. B. Hoffmann,
I. F. Sbalzarini
Abstract:
For $m,n \in \mathbb{N}$, $m\geq 1$ and a given function $f : \mathbb{R}^m\longrightarrow \mathbb{R}$ the polynomial interpolation problem (PIP) is to determine a \emph{generic node set} $P \subseteq \mathbb{R}^m$ and the coefficients of the uniquely defined polynomial $Q\in\mathbb{R}[x_1,\dots,x_m]$ in $m$ variables of degree $\mathrm{deg}(Q)\leq n \in \mathbb{N}$ that fits $f$ on $P$, i.e.,…
▽ More
For $m,n \in \mathbb{N}$, $m\geq 1$ and a given function $f : \mathbb{R}^m\longrightarrow \mathbb{R}$ the polynomial interpolation problem (PIP) is to determine a \emph{generic node set} $P \subseteq \mathbb{R}^m$ and the coefficients of the uniquely defined polynomial $Q\in\mathbb{R}[x_1,\dots,x_m]$ in $m$ variables of degree $\mathrm{deg}(Q)\leq n \in \mathbb{N}$ that fits $f$ on $P$, i.e., $Q(p) = f(p)$, $\forall\, p \in P$. We here show that in general, i.e., for arbitrary $m,n \in \mathbb{N}$, $m \geq 1$, there exists an algorithm that determines $P$ and computes the $N(\mbox{m,n})=\#P$ coefficients of $Q$ in $\mathcal{O}\big(N(\mbox{m,n})^2\big)$ time using $\mathcal{O}\big(\mbox{m}N(\mbox{m,n})\big)$ storage, without inverting the occurring Vandermonde matrix. We provide such an algorithm, termed PIP-SOLVER, based on a recursive decomposition of the problem and prove its correctness. Since the present approach solves the PIP without matrix inversion, it is computationally more efficient and numerically more robust than previous approaches. We demonstrate this in numerical experiments and compare with previous approaches based on matrix inversion and linear systems solving.
△ Less
Submitted 30 October, 2017;
originally announced October 2017.
-
A Dynamically Diluted Alignment Model Reveals the Impact of Cell Turnover on the Plasticity of Tissue Polarity Patterns
Authors:
Karl B. Hoffmann,
Anja Voss-Böhme,
Jochen C. Rink,
Lutz Brusch
Abstract:
The polarisation of cells and tissues is fundamental for tissue morphogenesis during biological development and regeneration. A deeper understanding of biological polarity pattern formation can be gained from the consideration of pattern reorganisation in response to an opposing instructive cue, which we here consider by example of experimentally inducible body axis inversions in planarian flatwor…
▽ More
The polarisation of cells and tissues is fundamental for tissue morphogenesis during biological development and regeneration. A deeper understanding of biological polarity pattern formation can be gained from the consideration of pattern reorganisation in response to an opposing instructive cue, which we here consider by example of experimentally inducible body axis inversions in planarian flatworms. Our dynamically diluted alignment model represents three processes: entrainment of cell polarity by a global signal, local cell-cell coupling aligning polarity among neighbours and cell turnover inserting initially unpolarised cells. We show that a persistent global orienting signal determines the final mean polarity orientation in this stochastic model. Combining numerical and analytical approaches, we find that neighbour coupling retards polarity pattern reorganisation, whereas cell turnover accelerates it. We derive a formula for an effective neighbour coupling strength integrating both effects and find that the time of polarity reorganisation depends linearly on this effective parameter and no abrupt transitions are observed. This allows to determine neighbour coupling strengths from experimental observations. Our model is related to a dynamic $8$-Potts model with annealed site-dilution and makes testable predictions regarding the polarisation of dynamic systems, such as the planarian epithelium.
△ Less
Submitted 9 October, 2017;
originally announced October 2017.
-
Testing the consistency of three-point halo clustering in Fourier and configuration space
Authors:
Kai Hoffmann,
Enrique Gaztanaga,
Roman Scoccimarro,
Martin Crocce
Abstract:
We compare reduced three-point correlations $Q$ of matter, haloes (as proxies for galaxies) and their cross correlations, measured in a total simulated volume of $\sim100 \ (h^{-1} \text{Gpc})^{3}$, to predictions from leading order perturbation theory on a large range of scales in configuration space. Predictions for haloes are based on the non-local bias model, employing linear ($b_1$) and non-l…
▽ More
We compare reduced three-point correlations $Q$ of matter, haloes (as proxies for galaxies) and their cross correlations, measured in a total simulated volume of $\sim100 \ (h^{-1} \text{Gpc})^{3}$, to predictions from leading order perturbation theory on a large range of scales in configuration space. Predictions for haloes are based on the non-local bias model, employing linear ($b_1$) and non-linear ($c_2$, $g_2$) bias parameters, which have been constrained previously from the bispectrum in Fourier space. We also study predictions from two other bias models, one local ($g_2=0$) and one in which $c_2$ and $g_2$ are determined by $b_1$ via an approximately universal relation. Overall, measurements and predictions agree when $Q$ is derived for triangles with $(r_1r_2r_3)^{1/3} \gtrsim 60 h^{-1}\text{Mpc}$, where $r_{1-3}$ are the sizes of the triangle legs. Predictions for $Q_{matter}$, based on the linear power spectrum, show significant deviations from the measurements at the BAO scale (given our small measurement errors), which strongly decrease when adding a damping term or using the non-linear power spectrum, as expected. Predictions for $Q_{halo}$ agree best with measurements at large scales when considering non-local contributions. The universal bias model works well for haloes and might therefore be also useful for tightening constraints on $b_1$ from $Q$ in galaxy surveys. Such constraints are independent of the amplitude of matter density fluctuation ($σ_8$) and hence break the degeneracy between $b_1$ and $σ_8$, present in galaxy two-point correlations.
△ Less
Submitted 4 February, 2018; v1 submitted 29 August, 2017;
originally announced August 2017.
-
On the propagation of singularities for a class of linearised hybrid inverse problems
Authors:
Guillaume Bal,
Kristoffer Hoffmann,
Kim Knudsen
Abstract:
For a general formulation of linearised hybrid inverse problems in impedance tomography, the qualitative properties of the solutions are analysed. Using an appropriate scalar pseudo-differential formulation, the problems are shown to permit propagating singularities under certain non-elliptic conditions, and the associated directions of propagation are precisely identified relative to the directio…
▽ More
For a general formulation of linearised hybrid inverse problems in impedance tomography, the qualitative properties of the solutions are analysed. Using an appropriate scalar pseudo-differential formulation, the problems are shown to permit propagating singularities under certain non-elliptic conditions, and the associated directions of propagation are precisely identified relative to the directions in which ellipticity is lost. The same result is found in the setting for the corresponding normal formulation of the scalar pseudo-differential equations. A numerical reconstruction procedure based of the least squares finite element method is derived, and a series of numerical experiments visualise exactly how the loss of ellipticity manifests itself as propagating singularities.
△ Less
Submitted 2 August, 2017;
originally announced August 2017.
-
Spin-transport, spin-torque and memory in antiferromagnetic devices: Part of a collection of reviews on antiferromagnetic spintronics
Authors:
J. Železný,
P. Wadley,
K. Olejník. A. Hoffmann,
H. Ohno
Abstract:
Ferromagnets are key materials for sensing and memory applications. In contrast, antiferromagnets that represent the more common form of magnetically ordered materials, have so far found less practical application beyond their use for establishing reference magnetic orientations via exchange bias. This might change in the future due to the recent progress in materials research and discoveries of a…
▽ More
Ferromagnets are key materials for sensing and memory applications. In contrast, antiferromagnets that represent the more common form of magnetically ordered materials, have so far found less practical application beyond their use for establishing reference magnetic orientations via exchange bias. This might change in the future due to the recent progress in materials research and discoveries of antiferromagnetic spintronic phenomena suitable for device applications. Experimental demonstrations of the electrical switching and electrical detection of the Néel order open a route towards memory devices based on antiferromagnets. Apart from the radiation and magnetic-field hardness, memory cells fabricated in antiferromagnets are inherently multilevel which could be used for neuromorphic computing. Switching speeds attainable in antiferromagnets far exceed those of the ferromagnetic and semiconductor memory technologies. Here we review the recent progress in electronic spin-transport and spin-torque phenomena in antiferromagnets that are dominantly of the relativistic quantum mechanics origin. We discuss their utility in pure antiferromagnetic or hybrid ferromagnetic/antiferromagnetic memory devices
△ Less
Submitted 30 May, 2017;
originally announced May 2017.
-
Linear and non-linear bias: predictions vs. measurements
Authors:
Kai Hoffmann,
Julien Bel,
Enrique Gaztanaga
Abstract:
We study the linear and non-linear bias parameters which determine the mapping between the distributions of galaxies and the full matter density fields, comparing different measurements and predictions. Associating galaxies with dark matter haloes in the MICE Grand Challenge N-body simulation we directly measure the bias parameters by comparing the smoothed density fluctuations of haloes and matte…
▽ More
We study the linear and non-linear bias parameters which determine the mapping between the distributions of galaxies and the full matter density fields, comparing different measurements and predictions. Associating galaxies with dark matter haloes in the MICE Grand Challenge N-body simulation we directly measure the bias parameters by comparing the smoothed density fluctuations of haloes and matter in the same region at different positions as a function of smoothing scale. Alternatively we measure the bias parameters by matching the probability distributions of halo and matter density fluctuations, which can be applied to observations. These direct bias measurements are compared to corresponding measurements from two-point and different third-order correlations, as well as predictions from the peak-background model, which we presented in previous articles using the same data. We find an overall variation of the linear bias measurements and predictions of $\sim 5 \%$ with respect to results from two-point correlations for different halo samples with masses between $\sim 10^{12} - 10^{15}$ $h^{-1}M_\odot$ at the redshifts $z=0.0$ and $0.5$. Variations between the second- and third-order bias parameters from the different methods show larger variations, but with consistent trends in mass and redshift. The various bias measurements reveal a tight relation between the linear and the quadratic bias parameters, which is consistent with results from the literature based on simulations with different cosmologies. Such a universal relation might improve constraints on cosmological models, derived from second-order clustering statistics at small scales or higher-order clustering statistics.
△ Less
Submitted 31 October, 2016; v1 submitted 4 July, 2016;
originally announced July 2016.
-
Holographic method for site-resolved detection of a 2D array of ultracold atoms
Authors:
Daniel Kai Hoffmann,
Benjamin Deissler,
Wolfgang Limmer,
Johannes Hecker Denschlag
Abstract:
We propose a novel approach to site-resolved detection of a 2D gas of ultracold atoms in an optical lattice. A near resonant laser beam is coherently scattered by the atomic array and its interference pattern is holographically recorded by superimposing it with a reference laser beam on a CCD chip. Fourier transformation of the recorded intensity pattern reconstructs the atomic distribution in the…
▽ More
We propose a novel approach to site-resolved detection of a 2D gas of ultracold atoms in an optical lattice. A near resonant laser beam is coherently scattered by the atomic array and its interference pattern is holographically recorded by superimposing it with a reference laser beam on a CCD chip. Fourier transformation of the recorded intensity pattern reconstructs the atomic distribution in the lattice with single-site resolution. The holographic detection method requires only a few hundred scattered photons per atom in order to achieve a high reconstruction fidelity. Therefore, additional cooling during detection might not be necessary even for light atomic elements such as lithium.
△ Less
Submitted 10 May, 2016;
originally announced May 2016.
-
What determines large scale galaxy clustering: halo mass or local density?
Authors:
Arnau Pujol,
Kai Hoffmann,
Noelia Jiménez,
Enrique Gaztañaga
Abstract:
Using dark matter simulations we show how halo bias is determined by local density and not by halo mass. This is not totally surprising, as according to the peak-background split model, local density is the property that constraints bias at large scales. Massive haloes have a high clustering because they reside in high density regions. Small haloes can be found in a wide range of environments whic…
▽ More
Using dark matter simulations we show how halo bias is determined by local density and not by halo mass. This is not totally surprising, as according to the peak-background split model, local density is the property that constraints bias at large scales. Massive haloes have a high clustering because they reside in high density regions. Small haloes can be found in a wide range of environments which determine their clustering amplitudes differently. This contradicts the assumption of standard Halo Occupation Distribution (HOD) models that the bias and occupation of haloes is determined solely by their mass. We show that the bias of central galaxies from semi-analytic models of galaxy formation as a function of luminosity and colour is not correctly predicted by the standard HOD model. Using local density instead of halo mass the HOD model correctly predicts galaxy bias. These results indicate the need to include information about local density and not only mass in order to correctly apply HOD analysis in these galaxy samples. This new model can be readily applied to observations and has the advantage that the galaxy density can be directly observed, in contrast with the dark matter halo mass.
△ Less
Submitted 15 September, 2016; v1 submitted 6 October, 2015;
originally announced October 2015.
-
Trapping in irradiated p-on-n silicon sensors at fluences anticipated at the HL-LHC outer tracker
Authors:
W. Adam,
T. Bergauer,
M. Dragicevic,
M. Friedl,
R. Fruehwirth,
M. Hoch,
J. Hrubec,
M. Krammer,
W. Treberspurg,
W. Waltenberger,
S. Alderweireldt,
W. Beaumont,
X. Janssen,
S. Luyckx,
P. Van Mechelen,
N. Van Remortel,
A. Van Spilbeeck,
P. Barria,
C. Caillol,
B. Clerbaux,
G. De Lentdecker,
D. Dobur,
L. Favart,
A. Grebenyuk,
Th. Lenzi
, et al. (663 additional authors not shown)
Abstract:
The degradation of signal in silicon sensors is studied under conditions expected at the CERN High-Luminosity LHC. 200 $μ$m thick n-type silicon sensors are irradiated with protons of different energies to fluences of up to $3 \cdot 10^{15}$ neq/cm$^2$. Pulsed red laser light with a wavelength of 672 nm is used to generate electron-hole pairs in the sensors. The induced signals are used to determi…
▽ More
The degradation of signal in silicon sensors is studied under conditions expected at the CERN High-Luminosity LHC. 200 $μ$m thick n-type silicon sensors are irradiated with protons of different energies to fluences of up to $3 \cdot 10^{15}$ neq/cm$^2$. Pulsed red laser light with a wavelength of 672 nm is used to generate electron-hole pairs in the sensors. The induced signals are used to determine the charge collection efficiencies separately for electrons and holes drifting through the sensor. The effective trapping rates are extracted by comparing the results to simulation. The electric field is simulated using Synopsys device simulation assuming two effective defects. The generation and drift of charge carriers are simulated in an independent simulation based on PixelAV. The effective trapping rates are determined from the measured charge collection efficiencies and the simulated and measured time-resolved current pulses are compared. The effective trapping rates determined for both electrons and holes are about 50% smaller than those obtained using standard extrapolations of studies at low fluences and suggests an improved tracker performance over initial expectations.
△ Less
Submitted 7 May, 2015;
originally announced May 2015.
-
Non-local bias contribution to third-order galaxy correlations
Authors:
Julien Bel,
Kai Hoffmann,
Enrique Gaztañaga
Abstract:
We study halo clustering bias with second- and third-order statistics of halo and matter density fields in the MICE Grand Challenge simulation. We verify that two-point correlations deliver reliable estimates of the linear bias parameters at large scales, while estimations from the variance can be significantly affected by non-linear and possibly non-local contributions to the bias function. Combi…
▽ More
We study halo clustering bias with second- and third-order statistics of halo and matter density fields in the MICE Grand Challenge simulation. We verify that two-point correlations deliver reliable estimates of the linear bias parameters at large scales, while estimations from the variance can be significantly affected by non-linear and possibly non-local contributions to the bias function. Combining three-point auto- and cross-correlations we find, for the first time in configuration space, evidence for the presence of such non-local contributions. These contributions are consistent with predicted second-order non-local effects on the bias functions originating from the dark matter tidal field. Samples of massive haloes show indications of bias (local or non-local) beyond second order. Ignoring non-local bias causes $20-30$\% and $5-10$\% overestimation of the linear bias from three-point auto- and cross-correlations respectively. We study two third-order bias estimators which are not affected by second-order non-local contributions. One is a combination of three-point auto- and cross- correlation. The other is a combination of third-order one- and two-point cumulants. Both methods deliver accurate bias estimations of the linear bias. Furthermore their estimations of second-order bias agree mutually. Ignoring non-local bias causes higher values of the second-order bias from three-point correlations. Our results demonstrate that third-order statistics can be employed for breaking the growth-bias degeneracy.
△ Less
Submitted 16 October, 2015; v1 submitted 8 April, 2015;
originally announced April 2015.
-
Comparing halo bias from abundance and clustering
Authors:
Kai Hoffmann,
Julien Bel,
Enrique Gaztanaga
Abstract:
We model the abundance of haloes in the $\sim(3 \ \text{Gpc}/h)^3$ volume of the MICE Grand Challenge simulation by fitting the universal mass function with an improved Jack-Knife error covariance estimator that matches theory predictions. We present unifying relations between different fitting models and new predictions for linear ($b_1$) and non-linear ($c_2$ and $c_3$) halo clustering bias. Dif…
▽ More
We model the abundance of haloes in the $\sim(3 \ \text{Gpc}/h)^3$ volume of the MICE Grand Challenge simulation by fitting the universal mass function with an improved Jack-Knife error covariance estimator that matches theory predictions. We present unifying relations between different fitting models and new predictions for linear ($b_1$) and non-linear ($c_2$ and $c_3$) halo clustering bias. Different mass function fits show strong variations in their performance when including the low mass range ($M_h \lesssim 3 \ 10^{12} \ M_{\odot}/h$) in the analysis. Together with fits from the literature we find an overall variation in the amplitudes of around $10$% in the low mass and up to $50$% in the high mass (galaxy cluster) range ($M_h > 10^{14} \ M_{\odot}/h$). These variations propagate into a $10$% change in $b_1$ predictions and a $50$% change in $c_2$ or $c_3$. Despite these strong variations we find universal relations between $b_1$ and $c_2$ or $c_3$ for which we provide simple fits. Excluding low mass haloes, different models fitted with reasonable goodness in this analysis, show percent level agreement in their $b_1$ predictions, but are systematically $5-10$% lower than the bias directly measured with two-point halo-mass clustering. This result confirms previous findings derived from smaller volumes (and smaller masses). Inaccuracies in the bias predictions lead to $5-10$% errors in growth measurements. They also affect any HOD fitting or (cluster) mass calibration from clustering measurements.
△ Less
Submitted 14 April, 2015; v1 submitted 1 March, 2015;
originally announced March 2015.
-
Non-local kinetic and macroscopic models for self-organised animal aggregations
Authors:
José A. Carrillo,
Raluca Eftimie,
Franca K. O. Hoffmann
Abstract:
The last two decades have seen a surge in kinetic and macroscopic models derived to investigate the multi-scale aspects of self-organised biological aggregations. Because the individual-level details incorporated into the kinetic models (e.g., individual speeds and turning rates) make them somewhat difficult to investigate, one is interested in transforming these models into simpler macroscopic mo…
▽ More
The last two decades have seen a surge in kinetic and macroscopic models derived to investigate the multi-scale aspects of self-organised biological aggregations. Because the individual-level details incorporated into the kinetic models (e.g., individual speeds and turning rates) make them somewhat difficult to investigate, one is interested in transforming these models into simpler macroscopic models, by using various scaling techniques that are imposed by the biological assumptions of the models. Here, we consider three scaling approaches (parabolic, hydrodynamic and grazing collision limits) that can be used to reduce a class of non-local 1D and 2D models for biological aggregations to simpler models existent in the literature. Next, we investigate how some of the spatio-temporal patterns exhibited by the original kinetic models are preserved via these scalings. To this end, we focus on the parabolic scaling for non-local 1D models and apply asymptotic preserving numerical methods, which allow us to analyse changes in the patterns as the scaling coefficient $ε$ is varied from $ε=1$ (for 1D transport models) to $ε=0$ (for 1D parabolic models). We show that some patterns (describing stationary aggregations) are preserved in the limit $ε\to 0$, while other patterns (describing moving aggregations) are lost in this limit. To understand the loss of these patterns, we construct bifurcation diagrams.
△ Less
Submitted 8 July, 2014;
originally announced July 2014.
-
Measuring the growth of matter fluctuations with third-order galaxy correlations
Authors:
Kai Hoffmann,
Julien Bel,
Enrique Gaztanaga,
Martin Crocce,
Pablo Fosalba,
Francisco J. Castander
Abstract:
Measurements of the linear growth factor $D$ at different redshifts $z$ are key to distinguish among cosmological models. One can estimate the derivative $dD(z)/d\ln(1+z)$ from redshift space measurements of the 3D anisotropic galaxy two-point correlation $ξ(z)$, but the degeneracy of its transverse (or projected) component with galaxy bias $b$, i.e. $ξ_{\perp}(z) \propto\ D^2(z) b^2(z)$, introduc…
▽ More
Measurements of the linear growth factor $D$ at different redshifts $z$ are key to distinguish among cosmological models. One can estimate the derivative $dD(z)/d\ln(1+z)$ from redshift space measurements of the 3D anisotropic galaxy two-point correlation $ξ(z)$, but the degeneracy of its transverse (or projected) component with galaxy bias $b$, i.e. $ξ_{\perp}(z) \propto\ D^2(z) b^2(z)$, introduces large errors in the growth measurement. Here we present a comparison between two methods which break this degeneracy by combining second- and third-order statistics. One uses the shape of the reduced three-point correlation and the other a combination of third-order one- and two-point cumulants. These methods use the fact that, for Gaussian initial conditions and scales larger than $20$ $h^{-1}$Mpc, the reduced third-order matter correlations are independent of redshift (and therefore of the growth factor) while the third-order galaxy correlations depend on $b$. We use matter and halo catalogs from the MICE-GC simulation to test how well we can recover $b(z)$ and therefore $D(z)$ with these methods in 3D real space. We also present a new approach, which enables us to measure $D$ directly from the redshift evolution of second- and third-order galaxy correlations without the need of modelling matter correlations. For haloes with masses lower than $10^{14}$ $h^{-1}$M$_\odot$, we find $10%$ deviations between the different estimates of $D$, which are comparable to current observational errors. At higher masses we find larger differences that can probably be attributed to the breakdown of the bias model and non-Poissonian shot noise.
△ Less
Submitted 17 December, 2014; v1 submitted 5 March, 2014;
originally announced March 2014.
-
Subhaloes gone Notts: Subhaloes as tracers of the dark matter halo shape
Authors:
Kai Hoffmann,
Susana Planelles,
Enrique Gaztanaga,
Alexander Knebe,
Frazer R. Pearce,
Hanni Lux,
Julian Onions,
Stuart I. Muldrew,
Pascal Elahi,
Peter Behroozi,
Yago Ascasibar,
Jiaxin Han,
Michal Maciejewski,
Manuel E. Merchan,
Mark Neyrinck,
Andrés N. Ruiz,
Mario A. Sgro
Abstract:
We study the shapes of subhalo distributions from four dark-matter-only simulations of Milky Way type haloes. Comparing the shapes derived from the subhalo distributions at high resolution to those of the underlying dark matter fields we find the former to be more triaxial if theanalysis is restricted to massive subhaloes. For three of the four analysed haloes the increased triaxiality of the dist…
▽ More
We study the shapes of subhalo distributions from four dark-matter-only simulations of Milky Way type haloes. Comparing the shapes derived from the subhalo distributions at high resolution to those of the underlying dark matter fields we find the former to be more triaxial if theanalysis is restricted to massive subhaloes. For three of the four analysed haloes the increased triaxiality of the distributions of massive subhaloes can be explained by a systematic effect caused by the low number of objects. Subhaloes of the fourth halo show indications for anisotropic accretion via their strong triaxial distribution and orbit alignment with respect to the dark matter field. These results are independent of the employed subhalo finder. Comparing the shape of the observed Milky Way satellite distribution to those of high-resolution subhalo samples from simulations, we find an agreement for samples of bright satellites, but significant deviations if faint satellites are included in the analysis. These deviations might result from observational incompleteness.
△ Less
Submitted 10 September, 2014; v1 submitted 9 January, 2014;
originally announced January 2014.
-
Continuity in time of solutions of a phase-field model
Authors:
Thomas G. Amler,
Nikolai D. Botkin,
Karl-Heinz Hoffmann,
Ibrahim Hoteit
Abstract:
A phase field model proposed by G. Caginalp for the description of phase changes in materials is under consideration. It is assumed that the medium is located in a container with heat conductive walls that are not subjected to phase changes. Therefore, the temperature variable is defined both in the medium and wall regions, whereas the phase variable is only considered in the medium part. The case…
▽ More
A phase field model proposed by G. Caginalp for the description of phase changes in materials is under consideration. It is assumed that the medium is located in a container with heat conductive walls that are not subjected to phase changes. Therefore, the temperature variable is defined both in the medium and wall regions, whereas the phase variable is only considered in the medium part. The case of Lipschitz domains in two and three dimensions is studied. We show that the temperature and phase variables are continuous in time functions with values in $L^2$ and $H^1$, respectively, provided that the initial values of them are from $L^2$ and $H^1$, respectively. Moreover, continuous dependence of solutions on the initial data and boundary conditions is proved.
△ Less
Submitted 1 August, 2013;
originally announced August 2013.
-
Measuring the Cosmic Web
Authors:
Volker Mueller,
Kai Hoffmann,
Sebastian E. Nuza
Abstract:
A quantitative study of the clustering properties of the cosmic web as a function of absolute magnitude and colour is presented using the SDSS Data Release 7 galaxy survey. Mark correlations are included in the analysis. We compare our results with mock galaxy samples obtained with four different semi-analytical models of galaxy formation imposed on the merger trees of the Millenium simulation. Th…
▽ More
A quantitative study of the clustering properties of the cosmic web as a function of absolute magnitude and colour is presented using the SDSS Data Release 7 galaxy survey. Mark correlations are included in the analysis. We compare our results with mock galaxy samples obtained with four different semi-analytical models of galaxy formation imposed on the merger trees of the Millenium simulation. The clustering of both red and blue galaxies is studied separately.
△ Less
Submitted 7 September, 2011;
originally announced September 2011.
-
Lorentz shift measurements in heavily irradiated silicon detectors in high magnetic fields
Authors:
Wim de Boer,
Karl-Heinz Hoffmann,
Andreas Sabellek,
Mike Schmanau,
Michael Schneider,
Valery Zhukov,
Theo Schneider
Abstract:
An external magnetic field exerts a Lorentz force on drifting electric charges inside a silicon strip sensor and thus shifts the cluster position of the collected charge. The shift can be related to the Lorentz angle which is typically a few degrees for holes and a few tens of degrees for electrons in a 4 T magnetic field. The Lorentz angle depends upon magnetic field, electric field inside the se…
▽ More
An external magnetic field exerts a Lorentz force on drifting electric charges inside a silicon strip sensor and thus shifts the cluster position of the collected charge. The shift can be related to the Lorentz angle which is typically a few degrees for holes and a few tens of degrees for electrons in a 4 T magnetic field. The Lorentz angle depends upon magnetic field, electric field inside the sensor and temperature. In this study the sensitivity to radiation for fluences up to 10^16 n/cm^2 has been studied. The Lorentz shift has been measured by inducing ionization with 670 nm red or 1070 nm infrared laser beams injected into the back side of the irradiated silicon sensor operated in magnetic fields up to 8 T. For holes the shift as a function of radiation is increasing, while for electrons it is decreasing and even changes sign. The fact that for irradiated sensors the Lorentz shift for electrons is smaller than for holes, in contrast to the observations in non-irradiated sensors, can be qualitatively explained by the structure of the electric field in irradiated sensors.
△ Less
Submitted 20 May, 2010;
originally announced May 2010.
-
A Hadronization Model for Few-GeV Neutrino Interactions
Authors:
T. Yang,
C. Andreopoulos,
H. Gallagher,
K. Hoffmann,
P. Kehayias
Abstract:
We present a detailed description of a new hadronic multiparticle production model for use in neutrino interaction simulations. Its validity spans a wide invariant mass range starting from pion production threshold. This model focuses on the low invariant mass region which is probed in few-GeV neutrino interactions and is of particular importance to neutrino oscillation experiments using acceler…
▽ More
We present a detailed description of a new hadronic multiparticle production model for use in neutrino interaction simulations. Its validity spans a wide invariant mass range starting from pion production threshold. This model focuses on the low invariant mass region which is probed in few-GeV neutrino interactions and is of particular importance to neutrino oscillation experiments using accelerator and atmospheric fluxes. It exhibits reasonable agreement with a wide variety of experimental data. We also describe measurements that can be made in upcoming experiments that can improve modeling in areas where uncertainties are currently large.
△ Less
Submitted 26 April, 2009;
originally announced April 2009.
-
The Quantum Refrigerator: The quest for absolute zero
Authors:
Yair Rezek,
Peter Salamon,
Karl Heinz Hoffmann,
Ronnie Kosloff
Abstract:
The scaling of the optimal cooling power of a reciprocating quantum refrigerator is sought as a function of the cold bath temperature as $T_c \to 0$. The working medium consists of noninteracting particles in a harmonic potential. Two closed-form solutions of the refrigeration cycle are analyzed, and compared to a numerical optimization scheme, focusing on cooling toward zero temperature. The op…
▽ More
The scaling of the optimal cooling power of a reciprocating quantum refrigerator is sought as a function of the cold bath temperature as $T_c \to 0$. The working medium consists of noninteracting particles in a harmonic potential. Two closed-form solutions of the refrigeration cycle are analyzed, and compared to a numerical optimization scheme, focusing on cooling toward zero temperature. The optimal cycle is characterized by linear relations between the heat extracted from the cold bath, the energy level spacing of the working medium and the temperature. The scaling of the optimal cooling rate is found to be proportional to $T_c^{3/2}$ giving a dynamical interpretation to the third law of thermodynamics.
△ Less
Submitted 4 August, 2008;
originally announced August 2008.
-
Desiccation of a clay film: Cracking versus peeling
Authors:
Supti Sadhukhan,
Janett Prehl,
Peter Blaudeck,
K. H. Hoffmann,
Tapati Dutta,
Sujata Tarafdar
Abstract:
Cracking and peeling of a layer of clay on desiccation has been simulated using a spring model. A vertical section through the layer with finite thickness is represented by a rectangular array of nodes connected by linear springs on a square lattice. The effect of reduction of the natural length of the springs, which mimics the drying is studied. Varying the strength of adhesion between sample a…
▽ More
Cracking and peeling of a layer of clay on desiccation has been simulated using a spring model. A vertical section through the layer with finite thickness is represented by a rectangular array of nodes connected by linear springs on a square lattice. The effect of reduction of the natural length of the springs, which mimics the drying is studied. Varying the strength of adhesion between sample and substrate and the rate of penetration of the drying front produces an interesting phase diagram, showing cross-over from peeling to cracking behavior. Changes in the number and width of cracks on varying the layer thickness is observed to reproduce experimental reports.
△ Less
Submitted 18 May, 2008;
originally announced May 2008.
-
Intermittent relaxation in hierarchical energy landscapes
Authors:
Andreas Fischer,
Paolo Sibani,
Karl Heinz Hoffmann
Abstract:
We numerically simulate a thermalization process in an energy landscape with hierarchically organized metastable states. The initial configuration is chosen to have a large energy excess, relative to the thermal equilibrium value at the running temperature. We show that the initial energy surplus is dissipated in a series of intermittent bursts, or quakes, whose rate decreases as the inverse of…
▽ More
We numerically simulate a thermalization process in an energy landscape with hierarchically organized metastable states. The initial configuration is chosen to have a large energy excess, relative to the thermal equilibrium value at the running temperature. We show that the initial energy surplus is dissipated in a series of intermittent bursts, or quakes, whose rate decreases as the inverse of the age of the system. In addition, one observes energy fluctuations with a zero centered Gaussian distribution. These pertain to the pseudo equilibrium dynamics within a single metastable state, and do not contribute to the energy dissipation. The derivative of the thermal energy with respect to the logarithm of time is asymptotically constant, and comprises a temperature independent part, and a part with an Arrhenius temperature dependence. The findings closely mirror recent numerical simulation results obtained for microscopic glassy models. For these models, record-sized energy fluctuations have been claimed to trigger intermittent events during low temperature thermalization. In the present model record-sized fluctuations are by construction needed to trigger changes from one metastable state to another. This property thus suffices to explain the statistical property of intermittent energy flow in complex metastable systems.
△ Less
Submitted 3 April, 2008; v1 submitted 22 September, 2007;
originally announced September 2007.
-
Optimization by thermal cycling
Authors:
A. Mobius,
A. Neklioudov,
A. Diaz-Sanchez,
K. H. Hoffmann,
A. Fachat,
M. Schreiber
Abstract:
An optimization algorithm is presented which consists of cyclically heating and quenching by Metropolis and local search procedures, respectively. It works particularly well when it is applied to an archive of samples instead of to a single one. We demonstrate for the traveling salesman problem that this algorithm is far more efficient than usual simulated annealing; our implementation can compe…
▽ More
An optimization algorithm is presented which consists of cyclically heating and quenching by Metropolis and local search procedures, respectively. It works particularly well when it is applied to an archive of samples instead of to a single one. We demonstrate for the traveling salesman problem that this algorithm is far more efficient than usual simulated annealing; our implementation can compete concerning speed with recent, very fast genetic local search algorithms.
△ Less
Submitted 9 February, 2005;
originally announced February 2005.
-
Optimization by thermal cycling
Authors:
A. Mobius,
K. H. Hoffmann,
C. Schon
Abstract:
Thermal cycling is an heuristic optimization algorithm which consists of cyclically heating and quenching by Metropolis and local search procedures, respectively, where the amplitude slowly decreases. In recent years, it has been successfully applied to two combinatorial optimization tasks, the traveling salesman problem and the search for low-energy states of the Coulomb glass. In these cases,…
▽ More
Thermal cycling is an heuristic optimization algorithm which consists of cyclically heating and quenching by Metropolis and local search procedures, respectively, where the amplitude slowly decreases. In recent years, it has been successfully applied to two combinatorial optimization tasks, the traveling salesman problem and the search for low-energy states of the Coulomb glass. In these cases, the algorithm is far more efficient than usual simulated annealing. In its original form the algorithm was designed only for the case of discrete variables. Its basic ideas are applicable also to a problem with continuous variables, the search for low-energy states of Lennard-Jones clusters.
△ Less
Submitted 9 February, 2005;
originally announced February 2005.
-
Nonlinear acousto-electric transport in a two-dimensional electron system
Authors:
A. O. Govorov,
A. V. Kalameitsev,
M. Rotter,
A. Wixforth,
J. P. Kotthaus,
K. H. Hoffmann,
N. Botkin
Abstract:
We study both theoretically and experimentally the nonlinear interaction between an intense surface acoustic wave and a two-dimensional electron plasma in semiconductor-piezocrystal hybrid structures. The experiments on hybrid systems exhibit strongly nonlinear acousto-electric effects. The plasma turns into moving electron stripes, the acousto-electric current reaches its maximum, and the sound…
▽ More
We study both theoretically and experimentally the nonlinear interaction between an intense surface acoustic wave and a two-dimensional electron plasma in semiconductor-piezocrystal hybrid structures. The experiments on hybrid systems exhibit strongly nonlinear acousto-electric effects. The plasma turns into moving electron stripes, the acousto-electric current reaches its maximum, and the sound absorption strongly decreases. To describe the nonlinear phenomena, we develop a coupled-amplitude method for a two-dimensional system in the strongly nonlinear regime of interaction. At low electron densities the absorption coefficient decreases with increasing sound intensity, whereas at high electron density the absorption coefficient is not a monotonous function of the sound intensity. High-harmonic generation coefficients as a function of the sound intensity have a nontrivial behavior. Theory and experiment are found to be in a good agreement.
△ Less
Submitted 6 May, 2000;
originally announced May 2000.
-
The state space of short-range Ising spin glasses: the density of states
Authors:
T. Klotz,
S. Schubert,
K. H. Hoffmann
Abstract:
The state space of finite square and cubic Ising spin glass models is analysed in terms of the global and the local density of states. Systems with uniform and gaussian probability distribution of interactions are compared. Different measures for the local state density are presented and discussed. In particular the question whether the local density of states grows exponentially or not is consi…
▽ More
The state space of finite square and cubic Ising spin glass models is analysed in terms of the global and the local density of states. Systems with uniform and gaussian probability distribution of interactions are compared. Different measures for the local state density are presented and discussed. In particular the question whether the local density of states grows exponentially or not is considered. The direct comparison of global and local densities leads to consequences for the structure of the state space.
△ Less
Submitted 15 October, 1997;
originally announced October 1997.