Stack Exchange Network

Stack Exchange network consists of 175 Q&A communities including Stack Overflow, the largest, most trusted online community for developers to learn, share their knowledge, and build their careers.

Visit Stack Exchange

Questions tagged [chemistry]

For questions related to chemical concepts such as solving equations related to reaction kinetics, graphically visualizing molecules, and using the chemistry-related curated data sets.

16
votes
4answers
622 views

How to apply different equations to different parts of a geometry in PDE?

I want to solve two coupled partial differential equations on two dimension. There are two variables v and m. The geometry is a ...
3
votes
1answer
83 views

Producing MoleculePlots of pigments

Of course no chemical database is complete, but I would like to create as large a list of MoleculePlots of pigments (organic or inorganic) as is available through ...
0
votes
1answer
55 views

How do I obtain the characteristics and properties of an alloy? [closed]

I would like to obtain the characteristics and properties of the Alloy Chromel. How do I download that information from the database of Wolfram Alpha and then print it in a Mathematica Notebook? ...
6
votes
1answer
164 views

Maxwell construction problem with FindRoot

A Maxwell construction is just an integral of area and then solving for a root to find vapor pressure. This must be exactly the same value as where the Gibbs energy plot crosses at a point. I found ...
11
votes
3answers
439 views

Minimum energy path of a potential energy surface

I have calculated the energies of a molecule by varying two parameters (Data points are here data). I want to find the minimum energy path from the point (180.0, 179.99, 21.132) to (124.5, 124.49, 0) ...
4
votes
2answers
120 views

Is it possible to use chemical data in polymers?

I would like to plot the ABS polymer. It is made of Styrene, Acrylonitrile and 1,3 Butadine. I can plot those components separately But the command ...
9
votes
3answers
372 views

Chemical data of Lactic acid [duplicate]

I'm using the next Mathematica command in order to plot molecules: ChemicalData["name", "MoleculePlot"] But I'm not sure if the Mathematica database of molecules ...
1
vote
0answers
50 views

ProteinData Error: “Cannot import data as PDB format” [closed]

I'm trying to use the Mathematica examples for plotting proteins: But when I type the command I get the next error: How can I import data as PDB format?
13
votes
1answer
673 views

visualizing a normal mode vibration of a molecule in mathematica

How can I visualize the normal mode vibrations of a molecule in Mathematica? Something like the following GIF: (from here) I have the XYZ file of the benzene molecule and the normal mode ...
4
votes
2answers
246 views

Periodic table with color code

Is it possible to do this table (from here) in Mathematica? In this table, the color code indicates the element’s dipolar character $(~ m\mu_m^2)$ in atomic units.
2
votes
1answer
58 views

Visualizing a imported XYZ-file and graphic object at same time

I want visualize a imported xyz-file and a graphic object at same time. g = Import["H2O.xyz", "XYZ"]; Show[g, Graphics3D[{Opacity[0.5], Sphere[]}]] But this ...
2
votes
1answer
56 views

Convert a Java object to interface via J/Link

I'm trying to use a Java library called CDK with Mathematica's J/Link. ...
10
votes
3answers
417 views

DFT for Helium atom

I'm trying to solve this problem I'm not a Mathematica expert, but my program works perfect when vhartree={0,0,..,0} This is the program I wrote: ...
7
votes
1answer
140 views

Suggestions to increase the speed of a fugacity calculation

I'm writing code to simulate a chemical reaction. This requires the calculation of the fugacity of a mixture of gases. The underlying model is the SRK equation of state. First the physical ...
1
vote
0answers
73 views

How can i set up an If or Do loop to follow this calculation procedure?

The following picture depicts the calculation procedure I wish to set up. The inputs are $T_r$, the chemical composition of the feedstock and constraints defined by the equivalence ratio. This ...
5
votes
1answer
91 views

WolframAlpha function does not return structure diagram for chemicals

Neither one of WolframAlpha["Caffeine"] WolframAlpha["Caffeine","FullOutput"] returns a structure diagram for the molecule. It's consistent since ...
0
votes
0answers
63 views

