Young Disk Research Articles

Resolving the disc–halo degeneracy – I: a look at NGC 628

The decomposition of the rotation curve of galaxies into contribution from the disc and dark halo remains uncertain and depends on the adopted mass to light ratio (M/L) of the disc. Given the vertical velocity dispersion of stars and disc scale height, the disc surface mass density and hence the M/L can be estimated. We address a conceptual problem with previous measurements of the scale height and dispersion. When using this method, the dispersion and scale height must refer to the same population of stars. The scale height is obtained from near-IR studies of edge-on galaxies and is weighted towards older kinematically hotter stars, whereas the dispersion obtained from integrated light in the optical bands includes stars of all ages. We aim to extract the dispersion for the hotter stars, so that it can then be used with the correct scale height to obtain the disc surface mass density. We use a sample of planetary nebulae (PNe) as dynamical tracers in the face-on galaxy NGC 628. We extract two different dispersions from its velocity histogram -- representing the older and younger PNe. We also present complementary stellar absorption spectra in the inner regions of this galaxy and use a direct pixel fitting technique to extract the two components. Our analysis concludes that previous studies, which do not take account of the young disc, underestimate the disc surface mass density by a factor of ~ 2. This is sufficient to make a maximal disc for NGC 628 appear like a submaximal disc.

The Complex Morphology of the Young Disk MWC 758: Spirals and Dust Clumps around a Large Cavity

Abstract We present Atacama Large Millimeter Array observations at an angular resolution of 0.″1–0.″2 of the disk surrounding the young Herbig Ae star MWC 758. The data consist of images of the dust continuum emission recorded at 0.88 millimeter, as well as images of the 13CO and C18O J = 3–2 emission lines. The dust continuum emission is characterized by a large cavity of roughly 40 au in radius which might contain a mildly inner warped disk. The outer disk features two bright emission clumps at radii of ∼47 and 82 au that present azimuthal extensions and form a double-ring structure. The comparison with radiative transfer models indicates that these two maxima of emission correspond to local increases in the dust surface density of about a factor 2.5 and 6.5 for the south and north clumps, respectively. The optically thick 13CO peak emission, which traces the temperature, and the dust continuum emission, which probes the disk midplane, additionally reveal two spirals previously detected in near-IR at the disk surface. The spirals seen in the dust continuum emission present, however, a slight shift of a few au toward larger radii and one of the spirals crosses the south dust clump. Finally, we present different scenarios to explain the complex structure of the disk.

Investigating the Early Evolution of Planetary Systems with ALMA and the Next Generation Very Large Array

Abstract We investigate the potential of the Atacama Large Millimeter/submillimeter Array (ALMA) and the Next Generation Very Large Array (ngVLA) to observe substructures in nearby young disks which are due to the gravitational interaction between disk material and planets close to the central star. We simulate the gas and dust dynamics in the disk using the LA-COMPASS hydrodynamical code. We generate synthetic images for the dust continuum emission at submillimeter to centimeter wavelengths and simulate ALMA and ngVLA observations. We explore the parameter space of some of the main disk and planet properties that would produce substructures that can be visible with ALMA and the ngVLA. We find that ngVLA observations with an angular resolution of 5 milliarcsec at 3 mm can reveal and characterize gaps and azimuthal asymmetries in disks hosting planets with masses down to from a solar-like star in the closest star-forming regions, whereas ALMA can detect gaps down to planetary masses of at 5 au. Gaps opened by super-Earth planets with masses are detectable by the ngVLA in the case of disks with low viscosity ( ) and low pressure scale height (h ≈ 0.025 au at 5 au). The ngVLA can measure the proper motion of azimuthal asymmetric structures associated with the disk–planet interaction as well as possible circumplanetary disks on timescales as short as one to a few weeks for planets at 1–5 au from the star.

On the diversity and statistical properties of protostellar discs

