Search

Transiting planets with orbital periods longer than 40 d are extremely rare among the 5000+ planets discovered so far. The lack of discoveries of this population poses a challenge to research into planetary demographics, formation, and evolution. Here, we present the detection and characterization of HD88986b, a potentially transiting sub-Neptune, possessing the longest orbital period among known transiting small planets (< 4 R$_{\oplus}$) with a precise mass measurement ($\sigma_M/M$ > 25%). Additionally, we identified the presence of a massive companion in a wider orbit around HD88986. Our analysis reveals that HD88986b, based on two potential single transits on sector 21 and sector 48 which are both consistent with the predicted transit time from the RV model, is potentially transiting. The joint analysis of RV and photometric data show that HD88986b has a radius of 2.49$\pm$0.18 R$_{\oplus}$, a mass of 17.2$^{+4.0}_{-3.8}$ M$_{\oplus}$, and it orbits every 146.05$^{+0.43}_{-0.40}$ d around a subgiant HD88986 which is one of the closest and brightest exoplanet host stars (G2V type, R=1.543 $\pm$0.065 R$_{\odot}$, V=$6.47\pm 0.01$ mag, distance=33.37$\pm$0.04 pc). The nature of the outer, massive companion is still to be confirmed; a joint analysis of RVs, Hipparcos, and Gaia astrometric data shows that with a 3$\sigma$ confidence interval, its semi-major axis is between 16.7 and 38.8 au and its mass is between 68 and 284 M$_{Jup}$. HD88986b's wide orbit suggests the planet did not undergo significant mass loss due to extreme-ultraviolet radiation from its host star. Therefore, it probably maintained its primordial composition, allowing us to probe its formation scenario. Furthermore, the cold nature of HD88986b (460$\pm$8 K), thanks to its long orbital period, will open up exciting opportunities for future studies of cold atmosphere composition characterization., Comment: 37 pages, accepted to be published in A&A

This paper concerns structure learning or discovery of discrete generative models. It focuses on Bayesian model selection and the assimilation of training data or content, with a special emphasis on the order in which data are ingested. A key move - in the ensuing schemes - is to place priors on the selection of models, based upon expected free energy. In this setting, expected free energy reduces to a constrained mutual information, where the constraints inherit from priors over outcomes (i.e., preferred outcomes). The resulting scheme is first used to perform image classification on the MNIST dataset to illustrate the basic idea, and then tested on a more challenging problem of discovering models with dynamics, using a simple sprite-based visual disentanglement paradigm and the Tower of Hanoi (cf., blocks world) problem. In these examples, generative models are constructed autodidactically to recover (i.e., disentangle) the factorial structure of latent states - and their characteristic paths or dynamics.

Large-scale exoplanet surveys like the TESS mission are powerful tools for discovering large numbers of exoplanet candidates. Single-transit events are commonplace within the resulting candidate list due to the unavoidable limitation of observing baseline. These single-transit planets often remain unverified due to their unknown orbital period and consequent difficulty in scheduling follow up observations. In some cases, radial velocity (RV) follow up can constrain the period enough to enable a future targeted transit detection. We present the confirmation of one such planet: TOI-2010 b. Nearly three years of RV coverage determined the period to a level where a broad window search could be undertaken with the Near-Earth Object Surveillance Satellite (NEOSSat), detecting an additional transit. An additional detection in a much later TESS sector solidified our final parameter estimation. We find TOI-2010 b to be a Jovian planet ($M_P = 1.29 \ M_{\rm Jup}$, $R_P = 1.05 \ R_{\rm Jup}$) on a mildly eccentric orbit ($e = 0.21$) with a period of $P = 141.83403$ days. Assuming a simple model with no albedo and perfect heat redistribution, the equilibrium temperature ranges from about 360 K to 450 K from apoastron to periastron. Its wide orbit and bright host star ($V=9.85$) make TOI-2010 b a valuable test-bed for future low-insolation atmospheric analysis., Comment: 27 pages, 10 figures, 6 tables

Anderson localization, i.e., destructive quantum interference of multiple-scattering paths, halts transport entirely. Contrarily, time-dependent random forces expedite transport via Fermi acceleration, proposed as a mechanism for high-energy cosmic rays. Their competition creates interesting dynamics, but experimental observations are scarce. Here, we experimentally study the expansion of an ultracold Fermi gas inside time-dependent disorder and observe distinct regimes from sub- to superdiffusion. Unexpectedly, quantum interference counteracts acceleration in strong disorder. Our system enables the investigation of Fermi acceleration in the quantum-transport regime., Comment: Main: 7 pages, 3 figures; Supp. Mat.: 4 pages, 2 figures

