Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Exact solver translate to Julia #2

Open
PaulWAyers opened this issue Nov 14, 2022 · 9 comments
Open

Exact solver translate to Julia #2

PaulWAyers opened this issue Nov 14, 2022 · 9 comments

Comments

@PaulWAyers
Copy link
Member

https://github.com/QuantumElephant/NICE/blob/master/nice/exact.py

@PaulWAyers
Copy link
Member Author

@PaulWAyers
Copy link
Member Author

PaulWAyers commented Jan 23, 2023

The constraint is that there can never be a negative amount of a reactant. So for reaction R, $\xi_R$ is greater than the -1 times the smallest positive stoichiometric coefficient and larger than the -1 times the largest negative stoichiometric coefficient.

@PaulWAyers
Copy link
Member Author

Initial guess:

  • specified by user.
  • from net-event KMC run.
  • randomly selected in the allowed interval of $\xi_R$ for each reaction.
  • selected to be zero.
  • more sophisticated choices are possible but perhaps not worthwhile. E.g., you could choose $xi_R$ to be greater than zero when $K_R > 1$ and less than zero otherwise.

@PaulWAyers
Copy link
Member Author

I'd probably give two options:

  1. User specified (can include user passing result from a previous simulation).
  2. Zero (if user doesn't specify).

@PaulWAyers
Copy link
Member Author

PaulWAyers commented Jan 23, 2023

It may be worth using the optim.jl package, as it doesn't require building the sum-of-squares objective function.

https://github.com/JuliaNLSolvers/Optim.jl

@PaulWAyers
Copy link
Member Author

It looks like there is a high-level interface with the Neqproblem keyword
https://github.com/JuliaNLSolvers/NLSolvers.jl

@NamanKothari5
Copy link

https://github.com/QuantumElephant/NICE/blob/master/nice/exact.py

@PaulWAyers as of now the URL is not working I wanted to attempt on the exact solver implementation can you point me to some relevant materials to get started with this.

Also, I wish to connect with you mentors and other members of the community, so, it would be beneficial if you can provide Slack, Element, Discord or any other communication channel link.

@PaulWAyers
Copy link
Member Author

Some notes on simultaneous equilibria are here. The first is the general framework, then a set of notes that has a lot of examples.
General Notes
General + Specific Examples

The supplementary material from J. Mol. Modelling v23, 155, (2017) is also relevant.

@magnamad
Copy link

Hi,

I am Madhavan, currently pursuing my masters from IISc Bangalore, India in the department of Chemical Engineering. @PaulWAyers First of all, thank you for the resources shared regarding the setting up of nonlinear equations. Additionally, can you also provide some resources pointing to how the KMC simulations needs to be implemented.
It would be great if you could share the link for some contact channel so that I can connect with the members of your community.

Thanks again!
Madhavan

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

3 participants