AbstractWe present results from the first population synthesis study of protostellar discs. We analyse the evolution and properties of a large sample of protostellar discs formed in a radiation hydrodynamical simulation of star cluster formation. Due to the chaotic nature of the star formation process, we find an enormous diversity of young protostellar discs, including misaligned discs, and discs whose orientations vary with time. Star–disc interactions truncate discs and produce multiple systems. Discs may be destroyed in dynamical encounters and/or through ram-pressure stripping, but reform by later gas accretion. We quantify the distributions of disc mass and radii for protostellar ages up to ≈105 yr. For low-mass protostars, disc masses tend to increase with both age and protostellar mass. Disc radii range from of order 10 to a few hundred au, grow in size on time-scales ≲ 104 yr, and are smaller around lower mass protostars. The radial surface density profiles of isolated protostellar discs are flatter than the minimum mass solar nebula model, typically scaling as Σ ∝ r−1. Disc to protostar mass ratios rarely exceed two, with a typical range of Md/M* = 0.1–1 to ages ≲ 104 yr and decreasing thereafter. We quantify the relative orientation angles of circumstellar discs and the orbit of bound pairs of protostars, finding a preference for alignment that strengths with decreasing separation. We also investigate how the orientations of the outer parts of discs differ from the protostellar and inner disc spins for isolated protostars and pairs.

Young Stellar Variability of GM Cephei by Circumstellar Dust Clumps

UX Orionis stars are a sub-type of Herbig Ae/be or T Tauri stars exhibiting sporadic extinction of stellar light due to circumstellar dust obscuration. GM Cep is such an UX Orionis star in the young (∼ 4 Myr) open cluster Trumpler 37 at ∼ 900 pc, showing a prominent infrared access, H-alpha emission, and flare activity. Our multi-color photometric monitoring from 2009 to 2016 showed (i) sporadic brightening on a time scale of days due to young stellar accretion, (ii) cyclic, but not strictly periodical, occultation events, each lasting for a couple months, with a probable recurrence time of about two years, (iii) normal dust reddening as the star became redder when dimmer, (iv) the unusual ”blueing” phenomena near the brightness minima, during which the star appeared bluer when dimmer, and (v) a noticeable polarization, from 3 to 9 percent in g’, r’, and i’ -bands. The occultation events may be caused by dust clumps, signifying the density inhomogeneity in a young stellar disk from grain coagulation to planetesimal formation. The level of polarization was anti-correlated with the brightness in the bright state, when the dust clump backscattered stellar light. We discussed two potential hypotheses: orbiting dust clumps versus dust clumps along a spiral arm structure.

Light-curve Modulation of Low-mass Stars in K2. I. Identification of 481 Fast Rotators in the Solar Neighborhood

Abstract The K2 mission is targeting large numbers of nearby (d < 100 pc) GKM dwarfs selected from the SUPERBLINK proper motion survey (μ > 40 mas yr−1, V < 20). Additionally, the mission is targeting low-mass, high proper motion stars associated with the local (d < 500 pc) Galactic halo population also selected from SUPERBLINK. K2 campaigns 0 through 8 monitored a total of 26,518 of these cool main-sequence stars. We used the auto-correlation function to search for fast rotators by identifying short-period photometric modulations in the K2 light curves. We identified 481 candidate fast rotators with rotation periods <4 days that show light-curve modulations consistent with starspots. Their kinematics show low average transverse velocities, suggesting that they are part of the young disk population. A subset (13) of the fast rotators is found among those targets with colors and kinematics consistent with the local Galactic halo population and may represent stars spun up by tidal interactions in close binary systems. We further demonstrate that the M dwarf fast rotators selected from the K2 light curves are significantly more likely to have UV excess and discuss the potential of the K2 mission to identify new nearby young GKM dwarfs on the basis of their fast rotation rates. Finally, we discuss the possible use of local halo stars as fiducial, non-variable sources in the Kepler fields.

On the inner disc structure of MWC480: evidence for asymmetries?

Studying the physical conditions structuring the young circumstellar disks is required for understanding the onset of planet formation. Of particular interest is the protoplanetary disk surrounding the Herbig star MWC480. The structure and properties of the circumstellar disk of MWC480 are studied by infrared interferometry and interpreted from a modeling approach. New observations are driving this study, in particular some recent Very Large Telescope Interferometer (VLTI)/MIDI data acquired in December 2013. Our one-component disk model could not reproduce simultaneously all our data: the Spectral Energy Distribution, the near-infrared Keck Interferometer data and the mid-infrared data obtained with the MIDI instrument. In order to explain all measurements, one possibility is to add an asymmetry in our one-component disk model with the assumption that the structure of the disk of MWC480 has not varied with time. Several scenarios are tested, and the one considering the presence of an azimuthal bright feature in the inner component of the disk model provides a better fit of the data. (In this study, we assumed that the size of the outer disk of MWC480 to be 20 au since most of the near and mid-IR emissions come from below 20 au. In our previous study (Jamialahmadi et al. 2015), we adopted an outer radius of 80 au, which is consistent with the value found by previous studies based on mm observations).