Disorder can fundamentally modify the transport properties of a system. A striking example is Anderson localization, suppressing transport due to destructive interference of propagation paths. In inhomogeneous many-body systems, not all particles are localized for finite-strength disorder, and the system can become partially diffusive. Unravelling the intricate signatures of localization from such observed diffusion is a long-standing problem. Here, we experimentally study a degenerate, spin-polarized Fermi gas in a disorder potential formed by an optical speckle pattern. We record the diffusion in the disordered potential upon release from an external confining potential. We compare different methods to analyze the resulting density distributions, including a new method to capture particle dynamics by evaluating absorption-image statistics. Using standard observables, such as diffusion exponent and coefficient, localized fraction, or localization length, we find that some show signatures for a transition to localization above a critical disorder strength, while others show a smooth crossover to a modified diffusion regime. In laterally displaced disorder, we spatially resolve different transport regimes simultaneously which allows us to extract the subdiffusion exponent expected for weak localization. Our work emphasizes that the transition toward localization can be investigated by closely analyzing the system's diffusion, offering ways of revealing localization effects beyond the signature of exponentially decaying density distribution., Comment: 19 pages, 11 figures

Planetary systems around solar analogs inform us about how planets form and evolve in Solar System-like environments. We report the detection and characterization of two planetary systems around the solar analogs TOI-1736 and TOI-2141 using TESS photometry data and spectroscopic data obtained with the SOPHIE instrument on the 1.93 m telescope at the Observatoire de Haute-Provence (OHP). We performed a detailed spectroscopic analysis of these systems to obtain the precise radial velocities (RV) and physical properties of their host stars. TOI-1736 and TOI-2141 each host a transiting sub-Neptune with radii of $2.44\pm0.18$ R$_{\oplus}$ and $3.05\pm0.23$ R$_{\oplus}$, orbital periods of $7.073088(7)$ d and $18.26157(6)$ d, and masses of $12.8\pm1.8$ M$_{\oplus}$ and $24\pm4$ M$_{\oplus}$, respectively. TOI-1736 shows long-term RV variations that are consistent with a two-planet solution plus a linear trend of $-0.177$ ms$^{-1}$d$^{-1}$. We measured an RV semi-amplitude of $201.1\pm0.7$ ms$^{-1}$ for the outer companion, TOI-1736 c, implying a projected mass of $m_{c}\sin{i}=8.09\pm0.20$ M$_{\rm Jup}$. From the GAIA DR3 astrometric excess noise, we constrained the mass of TOI-1736 c at $8.7^{+1.5}_{-0.6}$ M$_{\rm Jup}$. This planet is in an orbit of $570.2\pm0.6$ d with an eccentricity of $0.362\pm0.003$ and a semi-major axis of $1.381\pm0.017$ au, where it receives a flux of $0.71\pm0.08$ times the bolometric flux incident on Earth, making it an interesting case of a supergiant planet that has settled into an eccentric orbit in the habitable zone of a solar analog. Our analysis of the mass-radius relation for the transiting sub-Neptunes shows that both TOI-1736 b and TOI-2141 b likely have an Earth-like dense rocky core and a water-rich envelope., Comment: Accepted for publication in A&A on October 6, 2023

Recent observations suggest the presence of clouds in exoplanet atmospheres but have also shown that certain chemical species in the upper atmosphere might not be in chemical equilibrium. The goal of this work is to calculate the two main cloud formation processes, nucleation and bulk growth, consistently from a non-equilibrium gas-phase. The aim is further to explore the interaction between a kinetic gas-phase and cloud micro-physics. The cloud formation is modeled using the moment method and kinetic nucleation which are coupled to a gas-phase kinetic rate network. Specifically, the formation of cloud condensation nuclei is derived from cluster rates that include the thermochemical data of (TiO$_2$)$_N$ from N = 1 to 15. The surface growth of 9 bulk Al/Fe/Mg/O/Si/S/Ti binding materials considers the respective gas-phase species through condensation and surface reactions as derived from kinetic disequilibrium. The effect of completeness of rate networks and the time evolution of the cloud particle formation is studied for an example exoplanet HD 209458 b. A consistent, fully time-dependent cloud formation model in chemical disequilibrium with respect to nucleation, bulk growth and the gas-phase is presented and first test cases are studied. This model shows that cloud formation in exoplanet atmospheres is a fast process. This confirms previous findings that the formation of cloud particles is a local process. Tests on selected locations within the atmosphere of the gas-giant HD 209458 b show that the cloud particle number density and volume reach constant values within 1s. The complex kinetic polymer nucleation of TiO$_2$ confirms results from classical nucleation models. The surface reactions of SiO[s] and SiO$_2$[s] can create a catalytic cycle that dissociates H$_2$ to 2 H, resulting in a reduction of the CH$_4$ number densities., Comment: 21 pages, 13 figures. Accepted for publication in Astronomy and Astrophysics (A&A)

