Slope stability analysis of the Pacaya Volcano, Guatemala, using limit equilibrium and finite element method

The Pacaya volcanic complex is part of the Central American volcanic arc, which is associated with the subduction of the Cocos tectonic plate under the Caribbean plate. Located 30 km south of Guatemala City, Pacaya is situated on the southern rim of the Amatitlan Caldera. It is the largest post-caldera volcano, and has been one of Central America’s most active volcanoes over the last 500 years. Between 400 and 2000 years B.P, the Pacaya volcano had experienced a huge collapse, which resulted in the formation of horseshoe-shaped scarp that is still visible. In the recent years, several smaller collapses have been associated with the activity of the volcano (in 1961 and 2010) affecting its northwestern flanks, which are likely to be induced by the local and regional stress changes. The similar orientation of dry and volcanic fissures and the distribution of new vents would likely explain the reactivation of the pre-existing stress configuration responsible for the old-collapse. This paper presents the first stability analysis of the Pacaya volcanic flank. The inputs for the geological and geotechnical models were defined based on the stratigraphical, lithological, structural data, and material properties obtained from field survey and lab tests. According to the mechanical characteristics, three lithotechnical units were defined: Lava, Lava-Breccia and Breccia-Lava. The Hoek and Brown’s failure criterion was applied for each lithotechnical unit and the rock mass friction angle, apparent cohesion, and strength and deformation characteristics were computed in a specified stress range. Further, the stability of the volcano was evaluated by two-dimensional analysis performed by Limit Equilibrium (LEM, ROCSCIENCE) and Finite Element Method (FEM, PHASE 2 7.0). The stability analysis mainly focused on the modern Pacaya volcano built inside the collapse amphitheatre of “Old Pacaya”. The volcanic instability was assessed based on the variability of safety factor using deterministic, sensitivity, and probabilistic analysis considering the gravitational instability and the effects of external forces such as magma pressure and seismicity as potential triggering mechanisms of lateral collapse. The preliminary results from the analysis provide two insights: first, the least stable sector is on the south-western flank of the volcano; second, the lowest safety factor value suggests that the edifice is stable under gravity alone, and the external triggering mechanism can represent a likely destabilizing factor.

Pupil phase apodization for achromatic imaging of extra-solar planets

Direct imaging of extra-solar planets in the visible and infrared region has generated great interest among scientists and the general public as well. However, this is a challenging problem. Diffculties of detecting a planet (faint source) are caused, mostly, by two factors: sidelobes caused by starlight diffraction from the edge of the pupil and the randomly scattered starlight caused by the phase errors from the imperfections in the optical system. While the latter diffculty can be corrected by high density active deformable mirrors with advanced phase sensing and control technology, the optimized strategy for suppressing the diffraction sidelobes is still an open question. In this thesis, I present a new approach to the sidelobe reduction problem: pupil phase apodization. It is based on a discovery that an anti-symmetric spatial phase modulation pattern imposed over a pupil or a relay plane causes diffracted starlight suppression sufficient for imaging of extra-solar planets. Numerical simulations with specific square pupil (side D) phase functions, such as ... demonstrate annulling in at least one quadrant of the diffraction plane to the contrast level of better than 10^12 with an inner working angle down to 3.5L/D (with a = 3 and e = 10^3). Furthermore, our computer experiments show that phase apodization remains effective throughout a broad spectrum (60% of the central wavelength) covering the entire visible light range. In addition to the specific phase functions that can yield deep sidelobe reduction on one quadrant, we also found that a modified Gerchberg-Saxton algorithm can help to find small sized (101 x 101 element) discrete phase functions if regional sidelobe reduction is desired. Our simulation shows that a 101x101 segmented but gapless active mirror can also generate a dark region with Inner Working Distance about 2.8L/D in one quadrant. Phase-only modulation has the additional appeal of potential implementation via active segmented or deformable mirrors, thereby combining compensation of random phase aberrations and diffraction halo removal in a single optical element.

Reverse Engineering the Structural and Acoustic Behavior of a Stradivari Violin

There is a tremendous amount of mystery that surrounds the instruments of Antonio Stradivari. There have been many studies done in the past, but no one completely understands exactly how he made his instruments, or why they are still considered the best in the world. This project is designed to develop an engineering model of one of Stradivari's violins that will accurately simulate the structural and acoustic behavior of the instrument. It also hopes to shine some light on what makes the instruments of Stradivari unique when compared to other violins. It will focus on geometry and material properties, utilizing several modern engineering tools, including CT scanning, experimental modal analysis, finite element analysis, correlation techniques, and acoustic synthesis.

