<div dir="ltr">Thanks for clarifying its purpose! I guess I do still have one question, though: are ccp4 files computed via phenix.maps (from e.g., a .pdb and .mtz) unsuitable for real-space refinement (e.g., some kind of structural optimization with a scoring function that includes an electron density term) because such a ccp4 file is itself influenced by the starting model?</div><div class="gmail_extra"><br><div class="gmail_quote">On Wed, May 25, 2016 at 1:57 PM, Dale Tronrud <span dir="ltr">&lt;<a href="mailto:detBB@daletronrud.com" target="_blank">detBB@daletronrud.com</a>&gt;</span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><br>
   An mtz file is a generic container for information that is<br>
represented as Fourier coefficients - amplitudes and/or phases.  The mtz<br>
you get from data reduction contains your observed diffraction pattern.<br>
There are quite a variety of things that might be in there, e.g.<br>
intensities, amplitudes, indications of anomalous scattering.<br>
<br>
   The &quot;map_coeffs&quot; mtz probably echos all of those data and includes<br>
the amplitudes/phases which can be used to calculate several maps.<br>
These coefficients are derived from both your observations and the<br>
coordinates of your model that matches the &quot;map_coeffs&quot; file.  The<br>
program could simply write a map file directly, but the Fourier<br>
representation takes less disk space and is quicker to read and write<br>
than a map.<br>
<br>
   The &quot;map_coeffs&quot; file is a &quot;throw away&quot; file.  You use it to look at<br>