We present an efficient cooling scheme for fermionic $^{40}{\text{K}}$ atoms, using laser light red and blue detuned with respect to the $\mathcal{D}2$ and $\mathcal{D}1$ principle flourescence lines, respectively. The cooling scheme is found to significantly increase the saturation level for loading of a $^{40}{\text{K}}$ magneto-optical trap (MOT), resulting in increased atom numbers or decreased cycle times. While the attainable $^{40}{\text{K}}$ atom number is approximately doubled if exclusively $^{40}{\text{K}}$ atoms are cooled, the scheme is particularly powerful for dual-species MOTs, for example, if $^{40}{\text{K}}$ and $^{87}{\text{Rb}}$ atoms are cooled simultaneously in the same MOT configuration. The typical atom losses due to light-assisted hetero-nuclear collisions between $^{40}{\text{K}}$ and $^{87}{\text{Rb}}$ seem to be reduced giving rise to a threefold improvement of the $^{40}{\text{K}}$ atom number as compared to that in a conventional dual-species MOT, operating merely with $\mathcal{D}2$ light. Our scheme can be a useful extension to most dual-species experiments, aiming to reach simultaneous degeneracy of both species., Comment: 8 pages, 5 figures

Sparse direct linear solvers are at the computational core of domain decomposition preconditioners and therefore have a strong impact on their performance. In this paper, we consider the Fast and Robust Overlapping Schwarz (FROSch) solver framework of the Trilinos software library, which contains a parallel implementations of the GDSW domain decomposition preconditioner. We compare three different sparse direct solvers used to solve the subdomain problems in FROSch. The preconditioner is applied to different model problems; linear elasticity and more complex fully-coupled deformation diffusion-boundary value problems from chemo-mechanics. We employ FROSch in fully algebraic mode, and therefore, we do not expect numerical scalability. Strong scalability is studied from 64 to 4096 cores, where good scaling results are obtained up to 1728 cores. The increasing size of the coarse problem increases the solution time for all sparse direct solvers.

To produce sounds, we adjust the tension of our vocal folds to shape their properties and control the pitch. This efficient mechanism offers inspiration for designing reconfigurable materials and adaptable soft robots. However, understanding how flexible structures respond to a significant static strain is not straightforward. This complexity also limits the precision of medical imaging when applied to tensioned organs like muscles, tendons, ligaments and blood vessels among others. In this article, we experimentally and theoretically explore the dynamics of a soft strip subject to a substantial static extension, up to 180%. Our observations reveal a few intriguing effects, such as the resilience of certain vibrational modes to a static deformation. These observations are supported by a model based on the incremental displacement theory. This has promising practical implications for characterizing soft materials but also for scenarios where external actions can be used to tune properties.

The dispersion curves of (elastic) waveguides frequently exhibit crossings and osculations (also known as veering, repulsion, or avoided crossing). Osculations are regions in the dispersion diagram where curves approach each other arbitrarily closely without ever crossing before veering apart. In semi-analytical (undamped) waveguide models, dispersion curves are obtained as solutions to discretized parameterized Hermitian eigenvalue problems. In the mathematical literature, it is known that such eigencurves can exhibit crossing points only if the corresponding matrix flow (parameter-dependent matrix) is uniformly decomposable. We discuss the implications for the solution of the waveguide problem. In particular, we make use of a simple algorithm recently suggested in the literature for decomposing matrix flows. We also employ a method for mode tracing based on approximating the eigenvalue problem for individual modes by an ordinary differential equation that can be solved by standard procedures., Comment: accepted version, published in Ultrasonics

