# Tag Info

6

You're trying to have your cake and it it too. This does not work. As a general rule, for problems with features on different length scales, you need meshes that are fine in at least some parts of the mesh. This results in many cells, and this results in long computations, small time steps, and many linear iterations. All of these implications are rather ...

5

With conforming triangular meshes, it will be difficult to make an isotropic mesh which adapts to multiple dramatically different length scales in such a short space without introducing extraneous triangles, some of which may have very large/small angles. I'm not very familiar with them so take this with a grain of salt, but you may have better luck using ...

4

So Is it neccesary to use Poisson - Boltzmann equation if I only need to build electrostatic potential from a PQR file? No. You can use Poisson. Since you know the positions of each point charge, you know the charge distribution $\rho$, which is a sum of delta functions. You can thus solve numerically the Poisson's equation that links the charge ...

3

Your model seems to have a single stable critical point that attracts all orbits. You're plotting $(S(t), I(t), R(t))$ as functions of time for specific initial conditions, but actually the easier way to see this directly is to plot the vector field $(\dot S, \dot I)$ in the $S$-$I$ plane like so: Immediately you can see that for this particular choice of ...

3

The answer as of 2017, to the best of my knowledge, is "not yet." Your question is treated in detail here and in the references therein.

3

The Van der Walls radius - last column in the output - is calculated from the force field. Probably pdb2pqr and editconf uses different force fields, hence different radius. I don't use pdb2pqr, but it seems (1) AMBER99 is the default force field, though CHARMM, PARSE and TYL06 are supported. Gromacs' editconf reads force field from topology file generated ...

3

I'm not particularly familiar with biofilm literature. But in most computational literature, the border of the entire domain of a problem are usually referred to as the boundary. Outside of a boundary, there are no nodes, elements, or anything else under consideration. The entire domain may also be subdivided into smaller regions. Some of these regions ...

2

The AIC function need an 'lm' or 'glm' object (linear models). See functions lm and glm So just do : AIC(lm(logCPK~dataPOW))

2

There is something called the ODD (Overview, Design, and Details) protocol, proposed by Volker Grimm and others for describing an agent based model. It consists of a list of elements that are needed for understanding the functioning of an ABM and aims at making descriptions of such models more standardised. The checklist of what has to be described consists ...

2

Is there a reason you can't just call that function and then plot your scatter plot on the figure it creates? Otherwise you could either open up the existing function that plots the Ramachandran regions and copy/modify the relevant parts to do what you want or you could copy the data out of the plotted objects (using the get function, see this link) and re-...

2

Molecular dynamics, in the computation of interactions inversely proportional to distance. Other N-body codes may share this property, but it's often an explicit goal to minimize or eliminate explicit division for exactly the reason you've given.

1

You have AlphaNa = 0.1 * f(0.1*(VoltageDelta-10.0)) where $$f(x)=\frac{x}{1-\exp(-x)}$$ The limit for $x\to 0$ can actually be interpreted as difference quotient, you might also appeal to l'Hopital, to find that $f(x)\to 1$. So if you want to play it safe, you can employ a 3-tier or 4-tier implementation of this function def f(x): if abs(x)<1e-12: ...

1

Resolving small features in FEM will always be costly, there is no getting away from that fact. Your problem seems to be framed in terms of computational burden. In my own case, I was looking at electric field problems in anatomical structures, so had a similar set of problems to your own. The question is usually how detailed a mesh is "good enough" for the ...

1

So I managed to create a spiral finally. The method I used to initiate spiral is by initializing the mesh with a set of initial values expect a small square patch which is initialized with $180°$ phase shift with others and some having $90°$ phase shift. The link to spiral video I created is here.

1

Pretty much any biology can be made computational. Now as it is called "computational" I hope you can do some programming, for that's what it is going to be about. For gene stuff I think databases and algorithms are important, for they deal with big data. This is an area I know little about, but it seems what other people have been suggesting are mainly ...

1

I am surprised no one told about LAMMPS. It is pretty nice, and do have GPU capabilities.

1

Depends on whether you want free or are happy to pay. Have used CHARMM in the past. Here is a decent list of options (https://en.wikipedia.org/wiki/List_of_software_for_molecular_mechanics_modeling)

Only top voted, non community-wiki answers of a minimum length are eligible