Cone model-based SEP event calculations for applications to multipoint observations

The problem of modeling solar energetic particle (SEP) events is important to both space weather research and forecasting, and yet it has seen relatively little progress. Most important SEP events are associated with coronal mass ejections (CMEs) that drive coronal and interplanetary shocks. These shocks can continuously produce accelerated particles from the ambient medium to well beyond 1 AU. This paper describes an effort to model real SEP events using a Center for Integrated Space weather Modeling (CISM) MHD solar wind simulation including a cone model of CMEs to initiate the related shocks. In addition to providing observation-inspired shock geometry and characteristics, this MHD simulation describes the time-dependent observer field line connections to the shock source. As a first approximation, we assume a shock jump-parameterized source strength and spectrum, and that scatter-free transport occurs outside of the shock source, thus emphasizing the role the shock evolution plays in determining the modeled SEP event profile. Three halo CME events on May 12, 1997, November 4, 1997 and December 13, 2006 are used to test the modeling approach. While challenges arise in the identification and characterization of the shocks in the MHD model results, this approach illustrates the importance to SEP event modeling of globally simulating the underlying heliospheric event. The results also suggest the potential utility of such a model for forcasting and for interpretation of separated multipoint measurements such as those expected from the STEREO mission.

The formation of large-scale current sheets within magnetic clouds

Magnetic clouds are a class of interplanetary coronal mass ejections (CME) predominantly characterised by a smooth rotation in the magnetic field direction, indicative of a magnetic flux rope structure. Many magnetic clouds, however, also contain sharp discontinuities within the smoothly varying magnetic field, suggestive of narrow current sheets. In this study we present observations and modelling of magnetic clouds with strong current sheet signatures close to the centre of the apparent flux rope structure. Using an analytical magnetic flux rope model, we demonstrate how such current sheets can form as a result of a cloud’s kinematic propagation from the Sun to the Earth, without any external forces or influences. This model is shown to match observations of four particular magnetic clouds remarkably well. The model predicts that current sheet intensity increases for increasing CME angular extent and decreasing CME radial expansion speed. Assuming such current sheets facilitate magnetic reconnection, the process of current sheet formation could ultimately lead a single flux rope becoming fragmented into multiple flux ropes. This change in topology has consequences for magnetic clouds as barriers to energetic particle propagation.

An overview of Earth’s global electric circuit and atmospheric conductivity