Encountering error messages with both DSolve and NDSolve

Reposting my question in a more legible and understandable way. The earlier one was flouting a few too many posting guidelines. I am trying to solve the heat conduction partial differential equation, ...
8
votes
2answers
404 views

ChemistryData — where is specific gravity?

I cannot seem to find the basic chemical property specific gravity with ChemicalData or other comparable functions.
1
vote
2answers
152 views

Assign a unit to a variable in a function

I have the following input: PVW[V_, T_] := (R*T/(V - b)) - (a/(V^2)) Where I have already defined R as the molecular gas constant with its corresponding units ...
1
vote
0answers
57 views

Dealing with Vector Outputs from ItoProcess, RandomProcess (Stochastic Differential Equations)

I'm modeling stochastic chemical kinetics and the ItoProcess[] function has served me well. I am trying to write an efficient code to analyze many (Paths) trajectories of several different reagents (...
6
votes
1answer
512 views

Radial distribution function (pair correlation function) for particles in rectangular image

Some time ago I asked the following question on how to calculate the radial distribution function. I have here an add-on question, but before asking it I try to summarize and visualize my problem. ...
4
votes
1answer
339 views

Map volumetric data onto an isosurface

this question originates from an attempt to visualize molecular electrostatic potentials, but hopefully should be of general interest. Suppose that you have two sets of data, S1 and S2. Each set ...
5
votes
2answers
124 views

Delete duplicates of isomerism

I written some code to generate all partens of $\ \mathrm{C_2H_5O}$, but theres some isomerisms need to be deleted. Here's the raw output after DeleteDuplicates. <...
7
votes
1answer
126 views

Identify chemical expressions

Background I written a chemical equation balancing program. ...
1
vote
2answers
246 views

Visualizing .xyz file in Mathematica

If i import a .xyz file in Mathematica that will show me the moelcule in the default way. But i want to visualise the .xyz file in the following manner.
5
votes
1answer
158 views

Michaelis–Menten kinetics using Tau-Leaping

I am trying to implement Michaelis–Menten kinetics using Tau-Leaping. It works fine as expected. Only issue here is when substrate reaches zero the algorithm stop since PoissonDistribution[0] gives ...
8
votes
1answer
230 views

Constructing Icosahedral Fullerene Graphs

I am interested in constructing graphs of icosahedral fullerenes, following the construction explained in the first few pages of this article (alternate link). Mathematica has built-in information ...
2
votes
2answers
349 views

Solve for several constants using a set of initial conditions

I have a set of differential equations that I am trying to solve for chemical reaction rate constants. So far I can get DSolve to solve the ODEs with initial conditions which provides the chemical ...
0
votes
2answers
311 views

Michaelis–Menten kinetics using Gillespie Algorithm

I am trying to implement Michaelis–Menten kinetics using Gillespie Algorithm. $$\text{Reaction 1:}\quad E+S\rightarrow C$$ $$\text{Reaction 2:}\quad C\rightarrow E+S$$ $$\text{Reaction 1:}\quad C\...
0
votes
2answers
106 views

Solving Catalyst/Enzyme Diffusion Equation