Twelve Years of Spectroscopic Monitoring in the Galactic Center: The Closest Look at S-stars near the Black Hole

Abstract We study the young S-stars within a distance of 0.04 pc from the supermassive black hole in the center of our Galaxy. Given how inhospitable the region is for star formation, their presence is more puzzling the younger we estimate their ages. In this study, we analyze the result of 12 years of high-resolution spectroscopy within the central arcsecond of the Galactic Center (GC). By co-adding between 55 and 105 hr of spectra we have obtained high signal-to-noise H- and K-band spectra of eight stars orbiting the central supermassive black hole. Using deep H-band spectra, we show that these stars must be high surface gravity (dwarf) stars. We compare these deep spectra to detailed model atmospheres and stellar evolution models to infer the stellar parameters. Our analysis reveals an effective temperature of 21,000–28,500 K, a rotational velocity of 60–170 km s−1, and a surface gravity of 4.1–4.2. These parameters imply a spectral type of B0–B3V for these stars. The inferred masses lie within 8–14 . We derive an age of Myr for the star S2, which is compatible with the age of the clockwise-rotating young stellar disk in the GC. We estimate the ages of all other studied S-stars to be less than 15 Myr, which is compatible with the age of S2 within the uncertainties. The relatively low ages for these S-stars favor a scenario in which the stars formed in a local disk rather than a field binary-disruption scenario that occurred over a longer period of time.

Linear growth of streaming instability in pressure bumps

Streaming instability is a powerful mechanism which concentrates dust grains in pro- toplanetary discs, eventually up to the stage where they collapse gravitationally and form planetesimals. Previous studies inferred that it should be ineffective in viscous discs, too efficient in inviscid discs, and may not operate in local pressure maxima where solids accumulate. From a linear analysis of stability, we show that streaming instability behaves differently inside local pressure maxima. Under the action of the strong differential advection imposed by the bump, a novel unstable mode develops and grows even when gas viscosity is large. Hence, pressure bumps are found to be the only places where streaming instability occurs in viscous discs. This offers a promising way to conciliate models of planet formation with recent observations of young discs.

Protoplanetary disc ‘isochrones’ and the evolution of discs in the M˙-Md plane