The discovery of planets orbiting at less than 1 au from their host star and less massive than Saturn in various exoplanetary systems revolutionized our theories of planetary formation. The fundamental question is whether these close-in low-mass planets could have formed in the inner disk interior to 1 au, or whether they formed further out in the planet-forming disk and migrated inward. Exploring the role of additional giant planets in these systems may help us to pinpoint their global formation and evolution. We searched for additional substellar companions by using direct imaging in systems known to host close-in small planets. The use of direct imaging complemented by radial velocity and astrometric detection limits enabled us to explore the giant planet and brown dwarf demographics around these hosts to investigate the potential connection between both populations. We carried out a direct imaging survey with VLT/SPHERE to look for outer giant planets and brown dwarf companions in 27 systems hosting close-in low-mass planets discovered by radial velocity. Our sample is composed of very nearby (<20pc) planetary systems, orbiting G-, K-, and M-type mature (0.5-10Gyr) stellar hosts. We performed homogeneous direct imaging data reduction and analysis to search for and characterize point sources, and derived robust statistical detection limits. Of 337 point-source detections, we do not find any new bound companions. We recovered the emblematic very cool T-type brown dwarf GJ229B. Our typical sensitivities in direct imaging range from 5 to 30 MJup beyond 2 au. The non-detection of massive companions is consistent with predictions based on models of planet formation by core accretion. Our pilot study opens the way to a multi-technique approach for the exploration of very nearby exoplanetary systems with future ground-based and space observatories., Comment: 49 pages including 31 pages of appendices and references, 31 figures, A&A, accepted

The nascent but rapidly growing field of Quantum Information Science and Technology has led to an increased demand for skilled quantum workers and an opportunity to build a diverse workforce at the outset. In order to meet this demand and encourage women and underrepresented minorities in STEM to consider a career in QIST, we have developed a curriculum for introducing quantum computing to teachers and students at the high school level with no prerequisites. In 2022, this curriculum was delivered over the course of two one-week summer camps, one targeting teachers and another targeting students. Here, we present an overview of the objectives, curriculum, and activities, as well as results from the formal evaluation of both camps and the outlook for expanding QCaMP in future years., Comment: 9 pages, to be published in IEEE Xplore as part of the IEEE Quantum Week 2023 proceedings

The magnetic field of the Sun is generated by internal dynamo process with a cyclic period of 11 years or a 22 year magnetic cycle. The signatures of the Sun's magnetic cycle are observed in the different layers of its atmosphere and in its internal layers. In this review, we use the same diagnostics to understand the magnetic cycles of other stars with the same internal structure as the Sun. We review what is currently known about mapping the surface magnetic fields, chromospheric and coronal indicators, cycles in photometry and asteroseismology. We conclude our review with an outlook for the future., Comment: accepted by Space Science Reviews

Molecular isomerization kinetics in liquid solvents are determined by a complex interplay between the friction acting on a rotating dihedral due to interactions with the solvent, internal dissipation effects (also known as internal friction), the viscosity of the solvent, and the free energy profile over which a dihedral rotates. Currently, it is not understood how these quantities are related at the molecular scale. Here, we combine molecular dynamics simulations of isomerizing n-alkane chains and dipeptide molecules in mixed water-glycerol solvents with memory-kernel extraction techniques to directly evaluate the frequency-dependent friction acting on a rotating dihedral. We extract the friction and isomerization times over a range of glycerol concentrations and accurately evaluate the relationships between solvent viscosity, isomerization kinetics, and dihedral friction. We show that the total friction acting on a rotating dihedral does not scale linearly with solvent viscosity, thus violating the Stokes-Einstein relation. Additionally, we demonstrate that the kinetics of isomerization are significantly faster compared to the Kramers prediction in the overdamped limit. We suggest that isomerization kinetics are determined by the multi-time-scale friction coupling between a rotating dihedral and its solvent environment, which results in non-Markovian kinetic speed-up effects., Comment: 9 pages, 4 figures

Although the latent spaces learned by distinct neural networks are not generally directly comparable, recent work in machine learning has shown that it is possible to use the similarities and differences among latent space vectors to derive "relative representations" with comparable representational power to their "absolute" counterparts, and which are nearly identical across models trained on similar data distributions. Apart from their intrinsic interest in revealing the underlying structure of learned latent spaces, relative representations are useful to compare representations across networks as a generic proxy for convergence, and for zero-shot model stitching. In this work we examine an extension of relative representations to discrete state-space models, using Clone-Structured Cognitive Graphs (CSCGs) for 2D spatial localization and navigation as a test case. Our work shows that the probability vectors computed during message passing can be used to define relative representations on CSCGs, enabling effective communication across agents trained using different random initializations and training sequences, and on only partially similar spaces. We introduce a technique for zero-shot model stitching that can be applied post hoc, without the need for using relative representations during training. This exploratory work is intended as a proof-of-concept for the application of relative representations to the study of cognitive maps in neuroscience and AI., Comment: 19 pages, 1 table, 6 figures. Accepted paper at the 4th International Workshop on Active Inference (Ghent, Belgium 2023)