Evening all, I am attempting to solve a second order differential equation to get an equation of the substrate concentration as a function of r (i.e. S(r)). I have developed this equation, $\frac{d^...
8
votes
3answers
768 views

Why are Dalton and Grams / Moles not equivalent units?

The following code throws an error: UnitConvert[Quantity["Daltons"], "Grams"/"Moles"] Why?
2
votes
1answer
118 views

How to use FindFit to do kinetic simulation?

I'd like to simulate the chemical kinetic model by fitting the experimental data. Here are my codes for simulation. ...
5
votes
0answers
75 views

Helping Mandy become a better speller [duplicate]

tl;dr Can you improve upon this code, which checks to see if a word can be spelled using the element symbols from the periodic table? I recently created a periodic table that lights up. The design ...
2
votes
2answers
106 views

Question about getting the table behind ElementData

I am interested in retrieving the entire table (preferably in text or csv format) that is behind the ElementData function for ALL the elements. The table (or perhaps multiple tables) record(s) many ...
1
vote
2answers
97 views

How can I add new columns to a Table after each evaluation?

I'm interested in simulating chemical reactions with perturbations. I can simulate a given reaction using NDSolve ("rxn"} with a given added noise component ("noise1"). Due to the noise, each ...
9
votes
2answers
310 views

Plot polygonal surface

I have a list of 120 elements of the form $\{q,x,y,z\}$. In fact, each element of the represents the nuclear charge and 3 coordinates of an atom of some large molecule. ...
15
votes
3answers
705 views

Proteins, ProteinData, WolframAlpha

Apologies in advance, this is the first I've messed with ProteinData and access through WolframAlpha, so clearly I may not know ...
1
vote
2answers
2k views

How to plot Compressibility factor Z vs Pressure P using ParametricPlot?

I am trying to get my plot to look like this: but I cannot seem to get it right when I try to plot each curve as a parametric equation. Here is what I have been trying: ...
1
vote
0answers
121 views

DNA Modeling In Three Dimensions? [closed]

I can find any number of base pair nucleotide sequences for any number of genes in the fasta format. (Check out the COSMIC Tumor Database, for example). I'm having trouble finding base pair nucleotide ...
2
votes
1answer
319 views

Connectivity in a molecule and permutations

What am I trying to do is to make a script that calculates all chemically reasonable torsion angles in any molecule. For one torsional angle I need 4 CONNECTED atoms, 1-2-3-4, 1-2-3 defining one ...
6
votes
1answer
381 views

Chemical Kinetics: fitting multiple sets of time series data using multiple differential reaction equations

Let me start by saying that this is my first post, so if there are suggestions on proper formatting that is annoying you, please let fly and I will do my best to modify. I was doing well solving the ...
3
votes
1answer
116 views

How do I get the MeanPredictionBands from the fitting of the numerical solution of coupled differential equations to data?

This is my first post, and I will preface my question by stating that I am relatively new to Mathematica. In any case, I would like to fit the solution of chemical kinetics equations to data and get ...
2
votes
0answers
57 views

Another question on Mma data sources/syntax

Let's say I want to have the vapor pressure of mercury, as a function of temperature. I know that this information is available, for any specific temperature, from Wolfram Alpha, and I can also get ...
2
votes
1answer
158 views

How can I plot all solutions for a system of equations across a given domain for a selected variable?

I have a coupled chemical equilibria system consisting of 12 equations in 12 unknowns. I have managed use the Solve function to write a command that will find a particular solution at any fixed value ...
1
vote
0answers
74 views

Representing vertices as spheres in Graphics 3D [closed]

I'm representing a portion of molecule, P3HT, in Mathematica as a pentagon that I build using Polygon and Graphics3D. I would like to represent the vertices of the pentagon as spheres(corresponding ...
23
votes
1answer
2k views

Gillespie Stochastic Simulation Algorithm

The Gillespie SSA is a Monte Carlo stochastic simulation algorithm to find the trajectory of a dynamic system described by a reaction (or interaction) network, e.g. chemical reactions or ecological ...
-1
votes
1answer
146 views

Problem with a chemical optimization problem [closed]

First of all, I want to say that I'm quite new to Mathematica. I searched a lot on this site to learn how to solve my problem, which this is the problem of my life! I want to minimize a function. It ...
7
votes
3answers
453 views

How can I draw ternary plot by these data?

I have these data and want to draw similar fifures in mathematica. any one can help me? ...
2
votes
3answers
126 views

Parsing/validating CAS registry number identifiers

Mathematica includes a lot of curated information on molecular structure and properties available through ChemicalData and ...
6
votes
2answers
220 views

Atom weight in SI units

In chemistry is very common to use the units grams/mol for the atomic weight, so the mass of one atom (of Tin, for example) is calculated as : massOfOneAtom = (0.11871 (Kg/mol))/(AvogadroConstant (1/...