In this paper, we compare simple viscous diffusion models for the disc evolution with the results of recent surveys of the properties of young protoplanetary discs. We introduce the useful concept of `disc isochrones' in the accretion rate - disc mass plane and explore a set of Montecarlo realization of disc initial conditions. We find that such simple viscous models can provide a remarkable agreement with the available data in the Lupus star forming region, with the key requirement that the average viscous evolutionary timescale of the discs is comparable to the cluster age. Our models produce naturally a correlation between mass accretion rate and disc mass that is shallower than linear, contrary to previous results and in agreement with observations. We also predict that a linear correlation, with a tighter scatter, should be found for more evolved disc populations. Finally, we find that such viscous models can reproduce the observations in the Lupus region only in the assumption that the efficiency of angular momentum transport is a growing function of radius, thus putting interesting constraints on the nature of the microscopic processes that lead to disc accretion.

Magnetic activity and radial velocity filtering of young Suns: The weak-line T Tauri stars Par 1379 and Par 2244

We report the results of our spectropolarimetric monitoring of the weak-line T-Tauri stars (wTTSs) Par 1379 and Par 2244, within the MaTYSSE (Magnetic Topologies of Young Stars and the Survival of close-in giant Exoplanets) programme. Both stars are of a similar mass (1.6 and 1.8 M$_{\odot}$) and age (1.8 and 1.1 Myr), with Par 1379 hosting an evolved low-mass dusty circumstellar disc, and with Par 2244 showing evidence of a young debris disc. We detect profile distortions and Zeeman signatures in the unpolarized and circularly polarized lines for each star, and have modelled their rotational modulation using tomographic imaging, yielding brightness and magnetic maps. We find that Par 1379 harbours a weak (250 G), mostly poloidal field tilted $65^{\circ}$ from the rotation axis. In contrast, Par 2244 hosts a stronger field (860 G) split 3:2 between poloidal and toroidal components, with most of the energy in higher order modes, and with the poloidal component tilted $45^{\circ}$ from the rotation axis. Compared to the lower mass wTTSs, V819 Tau and V830 Tau, Par 2244 has a similar field strength, but is much more complex, whereas the much less complex field of Par 1379 is also much weaker than any other mapped wTTS. We find moderate surface differential rotation of $1.4\times$ and $1.8\times$ smaller than Solar, for Par 1379 and Par 2244, respectively. Using our tomographic maps to predict the activity related radial velocity (RV) jitter, and filter it from the RV curves, we find RV residuals with dispersions of 0.017 kms$^{-1}$ and 0.086 kms$^{-1}$ for Par 1379 and Par 2244, respectively. We find no evidence for close-in giant planets around either star, with $3\sigma$ upper limits of 0.56 and 3.54 M$_{\text{jup}}$ (at an orbital distance of 0.1 au).

The structure of young embedded protostellar discs

Young protostellar discs provide the initial conditions for planet formation. The properties of these discs may be different from those of late-phase (T Tauri) discs due to continuing infall from the envelope and protostellar variability resulting from irregular gas accretion. We use a set of hydrodynamic simulations to determine the structure of discs forming in collapsing molecular clouds. We examine how radiative feedback from the host protostar affects the disc properties by examining three regimes: without radiative feedback, with continuous radiative feedback and with episodic feedback, similar to FU Ori-type outbursts. We find that the radial surface density and temperature profiles vary significantly as the disc accretes gas from the infalling envelope. These profiles are sensitive to the presence of spiral structure, induced by gravitational instabilities, and the radiative feedback provided by the protostar, especially in the case when the feedback is episodic. We also investigate whether mass estimates from position-velocity (PV) diagrams are accurate for early-phase discs. We find that the protostellar system mass (i.e. the mass of the protostar and its disc) is underestimated by up to 20%, due to the impact of an enhanced radial pressure gradient on the gas. The mass of early-phase discs is a significant fraction of the mass of the protostar, so position-velocity diagrams cannot accurately provide the mass of the protostar alone. The enhanced radial pressure gradient expected in young discs may lead to an increased rate of dust depletion due to gas drag, and therefore to a reduced dust-to-gas ratio.

Scattering-produced (sub)millimetre polarization in inclined discs: optical depth effects, near–far side asymmetry and dust settling

Disk polarization at (sub)millimeter wavelengths is being revolutionized by ALMA observationally, but its origin remains uncertain. Dust scattering was recently recognized as a potential contributor to polarization, although its basic properties have yet to be thoroughly explored. Here, we quantify the effects of optical depth on the scattering-induced polarization in inclined disks through a combination of analytical illustration, approximate semi-analytical modeling using formal solution to the radiative transfer equation, and Monte Carlo simulations. We find that the near-side of the disk is significantly brighter in polarized intensity than the far-side, provided that the disk is optically thick and that the scattering grains have yet to settle to the midplane. This asymmetry is the consequence of a simple geometric effect: the near-side of the disk surface is viewed more edge-on than the far-side. It is a robust signature that may be used to distinguish the scattering-induced polarization from that by other mechanisms, such as aligned grains. The asymmetry is weaker for a geometrically thinner dust disk. As such, it opens an exciting new window on dust settling. We find anecdotal evidence from dust continuum imaging of edge-on disks that large grains are not yet settled in the youngest (Class 0) disks, but become more so in older disks. This trend is corroborated by the polarization data in inclined disks showing that younger disks have more pronounced near-far side asymmetry and thus less grain settling. If confirmed, the trend would have far-reaching implications for grain evolution and, ultimately, the formation of planetesimals and planets.

Simulations of small solid accretion on to planetesimals in the presence of gas

The growth and migration of planetesimals in a young protoplanetary disc are fundamental to planet formation. In all models of early growth, there are several processes that can inhibit grains from reaching larger sizes. Nevertheless, observations suggest that growth of planetesimals must be rapid. If a small number of 100 km sized planetesimals do manage to form in the disc, then gas drag effects could enable them to efficiently accrete small solids from beyond their gravitationally focused cross-section. This gas drag-enhanced accretion can allow planetesimals to grow at rapid rates, in principle. We present self-consistent hydrodynamics simulations with direct particle integration and gas drag coupling to estimate the rate of planetesimal growth due to pebble accretion. Wind tunnel simulations are used to explore a range of particle sizes and disc conditions. We also explore analytic estimates of planetesimal growth and numerically integrate planetesimal drift due to the accretion of small solids. Our results show that, for almost every case that we consider, there is a clearly preferred particle size for accretion that depends on the properties of the accreting planetesimal and the local disc conditions. For solids much smaller than the preferred particle size, accretion rates are significantly reduced as the particles are entrained in the gas and flow around the planetesimal. Solids much larger than the preferred size accrete at rates consistent with gravitational focusing. Our analytic estimates for pebble accretion highlight the timescales that are needed for the growth of large objects under different disc conditions and initial planetesimal sizes.

The First Scattered-light Image of the Debris Disk around the Sco–Cen Target HD 129590

Abstract We present the first scattered-light image of the debris disk around HD 129590, a ∼1.3 M ⊙ G1V member of the Scorpius–Centaurus association with an age of ∼10–16 Myr. The debris disk is imaged with the high contrast imaging instrument SPHERE at the Very Large Telescope, and is revealed by both the IRDIS and IFS subsytems, operating in the H and YJ bands respectively. The disk has a high infrared luminosity of ∼ 5 × 10−3, and has been resolved in other studies using ALMA. We detect a nearly edge-on ring, with evidence of an inner clearing. We fit the debris disk using a model characterized by a single bright ring, with radius ∼60–70 au, in broad agreement with previous analyses of the target SED. The disk is vertically thin, and has an inclination angle of ∼75°. Along with other previously imaged edge-on disks in the Sco–Cen association such as HD 110058, HD 115600, and HD 111520, this disk image will allow of the structure and morphology of very young debris disks, shortly after the epoch of planet formation has ceased.

Gravitoturbulence and the excitation of small-scale parametric instability in astrophysical discs

Young protoplanetary discs and the outer radii of active galactic nucleii may be subject to gravitational instability and, as a consequence, fall into a `gravitoturbulent' state. While in this state, appreciable angular momentum can be transported. Alternatively, the gas may collapse into bound clumps, the progenitors of planets or stars. In this paper, we numerically characterize the properties of 3D gravitoturbulence, focussing especially on its dependence on numerical parameters (resolution, domain size) and its excitation of small-scale dynamics. Via a survey of vertically stratified shearing box simulations with PLUTO and RODEO, we find (a) evidence that certain gravitoturbulent properties are independent of horizontal box size only when the box is larger than $\simeq 40 H_0$, where $H_0$ is the height scale, (b) at high resolution, small-scale isotropic turbulence appears off the midplane around $z\simeq 0.5 -1 H_0$, and (c) this small-scale dynamics results from a parametric instability, involving the coupling of inertial waves with a large-scale axisymmetric epicyclic mode. This mode oscillates at a frequency close to $\Omega$ and is naturally excited by gravito-turbulence, via a nonlinear process to be determined. The small-scale turbulence we uncover has potential implications for a wide range of disc physics, e.g. turbulent saturation levels, fragmentation, turbulent mixing, and dust settling.