Two-Phase Flow In Microchannels: Morphology And Interface Phenomena

The existence and morphology, as well as the dynamics of micro-scale gas-liquid interfaces is investigated numerically and experimentally. These studies can be used to assess liquid management issues in microsystems such as PEMFC gas flow channels, and are meant to open new research perspectives in two-phase flow, particularly in film deposition on non-wetting surfaces. For example the critical plug volume data can be used to deliver desired length plugs, or to determine the plug formation frequency. The dynamics of gas-liquid interfaces, of interest for applications involving small passages (e.g. heat exchangers, phase separators and filtration systems), was investigated using high-speed microscopy - a method that also proved useful for the study of film deposition processes. The existence limit for a liquid plug forming in a mixed wetting channel is determined by numerical simulations using Surface Evolver. The plug model simulate actual conditions in the gas flow channels of PEM fuel cells, the wetting of the gas diffusion layer (GDL) side of the channel being different from the wetting of the bipolar plate walls. The minimum plug volume, denoted as critical volume is computed for a series of GDL and bipolar plate wetting properties. Critical volume data is meant to assist in the water management of PEMFC, when corroborated with experimental data. The effect of cross section geometry is assessed by computing the critical volume in square and trapezoidal channels. Droplet simulations show that water can be passively removed from the GDL surface towards the bipolar plate if we take advantage on differing wetting properties between the two surfaces, to possibly avoid the gas transport blockage through the GDL. High speed microscopy was employed in two-phase and film deposition experiments with water in round and square capillary tubes. Periodic interface destabilization was observed and the existence of compression waves in the gas phase is discussed by taking into consideration a naturally occurring convergent-divergent nozzle formed by the flowing liquid phase. The effect of channel geometry and wetting properties was investigated through two-phase water-air flow in square and round microchannels, having three static contact angles of 20, 80 and 105 degrees. Four different flow regimes are observed for a fixed flow rate, this being thought to be caused by the wetting behavior of liquid flowing in the corners as well as the liquid film stability. Film deposition experiments in wetting and non-wetting round microchannels show that a thicker film is deposited for wetting conditions departing from the ideal 0 degrees contact angle. A film thickness dependence with the contact angle theta as well as the Capillary number, in the form h_R ~ Ca^(2/3)/ cos(theta) is inferred from scaling arguments, for contact angles smaller than 36 degrees. Non-wetting film deposition experiments reveal that a film significantly thicker than the wetting Bretherton film is deposited. A hydraulic jump occurs if critical conditions are met, as given by a proposed nondimensional parameter similar to the Froude number. Film thickness correlations are also found by matching the measured and the proposed velocity derived in the shock theory. The surface wetting as well as the presence of the shock cause morphological changes in the Taylor bubble flow.

DEVELOPMENT OF CONFOCAL IMAGING TECHNIQUES FOR PROBING INTERFACIAL DYNAMICS IN MICROSCALE, GAS-LIQUID, TWO-PHASE FLOW

Micro-scale, two-phase flow is found in a variety of devices such as Lab-on-a-chip, bio-chips, micro-heat exchangers, and fuel cells. Knowledge of the fluid behavior near the dynamic gas-liquid interface is required for developing accurate predictive models. Light is distorted near a curved gas-liquid interface preventing accurate measurement of interfacial shape and internal liquid velocities. This research focused on the development of experimental methods designed to isolate and probe dynamic liquid films and measure velocity fields near a moving gas-liquid interface. A high-speed, reflectance, swept-field confocal (RSFC) imaging system was developed for imaging near curved surfaces. Experimental studies of dynamic gas-liquid interface of micro-scale, two-phase flow were conducted in three phases. Dynamic liquid film thicknesses of segmented, two-phase flow were measured using the RSFC and compared to a classic film thickness deposition model. Flow fields near a steadily moving meniscus were measured using RSFC and particle tracking velocimetry. The RSFC provided high speed imaging near the menisci without distortion caused the gas-liquid interface. Finally, interfacial morphology for internal two-phase flow and droplet evaporation were measured using interferograms produced by the RSFC imaging technique. Each technique can be used independently or simultaneously when.