The Earth’s global atmospheric electric circuit depends on the upper and lower atmospheric boundaries formed by the ionosphere and the planetary surface. Thunderstorms and electrified rain clouds drive a DC current (∼1 kA) around the circuit, with the current carried by molecular cluster ions; lightning phenomena drive the AC global circuit. The Earth’s near-surface conductivity ranges from 10−7 S m−1 (for poorly conducting rocks) to 10−2 S m−1 (for clay or wet limestone), with a mean value of 3.2 S m−1 for the ocean. Air conductivity inside a thundercloud, and in fair weather regions, depends on location (especially geomagnetic latitude), aerosol pollution and height, and varies from ∼10−14 S m−1 just above the surface to 10−7 S m−1 in the ionosphere at ∼80 km altitude. Ionospheric conductivity is a tensor quantity due to the geomagnetic field, and is determined by parameters such as electron density and electron–neutral particle collision frequency. In the current source regions, point discharge (coronal) currents play an important role below electrified clouds; the solar wind-magnetosphere dynamo and the unipolar dynamo due to the terrestrial rotating dipole moment also apply atmospheric potential differences. Detailed measurements made near the Earth’s surface show that Ohm’s law relates the vertical electric field and current density to air conductivity. Stratospheric balloon measurements launched from Antarctica confirm that the downward current density is ∼1 pA m−2 under fair weather conditions. Fortuitously, a Solar Energetic Particle (SEP) event arrived at Earth during one such balloon flight, changing the observed atmospheric conductivity and electric fields markedly. Recent modelling considers lightning discharge effects on the ionosphere’s electric potential (∼+250 kV with respect to the Earth’s surface) and hence on the fair weather potential gradient (typically ∼130 V m−1 close to the Earth’s surface. We conclude that cloud-to-ground (CG) lightning discharges make only a small contribution to the ionospheric potential, and that sprites (namely, upward lightning above energetic thunderstorms) only affect the global circuit in a miniscule way. We also investigate the effects of mesoscale convective systems on the global circuit.

Solar induced climate effects

Solar electromagnetic radiation powers Earth’s climate system and, consequently, it is often naively assumed that changes in this solar output must be responsible for changes in Earth’s climate. However, the Sun is close to a blackbody radiator and so emits according to its surface temperature and the huge thermal time constant of the outer part of the Sun limits the variability in surface temperature and hence output. As a result, on all timescales of interest, changes in total power output are limited to small changes in effective surface temperature (associated with magnetic fields) and potential, although as yet undetected, solar radius variations. Larger variations are seen in the UV part of the spectrum which is emitted from the lower solar atmosphere (the chromosphere) and which influences Earth’s stratosphere. There is interest in“top-down” mechanisms whereby solar UV irradiance modulates stratospheric temperatures and winds which, in turn, may influence the underlying troposphere where Earth’s climate and weather reside. This contrasts with “bottom-up” effects in which the small total solar irradiance (dominated by the visible and near-IR) variations cause surface temperature changes which drive atmospheric circulations. In addition to these electromagnetic outputs, the Sun modulates energetic particle fluxes incident on the Earth. Solar Energetic Particles (SEP) are emitted by solar flares and from the shock fronts ahead of supersonic (and super-Alfvenic) ejections of material from the solar atmosphere. These SEPs enhance the destruction of polar stratospheric ozone which could be an additional form of top-down climate forcing. Even more energetic are Galactic Cosmic Rays (GCRs). These particles are not generated by the Sun, rather they originate at the shock fronts emanating from violent galactic events such as supernovae explosions; however, the expansion of the solar magnetic field into interplanetary space means that the Sun modulates the number of GCRs reaching Earth. These play a key role in enabling Earth’s global electric (thunderstorm) circuit and it has been proposed that they also modulate the formation of clouds. Both electromagnetic and corpuscular solar effects are known to vary over the solar magnetic cycle which is typically between 10 and 14 yrs in length (with an average close to 11 yrs). The solar magnetic field polarity at any one phase of one of these activity cycles is opposite to that at the same phase of the next cycle and this influences some phenomena, for example GCRs, which therefore show a 22 yr (“Hale”) cycle on average. Other phenomena, such as irradiance modulation, do not depend on the polarity of the magnetic field and so show only the basic 11-yr activity cycle. However, any effects on climate are much more significant for solar drifts over centennial timescales. This chapter discusses and evaluates potential effects on Earth’s climate system of variations in these solar inputs. Because of the great variety of proposed mechanisms, the wide range of timescales studied (from days to millennia) and the many debates (often triggered by the application of inadequate statistical methods), the literature on this subject is vast, complex, divergent and rapidly changing: consequently the number of references cited in this review is very large (yet still only a small fraction of the total).

An integrative model of amino acid metabolism in the liver of the lactating dairy cow

The objective of this work was to construct a dynamic model of hepatic amino acid metabolism in the lactating dairy cow that could be parameterized using net flow data from in vivo experiments. The model considers 22 amino acids, ammonia, urea, and 13 energetic metabolites, and was parameterized using a steady-state balance model and two in vivo, net flow experiments conducted with mid-lactation dairy cows. Extracellular flows were derived directly from the observed data. An optimization routine was used to derive nine intracellular flows. The resulting dynamic model was found to be stable across a range of inputs suggesting that it can be perturbed and applied to other physiological states. Although nitrogen was generally in balance, leucine was in slight deficit compared to predicted needs for export protein synthesis, suggesting that an alternative source of leucine (e.g. peptides) was utilized. Simulations of varying glucagon concentrations indicated that an additional 5 mol/d of glucose could be synthesized at the reference substrate concentrations and blood flows. The increased glucose production was supported by increased removal from blood of lactate, glutamate, aspartate, alanine, asparagine, and glutamine. As glucose Output increased, ketone body and acetate release increased while CO2 release declined. The pattern of amino acids appearing in hepatic vein blood was affected by changes in amino acid concentration in portal vein blood, portal blood flow rate and glucagon concentration, with methionine and phenylalanine being the most affected of essential amino acids. Experimental evidence is insufficient to determine whether essential amino acids are affected by varying gluconeogenic demands. (C) 2004 Published by Elsevier Ltd.

Effect of replacing calcium salts of palm oil distillate with extruded linseeds on milk fatty acid composition in Jersey and Holstein cows

Clinical and biomedical studies have provided evidence for the critical role of n-3 fatty acids on the reduction of chronic disease risk in humans, including cardiovascular disease. In the current experiment, the potential to enhance milk n-3 content in two breeds with inherent genetic differences in mammary lipogenesis and de novo fatty acid synthesis was examined using extruded linseeds. Six lactating cows (three Holstein and three Jersey) were used in a two-treatment switchback design with 3 × 21-day experimental periods to evaluate the effect of iso-energetic replacement of calcium salts of palm oil distillate (CPO) in the diet (34 g/kg dry matter (DM)) with 100 g/kg DM extruded linseeds (LIN). For both breeds, replacing CPO with LIN had no effect (P > 0.05) on DM intake or milk yield, but reduced (P < 0.05) milk fat and protein yield (on average, from 760 to 706 and 573 to 552 g/day, respectively). Relative to CPO, the LIN treatment reduced (P < 0.01) total saturated fatty acid content and enhanced (P < 0.001) 18:3n-3 in milk, whereas breed by diet interactions were significant for milk fat 16:0, total trans fatty acid and conjugated linoleic acid concentrations. Increases in 18:3n-3 intake derived from LIN in the diet were transferred into milk with a mean marginal transfer efficiency of 1.8%. Proportionate changes in milk fatty acid composition were greater in the Jersey, highlighting the importance of diet–genotype interactions on mammary lipogenesis. More extensive studies are required to determine the role of genotype on milk fat composition responses to oilseeds in the diet.

Production and metabolic effects of site of starch digestion in lactating dairy cattle

Milk solids yield in modern dairy cows has increased linearly over the last 50 years, stressing the need for maximal dietary energy intake to allow genetic potential for milk energy yield to be realized with minimal negative effects on health and reproduction. Feeding supplemental starch is a common approach for increasing the energy density of the ration and supplying carbon for meeting the substantial glucose requirement of the higher yielding cow. In this regard, it is a long held belief that feeding starch in forms that increase digestion in the small intestine and glucose absorption will benefit the cow in terms of energetic efficiency and production response, but data supporting this dogma are equivocal. This review will consider the impact of supplemental starch and site of starch digestion on metabolic and production responses of lactating dairy cows, including effects on feed intake, milk yield and composition, nutrient partitioning, the capacity of the small intestine for starch digestion, and nutrient absorption and metabolism by the splanchnic tissues (the portal-drained viscera and liver). Whilst there appears to be considerable capacity for starch digestion and glucose absorption in the lactating dairy cow, numerous strategic studies implementing postruminal starch or glucose infusions have observed increases in milk yield, but decreased milk fat concentration such that there is little effect on milk energy yield, even in early lactation. Measurements of energy balance confirm that the majority of the supplemental energy arising from postruminal starch digestion is used with high efficiency to support body adipose and protein retention, even in early lactation. These responses may be mediated by changes in insulin status, and be beneficial to the cow in terms of reproductive success and well-being. However, shifting starch digestion from the rumen impacts the nitrogen economy of the cow as well by shifting the microbial protein gained from starch digestion from potentially absorbable protein to endogenous faecal loss.

The effects of infection by Tetracapsuloides bryosalmonae (Myxozoa) and temperature on Fredericella sultana (Bryozoa)

The myxozoan, Tetracapsuloides bryosalmonae, exploits freshwater bryozoans as definitive hosts, occurring as cryptic stages in bryozoan colonies during covert infections and as spore-forming sacs during overt infections. Spores released from sacs are infective to salmonid fish, causing the devastating Proliferative Kidney Disease (PKD). We undertook laboratory studies using mesocosm systems running at 10, 14 and 20 degrees C to determine how infection by T bryosalmonae and water temperature influence fitness of one of its most important bryozoan hosts, Fredericella sultana, over a period of 4 weeks. The effects of infection were context-dependent and often undetectable. Covert infections appear to pose very low energetic costs. Thus, we found that growth of covertly infected F. sultana colonies was similar to that of uninfected colonies regardless of temperature, as was the propensity to produce dormant resting stages (statoblasts). Production of statoblasts, however, was associated with decreased growth. Overt infections imposed greater effects on correlates of host fitness by: (i) reducing growth rates at the two higher temperatures: (ii) increasing mortality rates at the highest temperature: (iii) inhibiting statoblast production. Our results indicate that parasitism should have a relatively small effect on host fitness in the field as the negative effects of infection were mainly expressed in environmentally extreme conditions (20 degrees C for 4 weeks). The generally low virulence of T. bryosalmonae is similar to that recently demonstrated for another myxozoan endoparasite of freshwater bryozoans. The unique opportunity for extensive vertical transmission in these colonial invertebrate hosts couples the reproductive interests of host and parasite and may well give rise to the low virulence that characterises these systems. Our study implies that climate change can be expected to exacerbate PKD outbreaks and increase the geographic range of PKD as a result of the combined responses of T. bryosalmonae and its bryozoan hosts to higher temperatures. Crown Copyright (C) 2009 Published by Elsevier Ltd. All rights reserved.

Fighting strategies in two species of fig wasp

Although theory exists concerning the types of strategies that should be used in contests over resources, empirical work explicitly testing its predictions is relatively rare. We investigated male fighting strategies in two nonpollinating. g wasp species associated with Ficus rubiginosa figs. In Sycoscapter sp. A, males did not assess each other before or during fights over mating opportunities. Instead,fights continued until the loser reached an energetic cost threshold that was positively correlated with its body size (fighting ability) and retreated. In Philotrypesis sp. B, pre fight assessment was indicated, with males attacking competitively inferior rivals to remove them from the competitor pool ( they then continued to do so until they reached a cost threshold that was again positively correlated with body size). Using data on species ecology, we discuss our findings with respect to theory on when different fighting strategies should evolve. We argue that the type of strategy used by a. g wasp species is determined by its relative benefits in terms of inclusive fitness. (c) 2008 The Association for the Study of Animal Behaviour. Published by Elsevier Ltd. All rights reserved.

Temporal and spatial variability of mobile fauna on a submarine cliff and boulder scree complex: a community in flux

The mobile component of a community inhabiting a submarine boulder scree/cliff was investigated at Lough Hyne, Ireland at dawn, midday, dusk and night over a 1-week period. Line transects (50 m) were placed in the infralittoral (6 m) and circumlittoral (18 m) zones and also the interface between these two zones (12 m). The dominant mobile fauna of this cliff consisted of echinoderms (6 species), crustaceans (10 species) and fish (23 species). A different component community was identified at each time/depth interval using Multi-Dimensional Scaling (MDS) even though both species diversity (Shannon-Wiener indices) and richness (number of species) remained constant. These changes in community composition provided indirect evidence for migration by these mobile organisms. However, little evidence was found for migration between different zones with the exception of the several wrasse species. These species were observed to spend the daytime foraging in the deeper zone, but returned to the upper zone at night presumably for protection from predators. For the majority of species, migration was considered to occur to cryptic habitats such as holes and crevices. The number of organisms declined during the night, although crustacean numbers peaked, while fish and echinoderms were most abundant during day, possibly due to predator-prey interactions. This submarine community is in a state of flux, whereby, community characteristics, including trophic and energetic relationships, varied over small temporal (daily) and spatial (m) scales.

Ultrafast excited state dynamics controlling photochemical isomerization of N-methyl-4-[trans-2-(4-pyridyl)ethenyl]pyridinium coordinated to a {Re-1(CO)(3)(2,2 '-bipyridine)} chromophore

Two Multifunctional photoactive complexes [Re(Cl)(CO)(3)-(MeDpe(+))(2)](2+) and [Re(MeDpe(+))(CO)(3)(bpy)](2+) (MeDpe(+) = N-methyl-4-[trans-2-(4-pyridyl)ethenyl]pyridinium, bpy = 2,2'-bipyridine) were synthesized. characterized. and their redox and photonic properties were investigated by cyclic voltammetry: ultraviolet-visible-infrared (UV/Vis/IR) spectroelectrochemistry, stationary UV/Vis and resonance Raman spectroscopy; photolysis; picosecond time-resolved absorption spectroscopy in the visible and infrared regions: and time-resolved resonance Raman spectroscopy. The first reduction step of either complex Occurs at about -1.1 V versus Fc/Fc(+) and is localized at MeDpe(+). Reduction alone does not induce a trans -> cis isomerization of MeDpe(+). [Re(Cl)(CO)(3)(MeDPe(+))(2)](2+) is photostable, while [Re(MeDpe(+))(CO)(3)(bpy)](2+) and free MeDpe(+) isomerize under near-UV irradiation. The lowest excited state of [Re(Cl)(CO)(3)(MeDPe(+))(2)](2+) has been identified as the Re(Cl)(CO)(3) -> MeDpe(+) (MLCT)-M-3 (MLCT = metal-to-ligand charge transfer), decaying directly to the ground state with lifetimes of approximate to 42 (73%) and approximate to 430ps (27%). Optical excitation of [Re(MeDpe(+))(CO)(3)(bpy)](2+) leads to population of Re(CO)(3) -> MeDpe(+) and Re(CO)(3) -> bpy (MLCT)-M-3 states, from which a MeDpe(+) localized intraligand 3 pi pi* excited state ((IL)-I-3) is populated with lifetimes of approximate to 0.6 and approximate to 10 ps, respectively. The 3IL state undergoes a approximate to 21 ps internal rotation, which eventually produces the cis isomer on a much longer timescale. The different excited-state behavior of the two complexes and the absence of thermodynamically favorable interligand electron transfer in excited [Re(MeDpe(+))(CO)(3)(bpy)](2+) reflect the fine energetic balance between excited states of different orbital origin, which can be tuned by subtle Structural variations. The complex [Re(MeDpe+)(CO)(3)(bpy)](2+) emerges as a prototypical, multifunctional species with complementary redox and photonic behavior.

Self-nucleation and crystallization kinetics of double crystalline poly(p-dioxanone)-b-poly(epsilon-caprolactone) diblock copolymers

The crystallization kinetics of each constituent of poly(p-dioxanone)-b-poly(epsilon-caprolactone) diblock copolymers (PPDX-b-PCL) has been determined in a wide composition range by differential scanning calorimetry and compared to that of the equivalent homopolymers. Spherulitic growth rates were also measured by polarized optical microscopy while atomic force microscopy was employed to reveal the morphology of one selected diblock copolymer. It was found that crystallization drives structure formation and both components form lamellae within mixed spherulitic superstructures. The overall isothermal crystallization kinetics of the PPDX block at high temperatures, where the PCL is molten, was determined by accelerating the kinetics through a previous self-nucleation procedure. The application of the Lauritzen and Ho. man theory to overall growth rate data yielded successful results for PPDX and the diblock copolymers. The theory was applied to isothermal overall crystallization of previously self-nucleated PPDX ( where growth should be the dominant factor if self-nucleation was effective) and the energetic parameters obtained were perfectly matched with those obtained from spherulitic growth rate data of neat PPDX. A quantitative estimate of the increase in the energy barrier for crystallization of the PPDX block, caused by the covalently bonded molten PCL as compared to homo-PPDX, was thus determined. This energy increase can dramatically reduce the crystallization rate of the PPDX block as compared to homo-PPDX. In the case of the PCL block, both the crystallization kinetics and the self-nucleation results indicate that the PPDX is able to nucleate the PCL within the copolymers and heterogeneous nucleation is always present regardless of composition. Finally, preliminary results on hydrolytic degradation showed that the presence of relatively small amounts of PCL within PPDX-bPCL copolymers substantially retards hydrolytic degradation of the material in comparison to homo-PPDX. This increased resistance to hydrolysis is a complex function of composition and its knowledge may allow future prediction of the lifetime of the material for biomedical applications.

Consequences of N,C,N '- and C,N,N '-coordination modes on electronic and photophysical properties of cyclometalated aryl ruthenium(II) complexes

The effects of isoelectronic replacement of a neutral nitrogen donor atom by an anionic carbon atom in terpyridine ruthenium(II) complexes on the electronic and photophysical properties of the resulting N,C,N'- and C,N,N'-cyclometalated aryl ruthenium(II) complexes were investigated. To this end, a series of complexes was prepared either with ligands containing exclusively nitrogen donor atoms, that is, [Ru(R-1-tpy)(R-2-tpy)](2+) (R-1, R-2 = H, CO2Et), or bearing either one N,C,N'- or C,N,N'-cyclometalated ligand and one tpy ligand, that is, [Ru(R-1-(NCN)-C-Lambda-N-Lambda)(R-2-tpy)](+) and [Ru(R-1-(CNN)-N-Lambda-N-Lambda)(R-2-tpy)](+), respectively. Single-crystal X-ray structure determinations showed that cyclometalation does not significantly alter the overall geometry of the complexes but does change the bond lengths around the ruthenium(II) center, especially the nitrogen-to-ruthenium bond length trans to the carbanion. Substitution of either of the ligands with electron-withdrawing ester functionalities fine-tuned the electronic properties and resulted in the presence of an IR probe. Using trends obtained from redox potentials, emission energies, IR spectroelectrochemical responses, and the character of the lowest unoccupied molecular orbitals from DFT studies, it is shown that the first reduction process and luminescence are associated with the ester-substituted C,N,N'-cyclometalated ligand in [Ru(EtO2C-(CNN)-N-Lambda-N-Lambda)(tpy)](+). Cyclometalation in an N,C,N'-bonding motif changed the energetic order of the ruthenium d(zx), d(yz), and d(xy) orbitals. The red-shifted absorption in the N,C,N'-cyclometalated complexes is assigned to MLCT transitions to the tpy ligand. The red shift observed upon introduction of the ester moiety is associated with an increase in intensity of low-energy transitions, rather than a red shift of the main transition. Cyclometalation in the C,N,N'-binding motif also red-shifts the absorption, but the corresponding transition is associated with both ligand types. Luminescence of the cyclometalated complexes is relatively independent of the mode of cyclometalation, obeying the energy gap law within each individual series.

Polymer-mediated disruption of drug crystallinity

Ibuprofen (IB), a BCS Class II compound, is a highly crystalline substance with poor solubility properties. Here we report on the disruption of this crystalline structure upon intimate contact with the polymeric carrier cross-linked polyvinylpyrrolidone (PVP-CL) facilitated by low energy simple mixing. Whilst strong molecular interactions between APIs and carriers within delivery systems would be expected on melting or through solvent depositions, this is not the case with less energetic mixing. Simple mixing of the two compounds resulted in a significant decrease in the differential scanning calorimetry (DSC) melting enthalpy for IB, indicating that approximately 30% of the crystalline content was disordered. This structural change was confirmed by broadening and intensity diminution of characteristic IB X-ray powder diffractometry (PXRD) peaks. Unexpectedly, the crystalline content of the drug continued to decrease upon storage under ambient conditions. The molecular environment of the mixture was further investigated using Fourier transform infrared (FT-IR) and Fourier transform Raman (FT-Raman) spectroscopy. These data suggest that the primary interaction between these components of the physical mix is hydrogen bonding, with a secondary mechanism involving electrostatic/hydrophobic interactions through the IB benzene ring. Such interactions and subsequent loss of crystallinity could confer a dissolution rate advantage for IB. (C) 2006 Elsevier B.V. All rights reserved.

Methylxanthines are the psycho-pharmacologically active constituents of chocolate

Rationale: Liking, cravings and addiction for chocolate ("chocoholism") are often explained through the presence of pharmacologically active compounds. However, mere "presence" does not guarantee psycho-activity. Objectives: Two double-blind, placebo-controlled studies measured the effects on cognitive performance and mood of the amounts of cocoa powder and methylxanthines found in a 50 g bar of dark chocolate. Methods: In study 1, participants (n=20) completed a test battery once before and twice after treatment administration. Treatments included 11.6 g cocoa powder and a caffeine and theobromine combination (19 and 250 mg, respectively). Study 2 (n=22) comprised three post-treatment test batteries and investigated the effects of "milk" and "dark" chocolate levels of these methylxanthines. The test battery consisted of a long duration simple reaction time task, a rapid visual information processing task, and a mood questionnaire. Results: Identical improvements on the mood construct "energetic arousal" and cognitive function were found for cocoa powder and the caffeine+theobromine combination versus placebo. In chocolate, both "milk chocolate" and "dark chocolate" methylxanthine doses improved cognitive function compared with "white chocolate". The effects of white chocolate did not differ significantly from those of water. Conclusion: A normal portion of chocolate exhibits psychopharmacological activity. The identical profile of effects exerted by cocoa powder and its methylxanthine constituents shows this activity to be confined to the combination of caffeine and theobromine. Methylxanthines may contribute to the popularity of chocolate; however, other attributes are probably much more important in determining chocolate's special appeal and in explaining related self-reports of chocolate cravings and "chocoholism".