A Resolved and Asymmetric Ring of PAHs within the Young Circumstellar Disk of IRS 48

Abstract For one decade, the spectral type and age of the ρ Oph object IRS-48 were subject to debate and mystery. Modeling its disk with mid-infrared to millimeter observations led to various explanations to account for the complex intricacy of dust holes and gas-depleted regions. We present multi-epoch high-angular-resolution interferometric near-infrared data of spatially resolved emissions in the first 15 au of IRS-48, known to have very strong polycyclic aromatic hydrocarbon (PAH) emissions within this dust-depleted region. We make use of new Sparse-Aperture-Masking data to instruct a revised radiative-transfer model, where spectral energy distribution fluxes and interferometry are jointly fitted. Neutral and ionized PAH, very small grains (VSG), and classical silicates are incorporated into the model; new stellar parameters and extinction laws are explored. A bright (42 L ⊙) and hence large (2.5 R ⊙) central star with A v = 12.5 mag and R v = 6.5 requires less near-infrared excess: the inner-most disk at ≈1 au is incompatible with the interferometric data. The revised stellar parameters place this system on a 4 Myr evolutionary track, four times younger than the previous estimations, which is in better agreement with the surrounding ρ Oph region and disk-lifetime observations. The disk-structure solution converges to a classical-grain outer disk from 55 au combined with an unsettled and fully resolved VSG and PAH ring, between 11 and 26 au. We find two overluminosities in the PAH ring at color-temperatures consistent with the radiative transfer simulations; one follows a Keplerian circular orbit at 14 au. We show a depletion of a factor of ≈5 of classical dust grains up to 0.3 mm compared to very small particles: the IRS-48 disk is nearly void of dust grains in the first 55 au. A 3.5 M Jup planet on a 40 au orbit can qualitatively explain the new disk structure.

Planetesimal formation near the snowline: in or out?