Gas-phase Dissociation of homo-DNA Oligonucleotides

Synthetic modified oligonucleotides are of interest for diagnostic and therapeutic applications, as their biological stability, pairing selectivity, and binding strength can be considerably increased by the incorporation of unnatural structural elements. Homo-DNA is an oligonucleotide homologue based on dideoxy-hexopyranosyl sugar moieties, which follows the Watson-Crick A-T and G-C base pairing system, but does not hybridize with complementary natural DNA and RNA. Homo-DNA has found application as a bioorthogonal element in templated chemistry applications. The gas-phase dissociation of homo-DNA has been investigated by ESI-MS/MS and MALDI-MS/MS, and mechanistic aspects of its gas-phase dissociation are discussed. Experiments revealed a charge state dependent preference for the loss of nucleobases, which are released either as neutrals or as anions. In contrast to DNA, nucleobase loss from homo-DNA was found to be decoupled from backbone cleavage, thus resulting in stable products. This renders an additional stage of ion activation necessary in order to generate sequence-defining fragment ions. Upon MS(3) of the primary base-loss ion, homo-DNA was found to exhibit unspecific backbone dissociation resulting in a balanced distribution of all fragment ion series.

Finite element analysis for prediction of bone strength

Article preview View full access options BoneKEy Reports | Review Print Email Share/bookmark Finite element analysis for prediction of bone strength Philippe K Zysset, Enrico Dall'Ara, Peter Varga & Dieter H Pahr Affiliations Corresponding author BoneKEy Reports (2013) 2, Article number: 386 (2013) doi:10.1038/bonekey.2013.120 Received 03 January 2013 Accepted 25 June 2013 Published online 07 August 2013 Article tools Citation Reprints Rights & permissions Abstract Abstract• References• Author information Finite element (FE) analysis has been applied for the past 40 years to simulate the mechanical behavior of bone. Although several validation studies have been performed on specific anatomical sites and load cases, this study aims to review the predictability of human bone strength at the three major osteoporotic fracture sites quantified in recently completed in vitro studies at our former institute. Specifically, the performance of FE analysis based on clinical computer tomography (QCT) is compared with the ones of the current densitometric standards, bone mineral content, bone mineral density (BMD) and areal BMD (aBMD). Clinical fractures were produced in monotonic axial compression of the distal radii, vertebral sections and in side loading of the proximal femora. QCT-based FE models of the three bones were developed to simulate as closely as possible the boundary conditions of each experiment. For all sites, the FE methodology exhibited the lowest errors and the highest correlations in predicting the experimental bone strength. Likely due to the improved CT image resolution, the quality of the FE prediction in the peripheral skeleton using high-resolution peripheral CT was superior to that in the axial skeleton with whole-body QCT. Because of its projective and scalar nature, the performance of aBMD in predicting bone strength depended on loading mode and was significantly inferior to FE in axial compression of radial or vertebral sections but not significantly inferior to FE in side loading of the femur. Considering the cumulated evidence from the published validation studies, it is concluded that FE models provide the most reliable surrogates of bone strength at any of the three fracture sites.

An experimentally validated finite element method for augmented vertebral bodies

Background Finite element models of augmented vertebral bodies require a realistic modelling of the cement infiltrated region. Most methods published so far used idealized cement shapes or oversimplified material models for the augmented region. In this study, an improved, anatomy-specific, homogenized finite element method was developed and validated to predict the apparent as well as the local mechanical behavior of augmented vertebral bodies. Methods Forty-nine human vertebral body sections were prepared by removing the cortical endplates and scanned with high-resolution peripheral quantitative CT before and after injection of a standard and a low-modulus bone cement. Forty-one specimens were tested in compression to measure stiffness, strength and contact pressure distributions between specimens and loading-plates. From the remaining eight, fourteen cylindrical specimens were extracted from the augmented region and tested in compression to obtain material properties. Anatomy-specific finite element models were generated from the CT data. The models featured element-specific, density-fabric-based material properties, damage accumulation, real cement distributions and experimentally determined material properties for the augmented region. Apparent stiffness and strength as well as contact pressure distributions at the loading plates were compared between simulations and experiments. Findings The finite element models were able to predict apparent stiffness (R2 > 0.86) and apparent strength (R2 > 0.92) very well. Also, the numerically obtained pressure distributions were in reasonable quantitative (R2 > 0.48) and qualitative agreement with the experiments. Interpretation The proposed finite element models have proven to be an accurate tool for studying the apparent as well as the local mechanical behavior of augmented vertebral bodies.