your map, but when you continue refinement you should always use the<br>
original mtz.  You don&#39;t want to refine your model against an mtz that<br>
contains data calculated from your model.  While it might be okay, it is<br>
safer to stick with the original.<br>
<br>
Dale Tronrud<br>
<span class=""><br>
On 5/25/2016 10:38 AM, Andy Watkins wrote:<br>
&gt; Thanks--that certainly makes sense. In that case, I suppose I have a<br>
&gt; different question: what distinguishes the .mtz file suffixed<br>
&gt; &quot;map_coeffs,&quot; which is output by phenix.maps, from the input .mtz?<br>
&gt;<br>
&gt; On Wed, May 25, 2016 at 1:28 PM, Christian Roth<br>
</span><div><div class="h5">&gt; &lt;<a href="mailto:christianroth034@gmail.com">christianroth034@gmail.com</a> &lt;mailto:<a href="mailto:christianroth034@gmail.com">christianroth034@gmail.com</a>&gt;&gt; wrote:<br>
&gt;<br>
&gt;     You might look for the Real space correlation coefficient to see if<br>
&gt;     your model fit that particular map better than before, but thats<br>
&gt;     probably it.<br>
&gt;<br>
&gt;     There is no subset experimental data just influenced by one chain.<br>
&gt;     The whole complex contributes to each measurement and vice versa<br>
&gt;     your chain A model contribute to each measurement you made to<br>
&gt;     generate the maps. Sorry, but that&#39;s it<br>
&gt;<br>
&gt;     Christian<br>
&gt;<br>
&gt;     Am 25.05.2016 um 17:31 schrieb Andy Watkins:<br>
&gt;&gt;     Perhaps necessary preface: I am largely learning as I go, and<br>
&gt;&gt;     large parts of Phenix (and the details of refinement data file<br>
&gt;&gt;     formats) are still black boxes for me.<br>
&gt;&gt;<br>
&gt;&gt;     Suppose I have a model of one chain of a complex and an .mtz<br>
&gt;&gt;     describing the entire complex. I have performed some refinement<br>
&gt;&gt;     procedure altering the model of that single chain.<br>
&gt;&gt;<br>
&gt;&gt;     My aim is to determine how this refinement procedure has affected<br>
&gt;&gt;     the fit to the experimental data--or, really, the subset of the<br>
&gt;&gt;     experimental data expected to be relevant for my model.<br>
&gt;&gt;<br>
&gt;&gt;     For obvious reasons,<br>
&gt;&gt;<br>
&gt;&gt;     phenix.model_vs_data chain_A.pdb whole_complex.mtz<br>
&gt;&gt;     phenix.model_vs_data chain_A_refined.pdb whole_complex.mtz<br>
&gt;&gt;<br>
&gt;&gt;     are undesirable: they evaluate the model of a single chain against<br>
&gt;&gt;     the whole complex&#39;s data, and there&#39;s obviously a lot of<br>
</div></div>&gt;&gt;     unsatisfied electron density. So, while this command /works/, it<br>
<span class="">&gt;&gt;     reports an unrepresentative fit. (In theory, one solution to my<br>
&gt;&gt;     problem would be to re-combine each chain&#39;s refined version, then<br>
&gt;&gt;     run model_vs_data on the recombined, refined complex. I&#39;m<br>
&gt;&gt;     interested nonetheless in how it would work on the chains separately.<br>
&gt;&gt;<br>
&gt;&gt;     Now, during the refinement procedure, we of course generate .ccp4<br>
&gt;&gt;     density maps for the individual chain models:<br>
&gt;&gt;<br>
&gt;&gt;     phenix.maps chain_A.pdb whole_complex.mtz<br>
&gt;&gt;<br>
&gt;&gt;     which produces the ccp4 file as well as, crucially,<br>
&gt;&gt;     chain_A_map_coeffs.mtz. Attempting to employ that resulting .mtz<br>
&gt;&gt;     file, i.e.<br>
&gt;&gt;<br>
&gt;&gt;     phenix.model_vs_data chain_A.pdb whole_complex.mtz<br>
&gt;&gt;     phenix.model_vs_data chain_A_refined.pdb whole_complex.mtz<br>
&gt;&gt;<br>
&gt;&gt;     leads to an evocative error: Sorry_No_array_of_the_required_type:<br>
&gt;&gt;     No reflection arrays available.<br>
&gt;&gt;<br>
&gt;&gt;     My highest suspicion is that I need to alter maps.params in a<br>
&gt;&gt;     particular way so that reflection arrays are also output to<br>
&gt;&gt;     model_map_coeffs.mtz. I could also imagine that I need to be using<br>
&gt;&gt;     another program entirely!<br>
&gt;&gt;<br>
&gt;&gt;     Thanks in advance for whatever help you can provide;<br>
&gt;&gt;     unfortunately, I can&#39;t provide any input files.<br>
&gt;&gt;<br>
&gt;&gt;<br>
&gt;&gt;     _______________________________________________<br>
&gt;&gt;     phenixbb mailing list<br>
</span>&gt;&gt;     <a href="mailto:phenixbb@phenix-online.org">phenixbb@phenix-online.org</a> &lt;mailto:<a href="mailto:phenixbb@phenix-online.org">phenixbb@phenix-online.org</a>&gt;<br>
&gt;&gt;     <a href="http://phenix-online.org/mailman/listinfo/phenixbb" rel="noreferrer" target="_blank">http://phenix-online.org/mailman/listinfo/phenixbb</a><br>
&gt;&gt;     Unsubscribe: <a href="mailto:phenixbb-leave@phenix-online.org">phenixbb-leave@phenix-online.org</a><br>
&gt;&gt;     &lt;mailto:<a href="mailto:phenixbb-leave@phenix-online.org">phenixbb-leave@phenix-online.org</a>&gt;<br>
&gt;<br>
&gt;<br>
&gt;     _______________________________________________<br>
&gt;     phenixbb mailing list<br>
&gt;     <a href="mailto:phenixbb@phenix-online.org">phenixbb@phenix-online.org</a> &lt;mailto:<a href="mailto:phenixbb@phenix-online.org">phenixbb@phenix-online.org</a>&gt;<br>
&gt;     <a href="http://phenix-online.org/mailman/listinfo/phenixbb" rel="noreferrer" target="_blank">http://phenix-online.org/mailman/listinfo/phenixbb</a><br>
&gt;     Unsubscribe: <a href="mailto:phenixbb-leave@phenix-online.org">phenixbb-leave@phenix-online.org</a><br>
&gt;     &lt;mailto:<a href="mailto:phenixbb-leave@phenix-online.org">phenixbb-leave@phenix-online.org</a>&gt;<br>
<div class="HOEnZb"><div class="h5">&gt;<br>
&gt;<br>
&gt;<br>
&gt;<br>
&gt; _______________________________________________<br>
&gt; phenixbb mailing list<br>
&gt; <a href="mailto:phenixbb@phenix-online.org">phenixbb@phenix-online.org</a><br>
&gt; <a href="http://phenix-online.org/mailman/listinfo/phenixbb" rel="noreferrer" target="_blank">http://phenix-online.org/mailman/listinfo/phenixbb</a><br>
&gt; Unsubscribe: <a href="mailto:phenixbb-leave@phenix-online.org">phenixbb-leave@phenix-online.org</a><br>
&gt;<br>
<br>
</div></div><br>_______________________________________________<br>
phenixbb mailing list<br>
<a href="mailto:phenixbb@phenix-online.org">phenixbb@phenix-online.org</a><br>
<a href="http://phenix-online.org/mailman/listinfo/phenixbb" rel="noreferrer" target="_blank">http://phenix-online.org/mailman/listinfo/phenixbb</a><br>
Unsubscribe: <a href="mailto:phenixbb-leave@phenix-online.org">phenixbb-leave@phenix-online.org</a><br></blockquote></div><br></div>