The formation of planetesimals in protoplanetary disks is not well-understood. Streaming instability is a promising mechanism to directly form planetesimals from pebble-sized particles, provided a high enough solids-to-gas ratio. However, local enhancements of the solids-to-gas ratio are difficult to realize in a smooth disk, which motivates the consideration of special disk locations such as the snowline. In this article we investigate the viability of planetesimal formation by streaming instability near the snowline due to water diffusion and condensation. We adopt a viscous disk model, and numerically solve the transport equations for vapor and solids on a cylindrical, 1D grid. We take into account radial drift of solids, gas accretion on to the central star, and turbulent diffusion. We study the importance of the back-reaction of solids on the gas and of the radial variation of the mean molecular weight of the gas. We find that water diffusion and condensation can locally enhance the ice surface density by a factor 3--5 outside the snowline. Assuming that icy pebbles contain many micron-sized silicate grains that are released during evaporation, the enhancement is increased by another factor $\sim$$2$. In this `many-seeds' model, the solids-to-gas ratio interior to the snowline is enhanced as well, but not as much as just outside the snowline. In the context of a viscous disk, the diffusion-condensation mechanism is most effective for high values of the turbulence parameter $\alpha$ ($10^{-3}$--$10^{-2}$). Therefore, assuming young disks are more vigorously turbulent than older disks, planetesimals near the snowline can form in an early stage of the disk. In highly turbulent disks, tens of Earth masses can be stored in an annulus outside the snowline, which can be identified with recent ALMA observations.

Saving Super-Earths: Interplay between Pebble Accretion and Type I Migration

Abstract Overcoming type I migration and preventing low-mass planets from spiralling into the central star is a long-studied topic. It is well known that outward migration is possible in viscously heated disks relatively close to the central star because the entropy gradient can be sufficiently steep for the positive corotation torque to overcome the negative Lindblad torque. Yet efficiently trapping planets in this region remains elusive. Here we study disk conditions that yield outward migration for low-mass planets under specific planet migration prescriptions. In a steady-state disk model with a constant α-viscosity, outward migration is only possible when the negative temperature gradient exceeds ∼0.87. We derive an implicit relation for the highest mass at which outward migration is possible as a function of viscosity and disk scale height. We apply these criteria, using a simple power-law disk model, to planets that have reached their pebble isolation mass after an episode of rapid accretion. It is possible to trap planets with the pebble isolation mass farther than the inner edge of the disk provided that α crit ≳ 0.004 for disks older than 1 Myr. In very young disks, the high temperature causes the planets to grow to masses exceeding the maximum for outward migration. As the disk evolves, these more massive planets often reach the central star, generally only toward the end of the disk lifetime. Saving super-Earths is therefore a delicate interplay between disk viscosity, the opacity profile, and the temperature gradient in the viscously heated inner disk.

EPR Studies on the Properties of Model Photoreceptor Membranes Made of Natural and Synthetic Lipids

The membranes of retina photoreceptors have unique lipid composition. They contain a high concentration of polyunsaturated docosahexaenoic acid, with six double bonds, and are enriched in phosphatidylethanolamines. Based on their phospholipid composition and cholesterol content, membranes of photoreceptors can be divided into three types: plasma membrane, young disks membranes, and old disks membranes. High amount of docosahexaenoic acid, abundant illumination, and high respiratory demands make these membranes sensitive to oxidative stress and lipid peroxidation. Human retinas are not easily available for research, therefore most research is done on bovine retinas. However, to follow, in a controlled manner, the changes in membrane properties caused by different factors it seems advisable to apply carefully prepared models of photoreceptor membranes. Using synthetic lipids we prepared liposome models of three types of photoreceptor membranes, and by means of electron paramagnetic resonance spectroscopy and spin labeling technique we compared polarity and fluidity of those model membranes with the properties of membranes consisting of natural lipids extracted from photoreceptor outer segments of bovine retinas. Additionally, we studied the effect of oxidation on the membrane properties in the presence and in the absence of zeaxanthin, which is an antioxidant naturally present in the human retina. The results show that there are significant differences in polarity and fluidity between all investigated membranes, which reflect differences in their lipid composition. The properties of the membranes made of natural photoreceptor outer segment lipids are most similar to the ones of the models of old disks membranes. Oxidation did not change the membrane properties significantly; however, a slight ordering effect was observed in liposomes made of natural photoreceptor outer segment lipids and in the model of old disks membranes. Zeaxanthin affected polarity and fluidity mostly in the model of old disks membranes. The results show that by careful selection and appropriate proportions of lipid mixtures, it is possible to obtain synthetic membranes of the properties similar to the natural ones.