In situ dehydration behavior of zeolite-like pentagonite: A single-crystal X-ray study

The structural modifications upon heating of pentagonite, Ca(VO)(Si4O10)·4H2O (space group Ccm21, a=10.3708(2), b=14.0643(2), c=8.97810(10) Å, V=1309.53(3) Å3) were investigated by in situ temperature dependent single-crystal X-ray structure refinements. Diffraction data of a sample from Poona district (India) have been measured in steps of 25 up to 250 °C and in steps of 50 °C between 250 and 400 °C. Pentagonite has a porous framework structure made up by layers of silicate tetrahedra connected by V4+O5 square pyramids. Ca and H2O molecules are extraframework occupants. Room temperature diffraction data allowed refinement of H positions. The hydrogen-bond system links the extraframework occupants to the silicate layers and also interconnects the H2O molecules located inside the channels. Ca is seven-fold coordinated forming four bonds to O of the tetrahedral framework and three bonds to extraframework H2O. The H2O molecule at O9 showing a high displacement parameter is not bonded to Ca. The dehydration in pentagonite proceeds in three steps. At 100 °C the H2O molecule at O8 was released while O9 moved towards Ca. As a consequence the displacement parameter of H2O at O9 halved compared to that at room temperature. The unit-cell volume decreased to 1287.33(3) Å3 leading to a formula with 3H2O per formula unit (pfu). Ca remained seven-fold coordinated. At 175 °C Ca(VO)(Si4O10)·3H2O transformed into a new phase with 1H2O molecule pfu characterized by doubling of the c axis and the monoclinic space group Pn. Severe bending of specific TOT angles led to contraction of the porous three-dimensional framework. In addition, H2O at O9 was expelled while H2O at O7 approached a position in the center of the channel. The normalized volume decreased to 1069.44(9) Å3. The Ca coordination reduced from seven- to six-fold. At 225 °C a new anhydrous phase with space group Pna21 but without doubling of c had formed. Release of H2O at O7 caused additional contraction of TOT angles and volume reduction (V=1036.31(9) Å3). Ca adopted five-fold coordination. During heating excursion up to 400 °C this anhydrous phase remained preserved. Between room temperature and 225 °C the unit-cell volume decreased by 21% due to dehydration. The dehydration steps compare well with the thermo-gravimetric data reported in the literature.

In situ dehydration behavior of veszelyite (Cu,Zn)2Zn(PO4)(OH)3·2H2O: A single-crystal X-ray study

The rare mixed copper-zinc phosphate mineral veszelyite (Cu,Zn)2Zn(PO4)(OH)3·2H2O space group P21/c, a = 7.5096(2), b = 10.2281(2), c = 9.8258(2) Å, β = 103.3040(10)°, V = 734.45(3) Å3 was investigated by in situ temperature-dependent single-crystal X-ray structure refinements. The atomic arrangement of veszelyite consists of an alternation of octahedral and tetrahedral sheets. The Jahn-Teller distorted CuO6 octahedra form sheets with eight-membered rings. The tetrahedral sheet composed of PO4 and ZnO3(OH) tetrahedra shows strong topological similarities to that of cavansite, gismondine, and kipushite.Diffraction data of a sample from Zdravo Vrelo, near Kreševo (Bosnia and Herzegovina) have been measured in steps of 25 up to 225 °C. Hydrogen positions and the hydrogen-bond system were determined experimentally from the structure refinements of data collected up to 125 °C. At 200 °C, the hydrogen-bonding scheme was inferred from bond-valence calculations and donor-acceptor distances. The hydrogen-bond system connects the tetrahedral sheet to the octahedral sheet and also braces the Cu sheet.At 150 °C, the H2O molecule at H2O2 was released and the Cu coordination (Cu1 and Cu2) decreased from originally six- to fivefold. Cu1 has a square planar coordination by four OH groups and an elongate distance to O3, whereas Cu2 has the Jahn-Teller characteristic elongate bond to H2O1. The unit-cell volume decreased 7% from originally 734.45(3) to 686.4(4) Å3 leading to a formula with 1 H2O pfu. The new phase observed above 150 °C is characterized by an increase of the c axis and a shortening of the b axis. The bending of T-O-T angles causes an increasing elliptical shape of the eight-membered rings in the tetrahedral and octahedral sheets. Moreover a rearrangement of the hydrogen-bond system was observed.At 225 °C, the structure degrades to an X-ray amorphous residual due to release of the last H2O molecule at H2O1. The stronger Jahn-Teller distortion of Cu1 relative to Cu2 suggests that Cu1 is fully occupied by Cu, whereas Cu2 bears significant Zn. H2O1 is the fifth ligand of Cu2. Zn at Cu2 is not favorable to adopt planar fourfold coordination. Thus, if the last water molecule is expelled the structure is destabilized.This study contributes to understanding the dehydration mechanism and thermal stability of supergene minerals characterized by Jahn-Teller distorted octahedra with mixed Cu, Zn occupancy.

