Recent posts

#71

The MELTS family of algorithms, and scientific discussion / Re: setting oxygen fugacity fo...

Last post by Paula - May 23, 2024, 03:38:36 PM

Hi Melina,

I don't think there is an fO2 example in the MELTS-pMELTS directory. There is functionality like that in the newer Equilibrate package but I haven't looked at any of it lately. Your best bet might be to ask on the ENKI Slack channel (enki-portal.slack.com). It's the free version of Slack so you can't see old posts, unfortunately, but queries usually get answered by community members quickly.

Sanidine is an end member in the feldspar solid solution in (Rhyolite-)MELTS. So far as I know, you can't exclude it in Thermoengine without excluding the feldspar phase altogether, although you could might be able to  get around this by removing K2O from your bulk composition.

Paula

#72

The MELTS family of algorithms, and scientific discussion / setting oxygen fugacity for mo...

Last post by Melina - May 21, 2024, 05:49:04 AM

Hello everyone,

I have been trying to model the LLD for the volcanic system that I'm studying. I've been using the "MELTS v1.0.2 fractionation" Jupyter notebook on the ENKI portal sever for this and ran into a problem. I can't figure out how to change/set the oxygen fugacity? Can anyone give me a hint?

Also, I have excluded the phase sanidine from my calculations, but it is still included in my outputs? It doesn't happen with any other phase that I have excluded. Is there another way to exclude phases except with melts.set_phase_inclusion_status()?

#73

The MELTS family of algorithms, and scientific discussion / Re: How to explain the contino...

Last post by Paula - May 01, 2024, 05:49:59 PM

The main problem is that there aren't many solid phases in the MELTS models that take MnO or P2O5.

The only Mn-bearing solid phases are olivine and rhombohedral oxide (ilmenite). If ilmenite is not joining the assemblage early enough, you could try playing around with the fO2.

The only P-bearing solid phases are (hydroxy)apatite and whitlockite. MELTS underestimates the stability of apatite, mainly because of the lack of fluorine and chlorine in the system. There is a F, Cl-bearing version of MELTS under development. At the moment the most recent citation is this:

[CITATION] MELTS+ DEW: Modeling major element+ Cl+ F+ S phase equilibria, redox reactions and elemental partitioning in magmatic-hydrothermal systems
MS Ghiorso, S Matthews... - Goldschmidt 2023 ..., 2023 - conf.goldschmidt.info

Look out for more developments in the next few months.

Paula

#74

The MELTS family of algorithms, and scientific discussion / How to explain the continous i...

Last post by Wei-cheng Jiang - April 21, 2024, 04:53:14 PM

Hi everyone,

I tried to use Rhyolite-MELTS 1.2.0 to simulate the LLD evolution in my study area. It works quite well in most of the major oxides like MgO, FeO, Al2O3, CaO, Alkalis, and SiO2. Our phase-equilibria experiments also show a similar LLD pathway under similar conditions.

However, the main issue is that the concentration of MnO and P2O5 keep increasing throughout the magma evolution from my parental composition to rhyolite composition. This problem was also pointed out by some reviewers for my submitted paper. The phase assemblage does not contain any manganese- or phosphorus-rich minerals (the starting composition is a high-Mg andesite with plag + cpx + opx + ol + spinel + ilmenite). I am not sure if this is the reason.

Thank you!

#75

easyMelts / Re: Restarting from an earlier...

Last post by Paula - April 16, 2024, 02:31:33 PM

Hi Kevin,

I haven't managed to reproduce convergence failure, as opposed to just hanging, but there are two ways I would try to tackle it:

1. Write out a .melts file of the current state. Edit the file, then restart easyMelts and read in the file. This may be your only option for fractional melting.

2. Start from the top with "Set initial T & P", Save, "Set initial composition", Save. I think the Equilibrate button will reappear.

Note that, unlike the rhyolite-MELTS GUI, in easyMelts T, P, and composition are not updated in the input boxes as the calculation proceeds. So you may need to copy them from the last good step (if it's fractional crystallization) and/or the initial conditions (if it's a batch calculation). Drag the slider to the right to get the last good step.

Paul @asimow uses eayMelts more than me so may have a better workaround.

Paula

#76

easyMelts / Restarting from an earlier ste...

Last post by KevinFarries - April 13, 2024, 09:41:51 PM

Hi Paula,

I was previously using the rhyolite-Melts GUI and have switched across to easyMelts.  In rhyolite-Melts if the algorithm fails to converge, it is possible to tweak some of the parameters then continue with the analysis.  Is there a way to do this with easyMelts?  If the there is a convergence failure, the "equilibrate" button disappears and I can't see a way to continue the calculation.

Thanks!

Kevin

#77

alphaMELTS for MATLAB/Python / Re: How to set elements to be ...

Last post by Paula - April 02, 2024, 08:58:54 AM

Hi Isiah,

No, sorry, as we have decided to take a different route and have the (passive) trace elements handled in Python, so we can integrate with existing packages for trace element partitioning models etc. The water as a trace species (i.e. pHMELTS) is close to implemented. That will be integrated into the main alphaMELTS program more thoroughly, rather than as a tack on trace species engine. So it will feel similar to the interface to the fO2 buffering options.

Best wishes,
Paula

#78

The MELTS family of algorithms, and scientific discussion / Re: How to simulate isentropic...

Last post by Paula - April 02, 2024, 08:47:36 AM

Thanks for the video! I'm still not quite sure why it's working but I'm less surprised that it might work on Ubuntu 20.04 than 22.04, especially if that was upgraded from an earlier Ubuntu. Anyway I'm glad you have easyMelts working now - it's definitely the more convenient option these days.

Best wishes,
Paula

#79

The MELTS family of algorithms, and scientific discussion / Re: How to simulate isentropic...

Last post by Gary - April 01, 2024, 04:02:59 AM

Hi, Paula,

Thank you very much for your detailed answer! Followed your steps, my easyMELTS runs well. I will use easyMELTS next to simulate mantle partial melting because it is more convenient to operate. If I meet other difficulties, I will ask for your help.

I know you have some question about why the rhyolite-MELTS/pMELTS can run on my Ubuntu. I'm sorry I can't answer you well because I'm not familiar with these operations. I recorded my operating steps, maybe it will be helpful: https://youtu.be/_NW-fTJHNHA. If you have any questions about this, just tell me.

All the best
Gary

#80

alphaMELTS for MATLAB/Python / Re: How to set elements to be ...

Last post by IsaiahSpring - March 26, 2024, 11:34:01 AM

Hi Paula,

I was wondering if this Trace element capability was up and running, and if there were any materials or examples about its operation.

Thanks,
Isaiah