Quantum information can not be perfectly cloned, but approximate copies of quantum information can be generated. Quantum telecloning combines approximate quantum cloning, more typically referred as quantum cloning, and quantum teleportation. Quantum telecloning allows approximate copies of quantum information to be constructed by separate parties, using the classical results of a Bell measurement made on a prepared quantum telecloning state. Quantum telecloning can be implemented as a circuit on quantum computers using a classical co-processor to compute classical feed forward instructions using if statements based on the results of a mid-circuit Bell measurement in real time. We present universal, symmetric, optimal $1 \rightarrow M$ telecloning circuits, and experimentally demonstrate these quantum telecloning circuits for $M=2$ up to $M=10$, natively executed with real time classical control systems on IBM Quantum superconducting processors, known as dynamic circuits. We perform the cloning procedure on many different message states across the Bloch sphere, on $7$ IBM Quantum processors, optionally using the error suppression technique X-X sequence digital dynamical decoupling. Two circuit optimizations are utilized, one which removes ancilla qubits for $M=2, 3$, and one which reduces the total number of gates in the circuit but still uses ancilla qubits. Parallel single qubit tomography with MLE density matrix reconstruction is used in order to compute the mixed state density matrices of the clone qubits, and clone quality is measured using quantum fidelity. These results present one of the largest and most comprehensive NISQ computer experimental analyses on (single qubit) quantum telecloning to date. The clone fidelity sharply decreases to $0.5$ for $M > 5$, but for $M=2$ we are able to achieve a mean clone fidelity of up to $0.79$ using dynamical decoupling.

Atmospheric characterisation of exoplanets from the ground is an actively growing field of research. In this context we have created the ATMOSPHERIX consortium: a research project aimed at characterizing exoplanets atmospheres using ground-based high resolution spectroscopy. This paper presents the publicly-available data analysis pipeline and demonstrates the robustness of the recovered planetary parameters from synthetic data. Simulating planetary transits using synthetic transmission spectra of a hot Jupiter that were injected into real SPIRou observations of the non-transiting system Gl 15 A, we show that our pipeline is successful at recovering the planetary signal and input atmospheric parameters. We also introduce a deep learning algorithm to optimise data reduction which proves to be a reliable, alternative tool to the commonly used principal component analysis. We estimate the level of uncertainties and possible biases when retrieving parameters such as temperature and composition and hence the level of confidence in the case of retrieval from real data. Finally, we apply our pipeline onto two real transits of HD~189733 b observed with SPIRou and obtain similar results than in the literature. In summary, we have developed a publicly available and robust pipeline for the forthcoming studies of the targets to be observed in the framework of the ATMOSPHERIX consortium, which can easily be adapted to other high resolution instruments than SPIRou (e.g. VLT-CRIRES, MAROON-X, ELT-ANDES), Comment: Accepted in MNRAS

Large language models (LLMs) have recently gained popularity. However, the impact of their general availability through ChatGPT on sensitive areas of everyday life, such as education, remains unclear. Nevertheless, the societal impact on established educational methods is already being experienced by both students and educators. Our work focuses on higher physics education and examines problem solving strategies. In a study, students with a background in physics were assigned to solve physics exercises, with one group having access to an internet search engine (N=12) and the other group being allowed to use ChatGPT (N=27). We evaluated their performance, strategies, and interaction with the provided tools. Our results showed that nearly half of the solutions provided with the support of ChatGPT were mistakenly assumed to be correct by the students, indicating that they overly trusted ChatGPT even in their field of expertise. Likewise, in 42% of cases, students used copy & paste to query ChatGPT -- an approach only used in 4% of search engine queries -- highlighting the stark differences in interaction behavior between the groups and indicating limited reflection when using ChatGPT. In our work, we demonstrated a need to (1) guide students on how to interact with LLMs and (2) create awareness of potential shortcomings for users., Comment: Pre-print currently under review