Distortion-sensitive insight into the pretransitional behavior of 4- n -octyloxy-4′-cyanobiphenyl (8OCB)

Results of studies of the static and dynamic dielectric properties in rod-like 4-n-octyloxy-4'-cyanobiphenyl (8OCB) with isotropic (I)–nematic (N)–smectic A (SmA)–crystal (Cr) mesomorphism, combined with measurements of the low-frequency nonlinear dielectric effect and heat capacity are presented. The analysis is supported by the derivative-based and distortion-sensitive transformation of experimental data. Evidence for the I–N and N–SmA pretransitional anomalies, indicating the influence of tricritical behavior, is shown. It has also been found that neither the N phase nor the SmA phase are uniform and hallmarks of fluid–fluid crossovers can be detected. The dynamics, tested via the evolution of the primary relaxation time, is clearly non-Arrhenius and described via τ(T) = τc(T−TC)−phgr. In the immediate vicinity of the I–N transition a novel anomaly has been found: Δτ ∝ 1/(T − T*), where T* is the temperature of the virtual continuous transition and Δτ is the excess over the 'background behavior'. Experimental results are confronted with the comprehensive Landau–de Gennes theory based modeling.

Is there a common pattern of future gas-phase air pollution in Europe under diverse climate change scenarios?

Climate change alone influences future levels of tropospheric ozone and their precursors through modifications of gas-phase chemistry, transport, removal, and natural emissions. The goal of this study is to determine at what extent the modes of variability of gas-phase pollutants respond to different climate change scenarios over Europe. The methodology includes the use of the regional modeling system MM5 (regional climate model version)-CHIMERE for a target domain covering Europe. Two full-transient simulations covering from 1991–2050 under the SRES A2 and B2 scenarios driven by ECHO-G global circulation model have been compared. The results indicate that the spatial patterns of variability for tropospheric ozone are similar for both scenarios, but the magnitude of the change signal significantly differs for A2 and B2. The 1991–2050 simulations share common characteristics for their chemical behavior. As observed from the NO2 and α-pinene modes of variability, our simulations suggest that the enhanced ozone chemical activity is driven by a number of parameters, such as the warming-induced increase in biogenic emissions and, to a lesser extent, by the variation in nitrogen dioxide levels. For gas-phase pollutants, the general increasing trend for ozone found under A2 and B2 forcing is due to a multiplicity of climate factors, such as increased temperature, decreased wet removal associated with an overall decrease of precipitation in southern Europe, increased photolysis of primary and secondary pollutants as a consequence of lower cloudiness and increased biogenic emissions fueled by higher temperatures.

Finite Element Based Nonlinear Normalization of Human Lumbar Intervertebral Disc Stiffness to Account for Its Morphology