Context. Although more than one thousand sub-stellar companions have already been detected with the radial velocity (RV) method, many new companions remain to be detected in the public RV archives. Aims. We wish to use the archival data obtained with the ESO/HARPS spectrograph to search for sub-stellar companions. Methods. We use the astronomic acceleration measurements of stars obtained with the Hipparcos and Gaia satellites to identify anomalies that could be explained by the presence of a companion. Once hints for a companion are found, we combine the RV data with absolute astrometry and, when available, relative astrometry data, using a Markov Chain Monte Carlo (MCMC) algorithm to determine the orbital parameters and mass of the companion. Results. We find and characterize three new brown dwarfs (GJ660.1 C, HD73256 B, and HD165131 B) and six new planets (HD75302 b, HD108202 b, HD135625 b, HD185283 b, HIP10337 b, and HIP54597 b) with separations between 1 and 6 au and masses between 0.6 and 100 MJup. We also constrain the orbital inclination of ten known sub-stellar companions and determine their true mass. Finally, we identify twelve new stellar companions. This shows that the analysis of proper motion anomalies allows for optimizing the RV search for sub-stellar companions and their characterization., Comment: 72 pages, 61 figures, 13 tables, accepted for publication by A&A

We report the discovery and characterisation of a giant transiting planet orbiting a nearby M3.5V dwarf (d = 80.4 pc, $G$ = 15.1 mag, $K$=11.2 mag, R$_\star$ = 0.358 $\pm$ 0.015 R$_\odot$, M$_\star$ = 0.340 $\pm$ 0.009 M$_\odot$). Using the photometric time series from TESS sectors 10, 36, 46, and 63 and near-infrared spectrophotometry from ExTrA, we measured a planetary radius of 0.77 $\pm$ 0.03 R$_J$ and an orbital period of 1.52 days. With high-resolution spectroscopy taken by the CFHT/SPIRou and ESO/ESPRESSO spectrographs, we refined the host star parameters ([Fe/H] = 0.27 $\pm$ 0.12) and measured the mass of the planet (0.273 $\pm$ 0.006 M$_J$). Based on these measurements, TOI-4860 b joins the small set of massive planets ($>$80 M$_E$) found around mid to late M dwarfs ($<$0.4 R$_\odot$), providing both an interesting challenge to planet formation theory and a favourable target for further atmospheric studies with transmission spectroscopy. We identified an additional signal in the radial velocity data that we attribute to an eccentric planet candidate ($e=0.66\pm0.09$) with an orbital period of $427\pm7$~days and a minimum mass of $1.66\pm 0.26$ M$_J$, but additional data would be needed to confirm this., Comment: 16 pages, 14 figures, accepted for publication in A&A

We study the injectivity of the Kudla-Millson lift of genus 2 Siegel cusp forms, vector-valued with respect to the Weil representation associated to an even lattice L. We prove that if L splits off two hyperbolic planes and is of sufficiently large rank, then the lift is injective. As an application, we deduce that the image of the lift in the degree 4 cohomology of the associated orthogonal Shimura variety has the same dimension as the lifted space of cusp forms. Our results also cover the case of moduli spaces of quasi-polarized K3 surfaces. To prove the injectivity, we introduce vector-valued indefinite Siegel theta functions of genus 2 and of Jacobi type attached to L. We describe their behavior with respect to the split of a hyperbolic pane in L. This generalizes results of Borcherds to genus higher than 1., Comment: 65 pages

Elastic waves in anisotropic media can exhibit a power flux that is not collinear with the wave vector. This has notable consequences for waves guided in a plate. Through laser-ultrasonic experiments, we evidence remarkable phenomena due to slow waves in a single crystal silicon wafer. Waves exhibiting power flux orthogonal to their wave vector are identified. A pulsed line source that excites these waves reveals a wave packet radiated parallel to the line. Furthermore, there exist precisely eight plane waves with zero power flux. These so-called zero-group-velocity modes are oriented along the crystal's principal axes. Time acts as a filter in the wave vector domain that selects these modes. Thus, a point source leads to beating resonance patterns with moving nodal curves on the surface of the infinite plate. We observe this pattern as it emerges naturally after a pulsed excitation., Comment: 11 pages, 8 figures

This paper describes a modification of Ehlers' model for the inelastic behavior of granular media. The modified model can be applied for describing the inelastic behavior of porous media. The key feature is a subtle change of the yield potential, which allows the correct orientation of the triangular-shaped yield surface cross sections depending on the hydrostatic stress state. The model is incorporated into a general framework for isotropic plasticity. An elastic predictor/corrector algorithm is employed to solve the constitutive equations. The necessary derivatives for a Newton update are also given in detail. The model is calibrated using stress, and strain data obtained from finite element simulations of a generic highly porous open-cell Wheire-Phelan foam.