Disc degeneration, usually associated with low back pain and changes of intervertebral stiffness, represents a major health issue. As the intervertebral disc (IVD) morphology influences its stiffness, the link between mechanical properties and degenerative grade is partially lost without an efficient normalization of the stiffness with respect to the morphology. Moreover, although the behavior of soft tissues is highly nonlinear, only linear normalization protocols have been defined so far for the disc stiffness. Thus, the aim of this work is to propose a nonlinear normalization based on finite elements (FE) simulations and evaluate its impact on the stiffness of human anatomical specimens of lumbar IVD. First, a parameter study involving simulations of biomechanical tests (compression, flexion/extension, bilateral torsion and bending) on 20 FE models of IVDs with various dimensions was carried out to evaluate the effect of the disc's geometry on its compliance and establish stiffness/morphology relations necessary to the nonlinear normalization. The computed stiffness was then normalized by height (H), cross-sectional area (CSA), polar moment of inertia (J) or moments of inertia (Ixx, Iyy) to quantify the effect of both linear and nonlinear normalizations. In the second part of the study, T1-weighted MRI images were acquired to determine H, CSA, J, Ixx and Iyy of 14 human lumbar IVDs. Based on the measured morphology and pre-established relation with stiffness, linear and nonlinear normalization routines were then applied to the compliance of the specimens for each quasi-static biomechanical test. The variability of the stiffness prior to and after normalization was assessed via coefficient of variation (CV). The FE study confirmed that larger and thinner IVDs were stiffer while the normalization strongly attenuated the effect of the disc geometry on its stiffness. Yet, notwithstanding the results of the FE study, the experimental stiffness showed consistently higher CV after normalization. Assuming that geometry and material properties affect the mechanical response, they can also compensate for one another. Therefore, the larger CV after normalization can be interpreted as a strong variability of the material properties, previously hidden by the geometry's own influence. In conclusion, a new normalization protocol for the intervertebral disc stiffness in compression, flexion, extension, bilateral torsion and bending was proposed, with the possible use of MRI and FE to acquire the discs' anatomy and determine the nonlinear relations between stiffness and morphology. Such protocol may be useful to relate the disc's mechanical properties to its degree of degeneration.

Baclofen and gamma-hydroxybutyrate differentially altered behavior, EEG activity and sleep in rats.

OBJECTIVES Animal and human studies have shown that sleep may have an impact on functional recovery after brain damage. Baclofen (Bac) and gamma-hydroxybutyrate (GHB) have been shown to induce physiological sleep in humans, however, their effects in rodents are unclear. The aim of this study is to characterize sleep and electroencelphalogram (EEG) after Bac and GHB administration in rats. We hypothesized that both drugs would induce physiological sleep. METHODS Adult male Sprague-Dawley rats were implanted with EEG/electromyogram (EMG) electrodes for sleep recordings. Bac (10 or 20 mg/kg), GHB (150 or 300 mg/kg) or saline were injected 1 h after light and dark onset to evaluate time of day effect of the drugs. Vigilance states and EEG spectra were quantified. RESULTS Bac and GHB induced a non-physiological state characterized by atypical behavior and an abnormal EEG pattern. After termination of this state, Bac was found to increase the duration of non-rapid eye movement (NREM) and rapid eye movement (REM) sleep (∼90 and 10 min, respectively), reduce sleep fragmentation and affect NREM sleep episode frequency and duration (p<0.05). GHB had no major effect on vigilance states. Bac drastically increased EEG power density in NREM sleep in the frequencies 1.5-6.5 and 9.5-21.5 Hz compared to saline (p<0.05), while GHB enhanced power in the 1-5-Hz frequency band and reduced it in the 7-9-Hz band. Slow-wave activity in NREM sleep was enhanced 1.5-3-fold during the first 1-2 h following termination of the non-physiological state. The magnitude of drug effects was stronger during the dark phase. CONCLUSION While both Bac and GHB induced a non-physiological resting state, only Bac facilitated and consolidated sleep, and promoted EEG delta oscillations thereafter. Hence, Bac can be considered a sleep-promoting drug and its effects on functional recovery after stroke can be evaluated both in humans and rats.

Menstrual cycle phase affects discrimination of infant cuteness

Recent studies have shown that women are more sensitive than men to subtle cuteness differences in infant faces. It has been suggested that raised levels in estradiol and progesterone may be responsible for this advantage. We compared young women's sensitivity to computer-manipulated baby faces varying in cuteness. Thirty-six women were tested once during ovulation and once during the luteal phase of their menstrual cycle. In a two alternative forced-choice experiment, participants chose the baby which they thought was cuter (Task 1), younger (Task 2), or the baby that they would prefer to babysit (Task 3). Saliva samples to assess levels of estradiol, progesterone and testosterone were collected at each test session. During ovulation, women were more likely to choose the cuter baby than during the luteal phase, in all three tasks. These results suggest that cuteness discrimination may be driven by cyclic hormonal shifts. However none of the measured hormones were related to increased cuteness sensitivity. We speculate that other hormones than the ones measured here might be responsible for the increased sensitivity to subtle cuteness differences during ovulation.